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ABSTRACT: 

A  full  account  is  given  of  a  test  program  to  evaluate  a  prototype 
"Hybrid"  compressor,  wherein  the  flow  is  turned  through  180  (from  axial 
to  radial  to  axial)  within  the  rotor.   The  initial  stator,  or  diffuser, 
consisted  of  two  cascades  of  unswept  and  untwisted  blades  arranged  and 
designed  as  in  an  axial  compressor.  The  performance  was  measured  with 
two  diff users  of  different  solidities  and  with  changes  in  the  blade  stagger 
angles.  The  maximum  efficiency  approached  80$;  however,  the  flow  was 
separated  in  the  second  blade  row  for  all  diffuser  designs  tested.  Surveys 
from  hub  to  tip  at  the  rotor  exit  revealed  an  unusually  large  variation 
in  flow  angle.  Approximations  were  obtained  from  normalised  distributions 
of  velocity  and  angle  from  all  tests  which  will  allow  the  design  of  a 
diffuser  with  swept  or  twisted  blades.   It  is  recommended  that  such  a 
design  be  tested  to  provide  a  more  definitive  assessment  of  the  "Hybrid" 
concept. 

The  reported  tests  were  the  first  conducted  in  a  new  turbine -driven 
compressor  test  rig  at  the  Turbopropulsion  Laboratories.  An  account  is 
included  of  the  verification  of  flow,  torque  and  temperature  measurements 
and  of  a  method  devised  to  overcome  the  difficulty  of  establishing  the 
performance  with  highly  non-uniform  total  pressure  distributions  from 
multiple  passages  within  a  thin  annulus.  The  method,  which  was  analyti- 
cally correct  and  required  the  input  of  a  particular  number  of  sample 
measurements,  gave  results  almost  identical  to  the  conventional  mass- 
averaging  technique. 

This  research  program  was  carried  out  for  the  Department  of  the  Navy, 
Naval  Air  Systems  Command,  in  accordance  with  Airtask  No.  3303300/186B/ 
3F*+lU 32301  under  the  cognizance  of  Mr.  K.  H.  Guttmann,  Code  ATR-330C, 
NAVAIR  Project  Engineer. 
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1.  Introduction 

The  "Hybrid"  compressor  consists  of  a  centrifugal  rotor  in  which  the 
flow  is  turned  through  180  ,  from  axial  out  to  radial  and  back  to  axial, 
followed  by  an  axial  stator  or  diffuser.   This  arrangement  was  patented  by 
Mr.  George  Derderian  of  the  Naval  Air  Systems  Command  (U.S.  Patent  3)365,982 
of  January  30,  1968).   If  such  a  compressor  could  be  built  for  high  pressure 
ratios  at  good  efficiencies,  it  would  be  of  great  advantage  in  simple  gas 
turbines  for  many  applications. 

The  Naval  Air  Systems  Command  asked  Dr.  M.  H.  Vavra  of  the  Department 
of  Aeronautics  at  the  Naval  Postgraduate  School  to  evaluate  this  new  concept. 
Dr.  Vavra  pointed  out  that  the  flow  from  a  Hybrid  rotor  could  not  be  pre- 
dicted analytically,  and  in  practice  might  be  so  irregular  that  good  pressure 
recovery  in  the  axial  diffuser,  and  hence  good  overall  efficiency,  could  not 
be  achieved.  Subsequently  a  program  was  initiated  by  Naval  Air  Systems 
Command  to  test  the  feasibility  of  the  proposed  arrangement  at  the  Turbo- 
Propulsion  Laboratory  at  the  Naval  Postgraduate  School.  A  prototype  Hybrid 
compressor  was  designed  by  Dr.  M.  H.  Vavra  to  be  compatible  with  a  new  test 
rig  also  designed  by  him  and  then  under  construction.  A  complete  descrip- 
tion of  the  aerodynamic  and  mechanical  design  of  the  test  rig  and  the 
calculated  performance  of  the  Hybrid  compressor  are  given  in  Reference  1. 
The  compressor  was  built  and  was  installed  at  the  Turbo-Propulsion  Labora- 
tory during  1972 .  A  description  of  this  facility  is  given  in  Reference  2 . 
Testing  began  in  December  1971  when  the  author  joined  Dr.  Vavra  at  the 
Naval  Postgraduate  School  on  a  temporary  leave  of  absence  from  the  Boeing 
Company,  and  was  completed  while  remaining  in  the  Department  of  Aeronautics 
as  Naval  Air  Systems  Command  Visiting  Professor.   The  course  of  the  test 
program  was  guided  by  Dr.  Vavra  and  was  completed  in  May  1973.   The  initial 
tests  of  the  Hybrid  compressor  were  also  the  first  tests  conducted  with 


the  test  rig  itself.   The  development  of  on-line  reduction  and  the  verifi- 
cation of  primary  measurements  were  a  necessary  part  of  the  test  program. 
The  present  report  gives  a  detailed  account  of  the  complete  test  program 
and  presents  the  results  obtained  with  the  prototype  Hybrid  compressor. 

In  reporting  the  test  program  a  division  has  been  made  between  two 
phases  of  the  work.   In  the  'Preliminary  Test  Program'  reported  in  Section 
3,  changes  were  made  in  the  arrangement  of  components  and  in  instrumentation 
as  progressively  more  knowledge  of  the  machine  was  gained.   The  measurement 
of  a  large  variation  in  the  flow  angle  of  the  rotor  exit  required  that  the 
separation  between  the  rotor  and  stator  be  increased.  The  efficiency  of 
the  rotor  and  stator  alone  was  then  more  difficult  to  establish  because 
the  space  available  for  pressure  measurements  was  limited  and  the  flow  pro- 
files there  were  highly  distorted.   Consideration  was  then  given  to  the 
problem  of  establishing  efficiency  from  a  limited  number  of  fixed  impact 
pressure  probe  measurements,  and  this  is  described  in  Appendix  F.  Many  of 
the  details  involved  in  verifying  the  measurements  made  in  the  new  test  rig 
are  included  for  completeness  in  other  Appendices. 

The  second  phase  was  to  measure  and  analyze  the  performance  of  the 
compressor  for  two  different  diffusers  with  changes  in  diffuser  blade  angles. 
During  these  tests,  which  are  described  in  Section  U,  the  instrumentation 
and  methods  applied  were  not  changed.   Included  in  these  tests  were  surveys 
of  the  flow  at  the  exit  of  the  rotor.   In  Section  5  analytical  approximations 
to  a  large  number  of  data  are  obtained  which,  it  is  recommended  in  Section 
6,  should  be  used  to  design  diffusers  specifically  for  the  unique  flow 
generated  by  the  Hybrid  compressor  rotor. 


2.   Description  of  the  Test  Installation 

2.1  Hybrid  Compressor 

Figure  1  shows  the  construction  of  the  compressor.   The  design  is 
described  in  detail  in  reference  1,  and  only  the  features  of  significance 
in  the  test  program  will  be  mentioned  here.   The  arrangement  shown  in 
Figure  1  was  not  changed  throughout  the  preliminary  testing  which  is 
described  in  Section  3.   The  only  change  in  geometry  was  to  the  leading 
edges  of  the  inducer,  as  is  shown  in  Figure  2.   During  the  later  tests 
described  in  Section  h,   changes  were  made  in  the  axial  positions  of  the 
stator  blade  rows  and  in  the  blade  angles.   These  changes  are  given  in 
Section  k.l.     A  view  of  the  rotor  is  given  in  Figure  3?  and  other  views 
of  the  blade  passages  are  introduced  in  Section  U.5. 

The  flow  enters  the  rotor  axially.  A  honeycomb  section  ensures  that 
the  flow  leaves  the  compressor  with  zero  velocity  in  the  circumferential 
direction.   Thus  the  torque  supplied  to  the  rotor  is  equal  to  the  torque 
experienced  by  the  blading  and  honeycomb  of  the  stator  section.   The  stator 
section  is  mounted  on  ball  races  and  is  free  to  rotate  against  flexures 
instrumented  with  strain  gauges.   The  strain  gauge  reading  is  then  a  mea- 
sure of  the  torque  supplied  to  the  rotor.   Together  with  the  rotational 
speed  measured  by  a  magnetic  flux  cutter,  a  measure  of  the  power  supplied 
to  the  compressor  is  obtained. 

2.2  Compressor  Test  Rig 

The  compressor  test  rig  was  designed  initially  to  accomodate  a  single- 
stage  axial  transonic  compressor  at  speeds  up  to  31*500  r.p.m.   The  Hybrid 
compressor  was  subsequently  designed  to  adapt  to  the  existing  rig,  and 
much  of  the  difficulty  experienced  in  the  test  program  was  for  this  reason. 

The  arrangement  of  the  test  rig  (which  is  described  in  reference  l)  is 
shown  in  Figure  k.      The  test  compressor  was  driven  by  an  opposed-rotor 


single  stage  air  turbine  supplied  from  the  laboratory  compressed  air 
power  supply.   (A  12-stage  All is -Chalmers  axial  compressor  supplies  up  to 
11  lbs.  per  sec.  at  up  to  2  atmospheres  gauge  pressure.   These  facilities 
are  described  in  reference  2.)  The  power  to  the  drive  turbine  was  obtained 
from  the  flow  rate  metered  by  an  orifice  and  the  difference  between  the 
air  temperature  measured  into  and  out  of  the  turbine. 

The  Hybrid  compressor  rotor  was  connected  by  a  splined  shaft  to  the 
rotor  of  the  turbine  drive  unit.  Air  was  brought  into  the  compressor 
through  a  rotating-plate  throttle  followed  by  a  settling  chamber  containing 
eight  perforated  plates.  A  section,  3  inches  thick,  of  l/U-inch  aluminum 
honeycomb  restrained  by  a  #8  gauge  wire  mesh  screen  was  installed  in  the 
chamber  downstream  of  the  plates  during  the  preliminary  tests.  A  flow 
nozzle  was  contained  in  the  20  ft.  long  induction  pipe  between  the  settling 
chamber  and  the  compressor.   The  flow  nozzle  was  moved  upstream  during  the 
preliminary  tests  from  its  original  position,  and  the  inlet  to  the  compressor 
was  contoured  as  shown  in  Figure  5. 

The  measurement  of  flow  rate  given  by  the  flow  nozzle  and  the  difference 
in  air  temperature  measured  into  and  out  of  the  compressor  provided  a  second 
measure  of  the  power  absorbed. 

A  view  of  the  compressor  test  rig  is  shown  in  Figure  6  and  of  the 
Hybrid  compressor  installed  for  testing  in  Figure  7. 
2.3  Instrumentation  and  Data  Recording 

During  a  compressor  test,  operating  data  was  monitored  by  an  operator 
at  a  control  console.  Records  were  kept  of  oil  drop  rate  and  pressure  to 
the  bearing  lubrication  systems ,  of  temperatures  of  the  turbine  and 
compressor  bearings,  and  vibration  patterns  from  four  accelerometers  on 


the  turbine  drive  unit.  This  data  was  mainly  of  interest  for  the  safety 
of  operation  and  will  not  be  reported. 

Performance  data  included  pressure  probe  and  temperature  probe 
measurements,  wall  pressure  measurements  and  torque  transducer  measure- 
ments. Figure  8  shows  the  position  of  and  notation  for  the  temperature 
and  pressure  sensors.  A  data  logging  system  (described  in  reference  2) 
was  used  to  scan  and  to  record  performance  data.  A  schematic  of  the 
system  is  shown  in  Figure  9-   The  output  from  the  data  scan  was  a  punched - 
tape  and  teletype  listing  of  channel  numbers  and  data  value.   Preliminary 
data  reduction  was  carried  out  immediately  by  inputting  the  tape  into  a 
second  teletype  unit  which  was  audio-coupled  to  the  Naval  Postgraduate 
School  IBM  360  computer.  The  on-line  data  reduction  program  printed  out 
turbine  drive  and  test  compressor  performance  parameters  as  the  test 
progressed,  and  also  caused  the  primary  data  to  be  punched  onto  cards 
which  were  submitted  later  for  batch  processing. 

Pressures  were  exhibited  on  mercury  and  water  manometers.  Selected 
pressures  were  scanned  by  two  2k   port  Scanivalves  onto  Statham  absolute 
pressure  transducers.  Separate  differential  transducers  recorded  the 
turbine  drive  flow  orifice  and  the  compressor  flow  nozzle  differential 
water  columns.  An  absolute  transducer  recorded  the  pressure  upstream  of 
the  turbine  drive  flow  orifice.  Following  an  initial  pneumatic  calibra- 
tion of  the  transducers,  their  output  was  adjusted  before  each  test  by 
applying  known  electrical  calibration  signals  to  each  transducer  bridge. 
During  operation,  the  calibrations  are  routinely  reexamined  by  recording 
water  and  mercury  manometer  readings  at  intervals,  and  where  necessary  a 
correction  was  introduced  into  the  batch  reduction  program. 


Outputs  from  iron-constantan  thermocouple  probes  were  recorded  in 
millivolts  "by  the  scanning  system.  A  Pace  Engineering  Model  BRJW18S-12TP- 
80776  reference  junction  was  incorporated,  and  the  temperatures  were 
evaluated  using  an  analytical  representation  of  the  ISA  Standard  Tables, 
identified  as  Function  Subprogram  Tempic  in  Appendix  H.  Both  United 
Sensor  Corporation  and  in-house  temperature  probes  were  used,  and  where 
necessary  they  will  be  identified  in  later  sections. 

Torque  was  initially  measured  by  strain  gauges  on  flexures  between 
the  stator  and  the  main  frame  of  the  compressor.  As  reported  in  Section 
3,  a  revised  system  using  a  single  beam  flexure  was  built  external  to 
the  machine  (see  Figure  7).  As  is  shown  in  Appendix  C,  by  calibrating 
before  and  immediately  following  a  test  by  hanging  weights,  an  accuracy 
of  -  0.1$  in  the  torque  measurement  was  expected.   The  axial  force  on  the 
rotor,  kept  within  convenient  limits  by  a  balance  piston,  was  also  measured 
by  a  flexure  system.   These  measurements  are  not  reported  however  since 
the  axial  force  was  not  introduced  into  the  performance  analysis,  and 
since  the  calibration  was  not  repeated  following  the  initial  installation. 
2.k.      Data  Reduction 

A  computer  program,  designated  TCRIG1,  was  written  to  accept  the 
primary  data  input  from  punched  paper  tape  and  to  output  performance 
parameters  on  the  teletype  unit  in  the  laboratory.   The  program  also 
outputs  the  primary  data  on  cards  for  permanent  storage  and  reuse.   The 
final  version  of  the  program,  which  was  not  changed  significantly  during 
the  Hybrid  testing,  is  listed  in  Appendix  H. 

Based  on  TCRIG1,  but  continually  revised  during  the  testing  program, 
a  batch-processing  computer  program  was  also  written.   The  final  version 


of  this  program,  designated  TCRIG-U,  is  also  listed  in  Appendix  H.   The 
batch  program  accepted  data  cards  from  up  to  50  data  scans,  allowed 
changes  to  calibration  slopes  and  intercepts  to  be  introduced,  and  out- 
put reduced  data  on  cards  which  were  resubmitted  with  previous  test 
results  for  computer  plotting.   However,  the  major  difference  between 
TCRIG-U  and  TCRIG-1  was  the  inclusion  in  the  batch  program  of  the  calcu- 
lation of  machine  performance  according  to  the  methods  described  in 
Appendix  F. 

In  the  description  which  follows,  the  final  form  of  the  data 
analysis  (Program  TCRIG-U)  is  described.  The  results  given  in  Section  k 
were  reduced  using  this  program.  The  results  from  the  preliminary  test 
program  reported  in  Section  3  were  obtained  with  fewer  pressure  and 
temperature  probes  downstream  of  the  rotor,  and  evaluation  of  the 
performance  was  not  as  complete. 

The  weight  flow  to  the  turbine  drive  and  the  horsepower  delivered 
by  the  turbine  were  calculated  from  the  following  equations: 


r                 "i  r    °-027  ^T-i  /  pt  '  hwt 
WJ  =  6.7720  [l  +  0.000015  (Tt  -  60)  [l f — j  J  Jff   +  \2o)    C1) 
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HPT  =  1A137  Cp  WT  (TtT1  -  TtT2)  (3) 


Eq.  (l)  and  Eq.  (2)  were  derived  in  reference  3  and  describe  the  calibra- 
tion of  the  flow  orifice  using  vena  contracta  pressure  taps. 


The  first  measure  of  the  horsepower  to  the  compressor  was  obtained 
from  the  measurements  of  torque,  M  ,  and  speed,  N  :   thus 

HPC1  =  1.5865  x  10"5  .M.N.  (k) 

The  second  measure  was  obtained  from  the  flow  rate,  W   ,  and  temperature 
rise: 

HPC2  =  1.4137  Cp  W  (Ttl)  -  Tt0)  (5) 

The  weight  flow  to  the  compressor  was  given  by  Eq.B(7)in  Appendix  B  which 
resulted  from  a  calibration  of  the  flow  nozzle. 

The  performance  of  the  compressor  is  described  in  terms  of  referred 
quantities  that  retain  their  original  units: 

WREF  =  *  ^~/6 


NREF  =  *M   9 

_  (6) 

HPEEF  =  HP/5  ^~ 

^F  "  M/6 

where  9  =  T^A^p  and  6  =  Pq/Pj^j,   ,  and  where  T^  =  517. 8°R  and 
P^  -  29.92  ins.  Hg. 

The  performance  to  be  evaluated,  was  that  of  the  rotor  and  stator 
combination.   The  honeycomb  flow  straightener  in  the  outlet  duct 
(Figure  l)  was  necessary  in  order  for  the  torque  to  be  measured,  but  the 
losses  occuring  in  this  section  should  not  be  included  in  appraising  the 
compressor  performance.   Equally,  the  annular  90     bend  in  the  outlet 


flow  was  dictated  by  the  constraints  of  the  test  rig  and  losses  occurring 
here  should  also  be  excluded.   The  total -to-total  pressure  ratio  and 
efficiency  from  the  inlet  (station  0)  to  the  exit  of  the  stator  (station  k) 
were  therefore  the  performance  parameters  of  greatest  interest . 

The  compression  process  across-  rotor  and  stator  is  represented  in 
Figure  10.  The  total -to-total  efficiency  is  defined  as  the  ratio  of  the 
power  required  in  an  isentropic  compression  to  the  same  outlet  pressure, 
to  the  actual  power  required.  The  actual  power  is  given  by  Eq.  (k) ,  using 
measurements  of  torque  and  speed,  or  by  Eq.  (5)j  using  measurements  of  flow 
rate  and  temperatures.  The  ideal  horsepower,   IHP   ,  is  given  by 


-  1.1*137  ic[(lj   -T  ]  (7) 
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where   T   )    ,  in  Figure  10,  is  related  to  the  pressures  in  an 
V  Vis 


isentropic  process  by 
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which,  in  Eq.  (7)?  gives 
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IHPc  =  l.Ul37WCpT       [(-^)     "     -l]  (9) 


If  the  flow  were  uniform  at  station  (0)  and  at  station  (h)   and  if  all 
measurements  could  be  made  without  error,  a  single  evaluation  of  the 
efficiency  could  be  made  using  either  Eq.  (k)   or  Eq.  (5)  with  Eq.  (9). 
Because  it  was  found  in  the  present  tests  that  the  flow  was  significantly 
non-uniform  and  that  some  measurements  (torque,  for  example)  could  be 


made  with  greater  accuracy  than  others  (temperature  rise  for  example),  a 
total  of  six  different  evaluations  of  the  total -to-total  efficiency  were 
examined  in  the  course  of  the  test  program.  The  significance  of  the 
different  efficiencies  will  be  discussed  in  Sections  3  and  k.      The 
analytical  expressions  and  use  of  the  measurements  are  described  here. 

(i)  \ 

From  Eq.  (5)  and  Eq.  (9), 


Yzl 


T    was  measured  as  the  average  of  two  probes  ahead  of  the 
inlet  to  the  compressor.  T    was  measured,  as  described  in 
Appendix  D,  as  the  average  of  the  readings  of  four  selected 
probes  inserted  axially  through  the  struts  supporting  the 
honeycomb  section. 

P   was  measured  by  a  single  Kiel  probe  in  the  compressor 
inlet,  p    was  taken  as  the  mass  average  of  the  readings  of 
12  Kiel  probes  at  the  stator  exit.   (in  the  preliminary  test 
program,  the  average  of  only  four  impact  probes  was  used. ) 

(ii)  \ 

A  non-uniform  profile  at  the  compressor  inlet  was  found 
with  the  earliest  inlet  arrangement.   The  inlet  stagnation 
pressure  was  then  calculated  from  measurements  of  mass  flow 
and  static  pressure  using  relationships  which  are  derived  in 
Appendix  A. 
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and  Eq.  A(12),  it  follows  that 
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where  Xn  =  V  /V    ,  V   =  V2C  T    ,  and  B   is  defined 

0  t0      t0      P  0 
by  Eq.  A (13)  applied  at  station  0.  X   is  obtained  as  a 

solution  of  Eq.  A(l2),  or 
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(iii)     T|3 


evaluates  an  efficiency  which  involves  the  flow  rate  measure- 
ment, whereas  J\       does  not. 

From  Eq.  (k)    and  Eq.  (9), 
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so  that  j]       depends  on  the  measurement  of  torque  and  flow- 
rate  but,  unlike  T]   ,  does  not  depend  on  the  measurement 
of  temperature  rise.   In  the  results  reported  in  Section  U, 
the  mass  averaged  impact  pressure  at  the  exit,  given  by 
Eq.  F(3),  was  used  in  Eq.  (15). 

Analogously  to  the  use  of  the  inlet  static  pressure 
and  mass  flow  rate  in  the  definition  of  j\        ,  the  exit 
static  pressure  and  mass  flow  rate  could  be  used  instead  of 
measurements  of  exit  impact  pressure.  The  equations  are 
derived  in  Appendix  F2  where  Tk   is  given  by  Eq.  F(ll). 
However,  in  order  to  evaluate  7],   ,  a  value  must  be  assumed 
for  the  blockage  factor,  C,   ,  which  is  defined  in  Eq.  F(7) 
and  appears  in  the  continuity  equation,  Eq.  F(6).   Unlike  in 
the  inlet,  where  the  flow  was  shown  to  be  very  nearly  uniform, 
and  a  blockage  factor  of  near  unity  was  realistic  the  exit 
blockage  factor  was  totally  unknown,  and  would  be  expected 
to  vary  over  the  throttling  range.   Therefore  the  definition 
of  T|i   was  used  only  as  a  means  of  determining  a  value  for 
the  exit  blockage  factor.  Equating, 

\  =  h  (16) 

and  using  Eq.  (15)  and  Eq.  F(ll)  to  obtain  the  product 

2B.X.   ,  Eq.  F(8)  was  solved  for  X,   and  C,   was  then  given 

by  Eq.  F(9)- 
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The  above  four  efficiencies  axe  based  on  a  definition 
of  efficiency  for  uniform  flow.  Four  fixed  impact  probes 
at  the  stator  exit  showed  the  velocity  profile  to  be  highly 
non-uniform.  As  shown  in  Section  3>  different  methods  of 
averaging  the  measurements  gave  very  different  values  of 
efficiency.  By  assuming  an  analytical  form  for  the  distri- 
bution of  the  impact  pressure ,  involving  constants 
determined  from  the  available  measurements,  the  efficiency 
to  the  non -uniform  exit  condition  could  be  determined  by 
integration.  This  avoided  the  ad  hoc  definition  of  averages, 
for  the  use  of  which  there  is  no  analytical  justification. 

The  efficiency  T)e  w&s  obtained  using  the  distribution 
5 

of  impact  pressure  shown  in  Figure  F2.  The  result,  given  as 
Eq.  F(30)  in  Appendix  F3.2  is  repeated  here: 

IHPC  (Eq.  F(2l0) 


h   =      HPC1  (Eq.  (k))  <17> 


and  the  associated  blockage  factor  C   is  given  by  Eq.  F(25). 

The  evaluation  of  J\       and  C   depended  wholly  on  only 
four  impact  probe  measurements  on  the  center  line  of  the  exit 
annulus .  The  evaluation  was  continued  after  the  following 
improved  method  was  devised  in  order  that  a  comparison  could 
be  made  with  earlier  performance  measurements. 

(vi)  \ 

The  calculation  of  the  final  efficiency  is  described  in 

Appendix  F3-3-  A  distribution  of  impact  pressure  was  assumed 
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radially  between  the  hub  and  tip  in  each  of  six  segments  of 
a  blade  passage.  Each  distribution  was  determined  by  two 
Kiel  probe  measurements  at  the  corresponding  peripheral 
position.   The  total  of  twelve  Kiel  probes  however  were 
positioned  in  separate  blade  passages  chosen  at  random 
around  the  annulus.   In  this  way  variations  between  blade 
passages  were  partially  accounted  for,  and  the  distribution 
shown  in  Figure  F3  would  represent  an  average  for  a  single 
blade  passage. 

The  shape  of  the  assumed  radial  distribution  was  found 
to  be  realistic  in  probe  surveys  reported  in  Section  3.   -A 
reversed  flow  is  allowed  near  the  wall,  occur ing  when  the 
two  measurements  that  specify  the  distribution  are  widely 
different.   Continuity  is  satisfied,  and  serves  to  deter- 
mine a  blockage  factor  in  the  peripheral  direction. 

Thus  the  efficiency  given  in  Eq.  F(3*0  and  repeated 
here, 


IHPQ  (Eq.  F(33)) 
^6  =   HPC1  (Eq.  (10)  (l8) 


was  considered  to  be  a  good  measure  of  the  performance  of  the 
compressor  up  to  station  h.      Also,  the  pressure  ratio  given 
by  Eq.  F(38)>  which  is  defined  by  equating  the  integrated 
ideal  horsepower  (Eq.  F(33))  "to  the  expression  for  uniform 
flow  (Eq.  (9)) ?  was  considered  to  be  a  reasonable  basis  of 
comparison  when  evaluating  the  compressor  performance  with 
different  stator  geometries. 


1U 


The  calculation  of  an  efficiency  to  a  "fully  mixed" 
condition  was  considered,  "but  rejected.  This  would  involve 
a  constant  pressure  mixing  process  from  actual  conditions 
at  station  k   to  an  assumed  uniform  axial  flow  at,  say, 
kl     However,  the  unknown  flow  angle  at  station  k   makes  the 
calculation  impossible  without  making  assumptions  that  can 
not  be  supported.   In  the  above  methods,  no  attempt  was  made 
to  separate  the  area  blockage  factor  from  the  cosine  of  the 
flow  angle,  the  product  of  which  was  measured  with  reasonable 
assumptions.   However,  it  was  noted  that  CL-  was  calculated 
to  have  values  about  equal  to  the  cosine  of  the  expected 
average  flow  angle.   From  Eq.  F(37)?  the  area  blockage  factor 
assumed  in  specifying  the  impact  pressure  distribution, 
Cs      ,  was  therefore  acceptably  close  to  unity. 

The  overall  static -to-static  efficiency,  J\        ,  of  the 
compressor  was  given  by,  referring  to  Figure  10, 

no  -  ^jA  (19) 

Using  Eq.  (5),  Eq.  A(5),  Eq.  (15),  Eq.  (k)   and  the  relation- 
ship between  temperature  and  pressure  for  isentropic 
compression,  Eq.  (19)  may  be  written  in  terms  of  quantities 
already  evaluated  as 
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Since  T\       was  evaluated  using  Eq.  (15)  >  the  evaluation 
of  7]   does  not  require  a  knowledge  of  the  exit  flow  profile 

The  performance  of  the  diffuser  was  described  by  the 
pressure  recovery  factor,  R^   ,  defined  as 


Rn  --JL-1  (21) 

D   P^-Pg 


and  the  efficiency,      T]        ,   defined  as 


T+    -T, 
h 
^   =  m      -7m    \  (22) 


'4 


t,.    v   Vis 


Equation    (22)    can  be  written  as 
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In  calculating  R  and  7]  ,  the  static  pressures  were 
taken  to  be  equal  to  the  recorded  wall  pressures,  the  stagna- 
tion pressure  at  the  rotor  exit,  p    was  measured  on  the 

t2 
annulus  center line  by  a  cobra  probe,  and  p    was  obtained 

\ 
from  Eq.  P(38). 
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3.   Preliminary  Test  Program 

3.1.  Introduction 

The  first  tests  of  the  Hybrid  compressor  were  concurrently  the  first 
tests  of  the  compressor  test  rig  described  in  Section  2.  The  geometry  of 
the  stator  was  at  first  held  fixed  while  each  of  the  measurements  required 
to  establish  the  compressor  performance  was  verified.  Redundancy  in  the 
measures  of  actual  horsepower  resulted  in  the  detection  and  elimination  of 
discrepancies  in  the  individual  measurements .  Bearing  failure  was 
experienced  on  two  occasions,  but  an  error  in  the  installation  was  found 
and  corrected.  The  poor  indicated  performance  was  improved  by  an  axial 
movement  of  the  stator  rows  away  from  the  rotor  exit.  However,  because  of 
the  reduced  space ,  it  was  then  more  difficult  to  establish  the  average 
properties  of  the  flow  at  the  stator  exit.  Additional  pressure  probes 
were  installed  and  their  measurements  reduced  using  an  analysis  for  non- 
uniform exit  flow  which  is  given  in  Appendix  F. 

Infection  33  the  results  of  the  early  tests  which  are  of  interest, 
are  reported.  The  instrumentation  developed  during  these  tests  was  then 
held  fixed  during  the  tests  reported  in  Section  k. 

3.2.  Effect  of  Contouring  the  Inducer  Leading  Edges 

The  performance  of  the  compressor  was  measured  at  8,000  and  10,000 
r.p.m.  The  leading  edges  of  the  inducer  were  as  shown  in  Figure  2(a), 
the  intake  was  that  shown  in  Figure  5(a),  and  the  performance  was  based  on 
the  manifold  average  of  four  impact  pressures  measured  on  the  annulus 
centerline  at  station  k.     The  indicated  performance  was  poor  and  the 
inducer  leading  edges,  which  were  not  initially  contoured  in  order  to  test 
arguments  given  in  reference  1,  were  reworked  to  have  the  contour  shown  in 
Figure  2  (b  ). 
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The  performances  with  square  and  contoured  inducer  leading  edge 
geometries  are  compared  in  Figure  11.   There  was  an  improvement  in  the 
pressure  ratio  and  the  efficiency  at  a  given  flow  rate,  however  the 
maximum  efficiency  increased  by  only  1%   from  the  initial  low  value  of 
68.5$>.   The  accuracy  of  each  of  the  measurements  was  questioned  and  the 
torque  measurement  of  horsepower  (which  is  not  involved  in  T| ,  )  was  found 
not  to  agree  with  the  measurements  based  on  flow  rate  and  temperature 
rise.  The  torque  measuring  system  was  later  revised  as  described  in 
Appendix  C  and  an  explanation  was  sought  for  the  low  value  of  T]  n  • 
3.3«   Measurements  of  the  Flow  Into  and  Out  of  the  Rotor 

The  position  of  the  flow  nozzle  in  the  initial  tests  (Figure  5(a)), 
was  the  result  of  adapting  the  Hybrid  compressor  to  an  existing  test  rig. 
A  Kiel  probe  survey  was  made  at  station  0  and  showed  the  impact  pressure 
distribution  (Figure  12)  to  be  unacceptably  non-uniform.   The  flow  nozzle 
was  relocated  as  shown  in  Figure  5(^)  and  the  interim  arrangement  shown 
in  Figure  5(b)  was  used  for  several  tests.  The  entrance  to  the  compressor 
was  faired  with  red  wax.  As  can  be  seen  in  Figure  12,  both  the  interim 
and  final  arrangements  gave  nearly  uniform  flows  into  the  compressor.  The 
slight  asymmetry  is  the  effect  of  probe  blockage. 

Flow  angle  measurements  were  made  at  the  exit  of  the  rotor  using  a 
l/8   inch  5-hole  United  Sensor  probe  for  the  final  intake  arrangement. 
The  results  shown  in  Figure  13  are  for  a  condition  close  to  surge  at 
10,000  r.p.m.  The  curves  for  velocity  were  derived  from  the  measured  flow 
angles,  impact  pressure,  exit  static  pressure  (which  was  assumed  to  be 
uniform)  and  the  measured  stagnation  temperature.   The  angle  of  the  flow 
with  respect  to  the  axial  direction  measured  in  a  circumferential 
plane   (o  )   is  seen  to  vary  as  much  as  15  over  the  exit  blade  height. 
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When  the  flow  angle  is  considered  with  the  absolute  velocity,  the  profile 
of  the  calculated  time-averaged  flow  velocity  relative  to  the  rotor  is 
seen  to  be  peaked  toward  the  tip. 

The  flow  relative  to  the  rotor  in  this  unconventional  design  was  of 
particular  interest  since  the  design  of  the  diffuser  assumed,  in  the 
usual  way,  that  the  incident  flow  would  be  uniform.  A  knowledge  of  the 
actual  flow  at  the  rotor  exit  might  lead  to  a  better  diffuser  design. 
The  observation  on  dissassembly  of  particulate  deposits  in  the  pressure- 
side  outer  corners  of  the  rotor  suggested  that  some  separation  had  occured 
in  the  flow  in  the  rotor.   It  was  evident  however,  from  the  measured  flow 
angle  variation,  that  the  existing  diffuser  with  straight  blades  might  not 
be  suitable  for  this  type  of  compressor,  and  might  be  the  reason  for  the 
low  efficiency.   In  order  to  reduce  the  effect  of  the  non-uniform  flow 
from  the  rotor  on  the  performance  of  the  diffuser,  the  diffuser  blade  rows 
were  moved  away  from  the  rotor  as  far  as  the  physical  constraints  allowed 
(Figure  l).  This  change  was  preceded  by  tests  to  establish  a  preliminary 
performance  map  and  tests  to  verify  the  accuracy  of  the  measurements  of 
torque  and  temperature. 
3.^.   Performance  with  Small  Rotor -St ator  Clearance 

Performance  data  were  measured  at  10,000,  12,000  and  lU,000  r.p.m. 
The  leading  edges  of  the  rotor  were  contoured  as  in  Figure  2(b)  and  the 
intake  configuration  was  that  shown  in  Figure  5(c).  The  torque  was 
measured  with  the  original  flexures,  and  the  stagnation  pressure  at 
station  k   was  measured  as  the  manifold  (area)  average  of  four  probes  at 
the  face  of  the  honeycomb  flow  straightener . 

Figure  Ik   shows  the  efficiency  TU  ,  and  the  stagnation  pressure 
ratio  as  a  function  of  the  referred  flow  rate.  At  10,000  r.p.m.  there 
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was  a  3  1/2%  improvement  in  the  peak  efficiency  over  the  results  shown  in 
Figure  11  for  the  initial  intake  arrangement.   The  pressure  ratio  is 
everywhere  increased.   There  is  seen  to  he  no  measurable  effect  of 
compressor  speed  on  the  peak  efficiency  of  72.5%. 

TL   ,  from  Eq.  (10) ,  is  shown  rather  than  f]   ,  from  Eq.  (15), 
because  the  torque  measurement  was  unreliable  during  the  early  tests. 
This  is  seen  in  Figure  15  which  compares  three  separate  measures  of  the 
machine  horsepower  as  a  function  of  compressor  flow  rate  at  12,000  and 
1^,000  r.p.m.  HP    and  HP   which  depend  on  the  measurements  of  flow 
rate  through  and  temperature  change  across  the  compressor  and  drive 
turbine  respectively  (Eq.  (5)  and  Eq.  (3))?  are  in  very  good  agreement. 
HP-,,  which  depends  on  the  measurement  of  speed  and  torque  (Eq.  (*0),  is 
seen  to  disagree  by  as  much  as  25%,  and  its  use  in  Eq.  (ik)   gave  values 
of  efficiency  which  were  greater  than  unity. 
3.5«   Performance  with  Maximum  Rotor -St ator  Clearance 

The  blade  rows  of  the  stator  were  moved  away  from  the  rotor  by  a 
rearrangement  of  the  spacing  rings  shown  in  Figure  1.   The  increased 
separation,  which  placed  the  trailing  edges  of  the  second  row  of  blades  of 
Diffuser  A  at  the  exit  of  the  0.6l"  annular  duct,  was  intended  to  allow  as 
much  mixing  as  possible  to  occur  in  the  non-uniform  flow  leaving  the  rotor 
(Figure  13).   However,  the  performance  measured  for  this  configuration 
cannot  be  directly  compared  with  the  performance  described  in  Section  3«*+« 
The  axial  movement  of  the  blade  rows  placed  the  station  k   pressure  and 
temperature  sensors,  which  were  mounted  through  struts  supporting  the  exit 
honeycomb  (Figure  8),  very  close   (~  1/V)  to  the  exit  plane  of  the 
stator.  Also,  following  an  examination  of  the  temperature  measurements 
carried  out  during  the  preceding  tests,  eight  shielded  thermocouple  probes 
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were  installed  in  the  struts  to  replace  four  unshielded  probes .  An 
account  of  the  difficulties  found  in  the  measurement  of  temperatures  is 
given  in  Appendix  D.  The  measurement  of  torque  was  completely  revised  as 
described  in  Appendix  C  and  an  accuracy  of  0.1$  was  obtained.   The 
clearances  of  the  labyrinth  glands  in  the  stator  section  were  increased 
to  eliminate  interference  in  the  torque  measurement  from  differential 
heating  of  the  machine  supports.  The  effect  of  labyrinth  leakage  on  the 
efficiencies  however  was  analyzed  and  is  shown  in  Appendix  E  to  have  been 
small. 

Performance  measurements  at  four  speeds  are  shown  in  Figure  l6  and 
Figure  17.   In  Figure  16,  the  referred  torque  is  shown  as  a  function  of 
referred  flow  rate.   The  consistency  of  the  torque  measurement  is  evident 
in  these  results.  The  torque  varies  linearly  with  the  flow  rate  at  a 
given  speed  and  linearly  with  the  speed  at  a  given  flow  rate.  Both  the 
efficiencies  T]j  ,  and  T|   are  shown  in  Figure  17.  For  the  data  shown 
from  RUN  36  and  RUN  37?  the  exit  temperature  was  evaluated  from  the 
manifold  average  of  seven  thermocouple  probes  (one  of  the  eight  was  found 
to  be  broken) .  Superposed  are  the  data  from  RUN  39  in  which  the  outputs 
of  the  seven  good  probes  were  separately  recorded.   The  result  of  calcu- 
lating T]j  »  using  the  average  of  different  numbers  of  probes  is  shown. 
The  inclusion  of  only  the  four  probes  giving  the  highest  temperatures  was 
seen  to  give  reasonable  agreement  between  the  two  efficiencies  at  speeds 
up  to  lU,000  r.p.m.   It  was  again  observed  that  neither  the  peak 
efficiency  nor  the  shape  of  the  variation  of  the  efficiency  with  flow 
rate  depended  on  the  speed  for  speeds  up  to  1^,000  r.p.m.   Less 
consistency  is  shown  in  the  results  for  j\1   ,  than  for  T\       at  15,000  r.p. 
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In  the  above  tests,  the  stagnation  pressure  at  station  h   was  taken 
as  the  manifold  average  (equivalently  the  area  average)  of  four  fixed 
probes  as  had  been  done  previously.   For  the  reasons  given  in  the 
following  section,  the  above  measurements  of  efficiency  are  given  only  to 
illustrate  the  points  discussed,  and  no  results  are  shown  for  the 
pressure  ratio  itself. 

3.6.  Evaluation  of  the  Performance  with  Highly  Non-Uniform  Exit 
Conditions 

When  the  exit  impact  pressure  probe  readings  were  recorded 
individually,  the  presence  of  flow  separation  or  near  separation  at  the 
diffuser  exit  was  apparent.   Figure  18  shows  a  sketch  made  of  oil 
deposits  found  on  the  surfaces  of  the  diffuser  when  the  compressor  was 
dismantled.   Figure  19  shows  the  orientation  of  the  station  k   instru- 
mentation with  respect  to  the  diffuser  blades ,  and  clearly  the  two 
pressure  probes  near  a  second  row  pressure  surface  would  be  expected  to 
read  high  compared  to  the  two  probes  near  a  suction  surface,  as  was 
observed.   It  was  instructive  to  examine  different  methods  of  evaluating 
the  efficiency  from  the  available  measurements.   Figure  20  shows  the 
pressure  ratio  and  efficiency  T]   based  on  the  area  average  of  only  the 
two  highest  pressures,  on  the  area  average  of  all  four  probes,  and  on  the 
mass  average  of  all  four  probes.  Also  shown  is  the  efficiency  Re- 
evaluated using  an  assumed  peripheral  distribution  of  impact  pressure 
between  the  probes.   These  methods  are  described  in  Appendix  F.  While 
there  was  good  agreement  between  the  two  methods  in  which  the  continuity 
equation  was  satisfied  at  station  4,  there  was  an  unacceptable 
uncertainty  in  the  efficiency  evaluated  by  any  of  these  methods,  since 
the  flow  profiles  would  be  expected  to  change  as  changes  were  made  in 
diffuser  blade  settings.   The  overall  static-to-static  efficiency  of  the 
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compressor,  7]   (defined  by  Eq.  (19))  which  includes  the  losses  in  the 
honeycomb  section  and  in  the  exit,  is  seen  in  Figure  20  to  be  very  low. 
It  was  thought  (incorrectly)  that  TL  might  not  be  a  sensitive  indicator 
of  improvements  that  could  be  effected  in  the  diffuser  performance. 

In  order  to  obtain  an  accurate  measurement  of  the  total-to-total 
stage  efficiency,  twelve  holes  were  drilled  radially  between  the  stator 
exit  plane  and  the  flow  straightener  to  allow  Kiel  probe  surveys  to  be 
carried  out.  The  peripheral  displacements  of  the  holes  are  shown  in 
Figure  19.  The  holes  were  so  distributed  that  each  of  the  two  sets  of  six 
radial  surveys  gave  the  distribution  of  stagnation  pressure  at  the  stator 
exit  over  a  representative  average  blade  passage  as  shown  in  Figure  F  3  •  A 
machined  gauge  block  (seen  in  Figure  29)  allowed  the  probe  to  be  positioned 
accurately  in  0.050"  radial  increments.  A  comparison  of  results  from  the 
two  sets  of  six  surveys  indicated  the  success  of  this  process  of  averaging. 

Surveys  were  carried  out  for  two  diffuser  configurations .  The 
integrals  in  Eq.  F(32)  were  evaluated  using  subroutine  DATINT  (reference  k) 
and  the  efficiency  was  calculated  using  Eq.  F(35).   The  distributions  of 
impact  pressure  from  the  two  sets  of  six  surveys  were  observed  to  be 
qualitatively  similar,  and  the  two  values  calculated  for  the  efficiency 
were  different  by  less  than  1%.     However,  the  efficiencies  calculated  from 
the  surveys  for  both  configurations  were  lower  than  that  given  by  any  of  the 
three  methods  applied  to  the  four  fixed  probe  readings .  An  example  is  shown 
in  Figure  21.  The  improved  agreement  between  the  area  and  mass -averaged 
methods  in  Figure  21  compared  to  the  results  shown  in  Figure  20  was  thought 
to  be  due  to  the  change  in  profile  shape,  and  this  further  illustrated  the 
uncertainty  in  the  use  of  only  four  probes. 

The  shape  of  the  impact  pressure  profiles  measured  in  these  tests 
(Figure  21)  suggested  that  the  analytic  representation  of  the  profile  by  a 
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cubic  function  might  reasonably  be  applied  in  the  radial  direction.   There- 
fore twelve  Kiel  probes  were  inserted  in  the  survey  holes  (Figure  19)  to 
depths  required  to  provide  the  pattern  of  coverage  shown  in  Figure  F  3- 
A  test  was  conducted  to  determine  whether  the  machine  performance  was 
affected  by  the  presence  of  the  fixed  probes,  and  to  compare  the  performance 
evaluated  from  the  twelve  fixed  probes  with  that  from  the  four  fixed  probes. 
The  results  are  shown  in  Figure  22.   During  part  1  of  the  test,  the  twelve 
Kiel  probes  and  a  cobra  probe  at  the  rotor  exit  were  removed  and  the  perfor- 
mance was  measured  based  on  the  four  fixed  probes.   During  part  2,  the 
twelve  Kiel  probes  at  the  stator  exit  and  the  cobra  probe  at  the  rotor 
exit  were  installed  and  the  performance  was  measured  again.   Figure  22  shows 
that  the  presence  of  the  extra  probes  did  not  measurably  affect  the  perfor- 
mance. Also,  consistent  with  the  result  shown  in  Figure  21  for  the  survey 
method,  the  efficiency  evaluated  from  the  twelve  Kiel  probes  according  to 
the  method  given  in  Appendix  F  3.3  (the  "radial  profile  method")  is  2-3$>  lower 
than  that  evaluated  from  the  four  fixed  probes.   It  was  also  observed  that 
the  efficiency  based  on  the  mass  averaged  stagnation  pressure  from  the 
twelve  probes  agreed  closely  with  that  given  by  the  radial  profile  method 
applied  to  the  same  readings . 

It  was  concluded  that  the  effect  of  changes  in  the  diffuser  geometry 
could  best  be  evaluated  from  an  examination  of  the  pressure  ratio  and 
efficiency  calculated  from  the  twelve  Kiel  probe  readings,  using  the  radial 
profile  method  described  in  Appendix  F  3-3' 
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k.      Performance  at  lU,000  R.P.M.  with  Different  Diffuser  Geometries 
k.l     Description  of  the  Tests 

The  performance  of  the  Hybrid  compressor  was  measured  at  lU,000  r.p.m. 
for  chosen  changes  in  the  geometry  of  the  stator  blading.   The  high  and  low 
solidity  bladings,  designated  here  and  in  reference  1  as  diffuser  A  and 
diffuser  B  respectively,  were  installed  and  tested  in  turn.  Data  were 
obtained  for  three  combinations  of  blade  angles  of  row  1  and  row  2  for 
diffuser  A  and  two  combinations  for  diffuser  B.   The  axial  separation  between 
rotor  and  stator  was  1.0  -  .05"  and  between  row  1  and  row  2  was  0.2  -  .05". 
The  leading  edges  of  the  second  blade  row  were  displaced  peripherally  with 
respect  to  the  trailing  edges  of  the  first  row  as  shown  in  Figure  23.   There 
was  consequently  a  greater  axial  separation  between  the  trailing  edge  plane 
of  diffuser  B  and  the  probes  at  station  k,   than  for  diffuser  A. 

The  instrumentation  was  held  fixed  and  the  same  procedure  was  followed 
in  each  test.   The  instrumentation  in  the  exit  flow  was  that  described  in 
section  3*6,  with  Kiel  probes  as  shown  in  Figure  19.  A  cobra  probe  measured 
the  flow  angle  in  the  peripheral  direction  and  the  impact  pressure  on  the 
centerline  of  the  annulus  at  the  rotor  exit  for  each  performance  data  point. 
The  static  pressure  at  station  2  was  measured  at  the  outer  wall  0.5"  downstream 
of  the  rotor. 

The  test  procedure  was  to  obtain  data  scans  at  lU,000  r.p.m.  from  open 
throttle  to  surge  conditions,  adjusting  the  cobra  probe  before  each  point 
was  recorded.  Radial  surveys  were  then  taken  using  the  cobra  probe  at  the 
rotor  exit  for  two  or  three  throttle  positions.  Finally,  with  the  compressor 
operating  close  to  the  highest  efficiency  measured,  a  small  amount  of  a 
lampblack  and  oil  mixture  was  bled  into  the  inlet  at  station  0.  After 
several  minutes  to  allow  drying,  the  compressor  was  shut  down  and  opened  for 
inspection. 
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k.2     Performance  Results 

The  results  from  six  tests  of  five  different  configurations,  from  which 
complete  data  were  obtained,  are  tabulated  in  Appendix  I.   In  Table  I-I 
performance  parameters  of  the  compressor  and  of  the  diffuser  separately 
are  listed.  Data  from  this  table  only  will  be  discussed.   Table  I-II  lists 
other  performance  parameters  of  the  compressor  and  compares  the  efficiencies 
evaluated  by  different  methods .   Table  I-III  is  a  listing  of  input  data  for 
these  tests.   It  should  be  noted  that  the  pressure  ratio  listed  in  Table  I-I 
(and  discussed  below)  was  obtained  using  Eq.  F(38),  whereas  the  pressure 
ratio  in  Table  I-II  is  the  mass  averaged  value  defined  by  Eq.  F(3)«   The 
difference  is  seen  to  be  very  small. 

The  performance  of  the  compressor  with  the  different  diffuser  geometries 
is  illustrated  in  Figure  2U  and  Figure  25,  which  show  the  total-to-total 
pressure  ratio  and  efficiency  respectively  as  a  function  of  referred  flow 
rate.   It  can  be  seen  that  the  pressure  increases  obtained  with  the  higher 
solidity  diffuser  A  were  larger,  by  as  much  as  10$,  than  with  diffuser  B  at 
correspondingly  higher  efficiencies.  With  the  knowledge  that  very  large 
variations  from  hub  to  tip  existed  in  the  angle  of  the  flow  leaving  the 
rotor,  this  result  was  not  unexpected,  and  accounts  for  the  brevity  of  the 
tests  conducted  with  diffuser  B. 

In  reference  1  the  blade  angles  selected  for  diffuser  A  were  56.17  and 
31.83  as  defined  in  Figure  23.  When  it  was  found  that  a  significant  flow 
separation  occured  in  the  second  row  at  these  settings,  the  loading  of  the 
first  row  was  increased  (by  decreasing  the  blade  angle  to  52.17  )  in  order 
to  reduce  the  loading  on  the  second  row.  Larger  flow  rates  were  obtained 
at  similar  efficiencies  over  a  wider  range  with  this  configuration.  A 
further  improvement  was  made  by  further  reducing  the  loading  on  the  second 
row  by  increasing  the  blade  angle  to  3^.5  •   The  latter  configuration,  for 
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which  the  data  from  RUN  1+9  are  shown,  give  a  maximum  total-total  stage  effi- 
ciency of  79$.   It  also  gave  the  highest  values  of  static -static  overall 
efficiency  for  the  compressor. 

Figure  26  shows  the  performance  of  the  diffuser  calculated  from  the 
available  measurements.   The  lower  figure  shows  the  pressure  recovery  factor 
defined  in  Eq.  (21 )  as  a  function  of  the  referred  flow  rate.   The  upper 
figure  shows  the  diffuser  efficiency  defined  by  Eq.  (23)  as  a  function  of 
the  recovery  factor.   It  was  observed  that  if  the  diffuser  efficiency  was 
plotted  as  a  function  of  the  flow  rate  it  deci eased  with  decreasing  flow 
rate  from  a  flat  maximum  close  to  open  throttle.   The  upper  curves  are 
typical  of  a  plot  of  (l  -  Q) ,  where  £  is  the  loss  coefficient,  against 
blade  loading  for  an  axial  cascade.   The  curves  show  that  the  losses 
in  all  diffusers  increase  rapidly  at  high  incidence  angles  from  the  rotor, 
corresponding  to  closing  throttle.   The  best  compressor  efficiencies  were 
obtained  with  diffusers  for  which  the  rapid  increase  in  losses  occured  at 
the  greatest  values  of  the  recovery  factor.  However,  the  levels  of  the 
curves  in  Figure  26  can  not  be  directly  interpreted  as  showing  the  relative 
merits  of  the  diffusers.   The  flow  from  the  rotor  is  a  complicated  one  and 
the  single  impact  pressure  on  the  centerline  at  Station  2  is  not  a  sufficient 
definition  of  the  flow  entering  the  diffuser.   There  was  a  measurable  effect 
of  the  diffuser  geometry  on  the  flow  from  the  rotor,  as  will  be  shown  in 
the  following  section. 

The  blockage  factors  listed  in  Table  I-I  are  the  result  of  applying  the 
analytical  methods  of  Appendix  F  to  the  probe  measurements .  C/-  and  C>  contain 
useful  physical  information.  C/-,  given  in  Eq.  F(37)?  is  the  product  of  the 
cosine  of  an  average  flow  angle  at  station  h,    &*,   and  an  area  blockage  factor 
C/r',  shown  in  Figure  F  3-   If  C/-'  is  assumed  not  to  change  with  throttling, 
the  variation  of  C^  corresponds  to  the  variation  of  the  average  flow  angle. 
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Indeed,  C/-  is  observed  to  decrease  with  decreasing  flow  rate  (and  increasing 
diffuser  pressure  recovery  factor)  and  to  have  values  corresponding  to 
expected  values  of  cos  a  .   Thus  C/-'  is  close  to  unity,  and  the  variation 
in  C/-  can  be  examined  to  find  the  variation  in  the  cosine  of  the  flow  angle 
at  station  h.      The  accuracy  of  the  procedure  would  increase  if  the  number 
of  probes  defining  the  exit  profile  were  increased.  C.  is  a  blockage  factor 
that  relates  the  average  impact  pressure  to  the  static  pressure  through  the 
continuity  equation  at  station  k   (Eq.  F(6)).   C.  depends  on  the  shape  of 
the  flow  profile  and  hence  can  be  examined  as  an  indicator  of  separation. 
k.3     Measurements  of  the  Flow  at  the  Exit  of  the  Stator 

Peripheral  distributions  of  impact  pressure  measured  by  the  twelve  Kiel 
probes  at  station  k   are  shown  in  Figure  27  for  what  was  indicated  to  be  the 
best  diffuser  configuration  tested  (Diffuser  A).   The  curves  are  plotted 
as  if  the  probes  were  located  in  the  same  blade  passage  (see  Figure  F3). 
Distributions  at  h  =  0.25  and  h  =  0.75  are  shown  for  five  flow  rates  from 
open  throttle  to  near  surge.   There  is  little  change  in  the  relative  shapes 
of  the  curves  over  the  mass  flow  range.   The  distribution  at  the  larger 
radius  is  reasonably  flat,  whereas  on  the  hub  side  of  the  centerline  a 
separating  profile  is  evident.  Similarly  shaped  profiles  were  obtained  with 
the  two  other  configurations  of  this  diffuser. 

Results  from  Diffuser  B  given  in  Figure  28  show  a  separated  profile 
near  the  hub  on  the  suction  side  of  the  blade,  and  the  profile  at  the  outer 
radius  is  no  longer  flat.   The  poorer  performance  obtained  with  this  diffuser 
is  evident  in  those  profiles. 
k .k     Surface  Flow  Visualization  Results 

Following  the  incidental  oil  flow  patterns  obtained  at  unknown  conditions 
during  the  preliminary  tests,  the  hub  surface  of  the  compressor  was  painted 
white  and  surface  patterns  were  obtained  at  chosen  throttle  settings  as 
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described  in  section  k.l.     Some  examples  of  the  results  are  shown  in  Figure 
29  and  Figure  30  which  are  photographs  taken  after  a  test  with  Diffuser  A 
with  blade  angles  of  52.17  and  31.83  .   In  the  upper  photograph  of  Figure 
29  the  oil  path  describes  clean  lines  over  the  hub  surface  that  show  no 
signs  of  flow  separation.   In  the  lower  figure,  it  can  be  seen  that  the 
flow  angle  at  the  surface  can  be  measured.   This  angle  was  found  to  agree 
well  with  the  flow  angle  indicated  by  cobra  probe  surveys.  A  separation 
at  the  hub  in  the  second  row  of  blades  is  evident,  and  the  patterns  on  the 
blades  themselves,  shown  in  Figure  30(a),  confirm  this.   The  clean  surface 
around  the  blade-hub  junction  is  presumably  the  result  of  inertial  separation 
of  particulates.  A  sketch  of  these  patterns  is  given  in  Figure  31.   The 
elaborate  pattern  in  the  region  between  the  blade  rows  suggests  that  the 
row  orientation  (here  set  at  5  x  blade  spacing)  might  have  a  strong  effect 
on  the  separation  pattern.  Controlled  experiments  with  this  parameter  as 
a  variable  were  not  carried  out,  but  for  all  configurations  of  all  diffusers 
tested  the  presence  of  separation  was  evidenced  in  both  oil  patterns  and 
in  exit  Kiel  probe  measurements . 

The  surface  patterns  obtained  in  the  rotor,  Figure  30(b),  can  not  easily 
be  interpreted,  however  one  feature  should  be  noted.   Particulate  matter 
collected  in  the  pressure-side  corner  of  the  blades  near  the  exit  plane. 
This  was  true  for  all  diffuser  configurations  tested  throughout  the  investi- 
gation of  the  Hybrid  compressor.   The  particulates  could  easily  be  brushed 
away,  and  so  the  presence  of  at  least  a  small  region  of  very  low  energy  flow 
in  the  corner  of  the  rotor  is  suspected. 
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5.   The  Flow  from  the  Rotor 

5.1  Introduction 

The  rotor  of  the  Hybrid  compressor  was  designed  using  available  methods 
for  radial  compressors.   The  diffuser  was  designed  using  available  methods 
for  axial  cascades,  taking  the  conditions  calculated  on  the  mean  radius  at 
the  exit  of  the  rotor  as  being  representative  of  the  conditions  along  the 
blade  height  (reference  l).  Early  surveys  at  the  rotor  exit  (Figure  13) 
had  indicated  that  large  variations  in  the  time -averaged  flow  properties 
existed.   Therefore,  during  the  performance  measurements  reported  in  Section 
h,   cobra  probe  surveys  were  carried  out  in  order  to  investigate  the  effect 
of  changes  in  the  diffuser  geometry  on  the  flow  leaving  the  rotor.   This 
information  would  be  required  as  an  input  to  the  design  of  an  improved 
diffuser. 

5.2  Probe  Measurements  and  Data  Reduction 

For  each  of  the  diffuser  configurations  tested  and  reported  in  Section 
U,  a  radial  cobra  probe  survey  was  carried  out  at  the  rotor  exit  for  a 
throttle  position  close  to  optimum  efficiency.   One  or  two  further  surveys 
were  made  at  other  throttle  settings  selected  arbitrarily.   The  cobra 
probe  gave  the  peripheral  flow  angle  (ou)   at  each  radial  position  and  the 
local  impact  pressure.   This  information  was  reduced  to  velocity  using  the 
assumptions  that  the  flow  was  in  cylindrical  planes  and  that  the  static 
pressure  at  the  outer  wall  was  constant  across  the  annulus .   The  data  reduc- 
tion is  described  and  the  Monroe  Calculator  Card  Program  No.  210  is  given 
in  Appendix  G.  Results  from  the  performance  tests  reported  in  Section  h 
are  listed  in  Table  I-IV.  Reference  can  be  made  by  RUN  and  POINT  numbers 
to  Tables  I-I,  I-II  and  I-III  for  the  test  conditions  represented  by  this 
data. 
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5.3  Results  for  the  Best  Diffuser  Configuration 

As  an  illustration  of  the  data  obtained  for  the  rotor,  plots  of  the 
results  from  the  tests  which  gave  the  best  performance  are  shown  in  Figure 
32  and  Figure  33- 

Figure  32  shows  the  variation  of  the  flow  angle  on  the  centerline,  6u    , 
as  a  function  of  the  referred  flow  rate.   Imposed  on  the  figure  are  the 
distributions  of  the  flow  angle  from  hub  to  tip  at  those  flow  rates.  What 
is  immediately  evident  in  this  figure  is  that  the  variation  of  the  flow 
angle  from  hub  to  tip  at  the  diffuser  entrance,  at  any  throttle  position, 
is  greater  than  the  variation  in  the  centerline  flow  angle  over  the  compressor 
operating  range.  As  can  be  seen  here  in  the  flow  angle  measurements,  the 
distribution  of  properties  generally  were  not  qualitatively  dependent  on 
the  throttle  position. 

The  distributions  of  properties  for  a  flow  rate  close  to  that  giving 
the  highest  efficiency  (Point  23  of  Run  U9)  are  shown  in  Figure  33-   (in 
this  and  subsequent  figures,  points  very  close  to  the  walls  have  been  omitted 
because  the  influence  of  the  wall  on  the  measurements  was  unknown.)  It 
can  be  seen  in  Figure  33(a),  which  shows  the  absolute  and  relative  velocity 
distributions,  that  relative  to  the  rotor,  the  velocity  is  much  larger  over 
the  outer  than  over  the  inner  half  of  the  rotor  exit.  Near  the  inner  surface 
the  velocity  is  reduced  but  the  profile  is  not  approaching  separation, 
confirming  the  evidence  of  the  oil  flow  observations.   The  outer  half  of 
the  absolute  velocity  profile  is  quite  flat,  as  is  the  outer  third  of  the 
flow  angle  profile  (c^)  shown  in  Figure  33(b).  From  the  slip  factor  varia- 
tion shown  in  Figure  33(c),  (or  equivalently  the  angle  [3p  variation  shown 
in  Figure  33(b))  the  flow  relative  to  the  rotor  is  found  to  be  within  10 
of  being  radial  except  over  the  inner  25%   of  the  blade  height.   From  these 
observations,  clearly  the  separation  observed  in  the  second  row  of  diffuser 
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"blades  is  the  result  of  the  depleted  velocity  profile  and  increased  flow 
angles  in  the  inner  half  of  the  annular  passage. 

Figure  33(d)  shows  the  distribution  of  a  quantity  (I  ,  defined  by 
Eq.  G(6)),  that  is  proportional  to  the  mass  flux  per  unit  passage  depth. 
From  the  integral  of  the  distribution  and  the  value  on  the  annulus  center- 
line  (using  Eq.  G(7)),  the  blockage  factor  k  ,  was  evaluated  to  have  the 
value  O.83.  A  similar  procedure  applied  to  data  from  two  other  configura- 
tions at  different  throttle  settings  gave  values  for  k  of  0.822  and  0.8^2. 

m 

Thus  k  =0.83  should  be  used  as  the  blockage  factor  if  centerline  values 
m 

are  taken  as  representative  over  the  passage  height.   It  is  noted  that  when 
the  mass  flux  distributions  were  integrated  using  Eq.  G(3),  the  integral 
was  found  to  agree  with  the  measured  mass  flow  rate  to  within  3^fo   and  hiff0 
in  two  different  cases.   This  reasonable  agreement  is  in  support  of  the 
assumptions  made  in  the  data  reduction. 
5.U  Design  Relations  From  the  Results  of  All  Diffuser  Configurations 

Comparison  of  the  rotor  exit  profiles  in  Figure  33  with  those  in  Figure 
13  suggested  that  there  was  no  qualitative  change  in  going  from  10,000  to 
lU,000  r.p.m.  with  a  movement  of  the  diffuser  blades  away  from  the  rotor 
and  with  a  change  in  the  first  row  blade  angles  of  h    .     Equally,  correspond- 
ing profiles  for  the  surveys  during  the  tests  of  Section  h   were  also  quali- 
tatively similar.   Therefore,  representations  of  the  profiles  at  the  rotor 
exit  were  sought  which  could  be  used  as  the  basis  for  the  design  of  a  new 
diffuser.   The  useful  result  of  this  effort  appears  in  the  ten  parts  of 
Figure  3^» 

The  design  of  an  axial  cascade  with  spanwise  gradients  requires  the 
specification  of  the  initial  velocity  in  magnitude  and  direction.   The  extent 
to  which  the  measured  radial  distributions  of  the  absolute  velocity  (Vp) 
and  peripheral  flow  angle  (a-,)  are  unaffected  by  changes  in  the  diffuser 
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geometry  can  be  seen  in  Figure  3^.  Figure  3^  shows  these  distributions 
normalized  with  repsect  to  the  centerline  values  (V   and  ou    )  which  were 
obtained  from  plots  similar  to  those  in  Figure  33.   It  should  be  noted 
that  the  normalized  profiles  of  velocity  and  flow  angle,  which  are  required 
for  the  diffuser  design,  were  closer  to  a  single  analytical  representation 
than,  for  example,  profiles  of  the  axial  velocity  component.  Shown  in 
Figure  3^?  in  broken  lines,  are  suggested  approximations  to  the  profiles 
which  can  be  used  as  input  to  the  design  a  diffuser  for  this  unique  rotor. 
The  approximations  shown  have  been  chosen  to  be  analytically  convenient, 
but  also  to  best  fit  the  data  obtained  with  the  different  configurations 
operating  near  optimum  efficiency.   (The  points  nearest  optimum  efficiency 
are  indicated  by  an  asterisk  in  each  legend.) 

The  suggested  approximations  for  the  distributions  of  velocity  (V9) 
and  flow  angle  (ou )  normalized  with  respect  to  their  centerline  values, 
the  normalized  quantities  being  denoted  by  V     and  ou   respectively,  are  as 


V 

V2   =  1   -   2(h   -  if 

0  <.  h   £  0.5 

v2  =  1.0 

0.5   £  h   £  1.0 

V 

S^  -  1.2  -  O.Uh 

0  £  h  <.  0.7 

Qb  =  °-92 

0.7  ^  h   <.  1.0 

(2k) 


(25) 


As  shown  in  Section  5«3>  the  blockage  factor  k  ,  which  appears  in  the 
continuity  equation  expressed  in  terms  of  centerline  properties,  should  be 
given  the  value  O.83.   The  value  which  would  be  obtained  using  the  above 
approximations  in  Eq.  G(7),  using  Eq.  G(6)  and  the  relations  given  in 
Table  Al,  would  not  include  the  effect  of  the  boundary  layers  on  the  inner 
and  outer  walls . 
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The  angular  deflection  through  the  diffuser  blading  A  a,  which  is 
required  to  produce  a  uniform  flow  angle  at  station  h,    a  ,  can  be  expressed 
in  terms  of  the  above  Eq.(25)  as  follows: 

Aa(h)  =  o^  -  o^.(h).  o^ 


Mr 


but  since 


then 


U*)m  =  ^  -  % 


Aa(h)  =  (Aff)m  +  o^m(l  -  a2)  (26) 

Using  Eq.(25),  the  required  radial  variation  in  the  total  deflection  is 
given  by. 


A  a   =   (A  a)      +  0.2  (2h   -  l) .    a^  for  0   ^  h   ^  0.7 

A  a   =   (  A  a)      +  0.08   cu     for  0.7   £  h   £  1.0 
m  <dra 


(27) 
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6.   Discussion  and  Recommendations 

The  prototype  Hybrid  compressor  was  designed  using  available  techniques 
for  both  radial  and  axial  compressors.   It  was  anticipated  however,  (refer- 
ence l),  that  the  diffuser  predicted  to  give  the  best  performance  in  a  uniform 
flow  (Diffuser  B)  might  not  be  the  best  design  to  follow  a  rotor  wherein 
secondary  flow  effects  were  expected  to  be  large,  but  could  not  be  predicted. 
Higher  solidity  tandem  cascades  were  therefore  provided  (Diffuser  A),  and 
the  program  of  tests  reported  here  has  confirmed  that  the  latter  design 
gave  the  better  performance.  A  maximum  efficiency  of  79$  was  obtained,  and 
it  is  possible  that  this  could  be  improved  slightly  by  a  change  in  blade 
angle  or  blade  row  orientation. 

However,  measurements  of  the  time  averaged  flow  at  the  exit  of  the 
rotor  have  shown  that  there  is  a  very  large  (15  -  20  )  variation  from  the 
hub  to  the  tip  in  the  angle  of  the  flow  approaching  the  diffuser.  Exit 
impact  pressure  profiles  and  oil  flow  patterns  on  the  blade  surfaces  have 
shown  that  for  all  configurations  examined,  the  flow  was  separated  in  the 
second  blade  row.   The  separation  would  be  expected  as  a  consequence  of 
too  high  an  incidence  near  the  hub  when  the  tip  was  operating  within  design 
limits.  Because  of  the  magnitude  of  the  variation  in  the  flow  angle,  which 
is  a  property  of  this  rotor,  there  is  no  angle  at  which  straight  blades 
can  be  set  in  order  to  eliminate  separation  altogether. 

The  flow  from  the  rotor,  while  containing  some  irregularities  evidenced 
by  dirt  deposits  in  the  pressure-side  corners  near  the  outlet,  had  well 
defined  time -averaged  profile  shapes  that  did  not  depend  significantly  on 
the  diffuser  geometry.  While  retarded  on  the  inner  wall  or  hub,  the  flow 
appeared  not  to  be  separated  on  the  hub  within  the  rotor  nor  downstream 
of  its  exit.   It  may  be  possible  therefore  to  design  a  system  of  blading 
which  would  turn  the  flow  axially  without  significant  separation.   The 
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first  row  of  blading  would  have  either  twisted  or  tapered  blades,  and  near 
the  hub  more  than  two  rows  would  probably  be  needed.   The  analytical 
representations  of  the  velocity  and  flow  angle  given  in  Section  5.k   are 
adequate  for  design  purposes.  Since  there  is  at  least  a  measurable  effect 
of  the  diffuser  geometry  on  the  rotor  exit  profiles,  further  refinement 
of  these  expressions  is  not  justified. 

It  is  the  author's  opinion  that  while  an  improved  diffuser  can  be 
designed  using  the  information  given  here,  a  meaningful  analytical  estimate 
of  its  efficiency  can  not  be  made.   The  performance  calculated  in  reference 
1  for  straight  diffuser  blades  in  a  flow  that  is  assumed  to  be  uniform  and 
free  of  separation  must  represent  an  upper  limit  to  the  performance  which 
would  be  calculated  for  a  flow  with  initial  gradients,  using  available  data. 
82. h°J0  was  calculated  in  reference  1  for  the  compressor  efficiency  operating 
on  design  with  the  blading  of  Diffuser  B.   79$  "was  measured  for  the  best 
configuration  of  Diffuser  A  that  was  tested.   The  latter  efficiency  was 
measured  with  exit  profiles  that  are  shown  in  Figure  27.   The  worst  diffuser 
configuration  gave  a  maximum  efficiency  of  only  75$  with  the  exit  profiles 
shown  in  Figure  28.   There  was  a  difference  of  10$  in  the  stagnation  pressure 
rise  between  the  two  extreme  cases  as  shown  in  Figure  2j+.   From  an  examina- 
tion of  these  differences  it  is  speculatively  concluded  that  an  improvement 
of  2$  in  efficiency  could  be  achieved  if  the  diffuser  separation  were 
eliminated  and  the  exit  pressure  profiles  correspondingly  flattened.  However, 
it  should  also  be  noted  that  the  efficiencies  measured  in  the  tests,  while 
evaluated  in  a  consistent  manner  for  the  various  diffusers  tested,  probably 
depart  by  some  small  unknown  amount  from  the  true  compressor  efficiency. 
If  1$  is  accepted  as  this  uncertainty,  based  on  all  available  information, 
an  efficiency  of  80-83$  would  be  expected  for  the  Hybrid  compressor  using 
the  present  rotor  and  a  new  diffuser  wherein  the  separation  was  eliminated. 
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A  complete  analysis  of  the  information  obtained  in  the  test  program 
has  not  been  possible  and  consequently  the  data  have  been  included  here  in 
the  tables  of  Appendix  I. 

A  continuation  of  the  Hybrid  compressor  development  should  include  the 
following: 

a)  Design  and  test  one  or  possibly  two  diffusers  using  the  design 
expressions  given  in  Section  ^.h.     Assume  cylindrical  stream  surfaces 

b)  Extend  the  region  between  the  stator  exit  and  the  entrance  to  the 
exit  honeycomb  to  reduce  non-uniformities  in  the  flow  where  the 
efficiency  is  to  be  established.   (This  might  be  necessary  in  order 
to  accommodate  a  third  blade  row). 

c)  Analyze  the  rotor  performance  using  the  methods  given  in  reference 
1.   This  will  establish  characteristics  of  this  rotor  to  compare 
with  more  conventional  radial  wheels.   The  required  data  are  given 
in  the  tables  of  Appendix  I. 
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and  Flow  Angle  at  the  Rotor  Exit 


76 


a 

B 

o 


I 


co 

W 
to 

I 


o 


1.0 

— _ , 

DATA  FROM  RUN  52  -  DIFFUSER  B 

£> 

(SEE  FIGURE  2k  AND  2$) 

PT.NO.       W/B  (LBS/SEC) 

O          12               U.662 

*      £          13              3.716 

3D 

0.5 

3 

qrf 

9a' 

• 

'& 

1.0 


0.5 


o.5 


1.0 


V  V2m 


& 
& 


O 
O 


&  o 


&\0 


aa 


sr 


«  -  DATA  POINT  CLOSEST  TO 
MAXIMUM  EFFICIENCY 


SUGGESTED 
APPROXIMATIONS 
(EQ.2U  AND  EQ.25) 


\   o 

\ 

\ 

to    o 
\  o 

V     o 

i jj 


0o8 


1.0 


1.2 


°V«**i 


1.5 


l.U 


FIGURE  3k.     Normalized  Distributions  of  Velocity 
(continued) 
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APPENDIX  A.  SOME  USEFUL  COMPRESSIBLE  FLOW  EXPRESSIONS 

In  the  design  and  analysis  of  turbomachines ,  and  in  many  propulsion 
problems,  it  is  desirable  to  deal  only  with  those  properties  of  the  flow 
which  of  necessity  enter  into  the  design,  and  those  which  can  be  measured. 
Therefore  we  would  like  to  have  equations  that  relate  the  velocity  to 
stagnation  pressure,  static  pressure  and  stagnation  temperature,  but  would 
choose  to  avoid  the  use  of  density  and  static  temperature. 

There  are  two  ways  of  non-dimensionalizing  velocity  that  avoid  the  use 
of  the  speed  of  sound,  which  depends  on  the  static  temperature  -  either  the 
speed  of  sound  at  the  stagnation  temperature  (a, ),  or  the  "stagnation"  or 
"limiting"  velocity  (V,  )  can  be  introduced.  While  the  choice  is  almost 
arbitrary,  the  limiting  velocity  has  a  clearly  defined  significance  even 
for  real  (non-perfect)  gases,  and  its  use  results  in  the  tidiest  form  of 
the  equations.  The  smaller  dimensionless  quantities  which  result,  present 
no  problem  when  a  digital  calculator  or  computer  is  used  to  perform  computations 

The  "limiting"  velocity  is  the  maximum  velocity  attainable  in  an  adiabatic 

expansion  from  the  prevailing  state  to  zero  temperature  and  pressure.  The 

steady-flow  energy  equation  gives,  for  any  gas  where  gravitational  effects 

are  neglected, 

2. 
\   =  h  +  £■  A(l) 

where  h  is  enthalpy  per  unit  mass,  V  is  velocity  and  subscript  t  , 
here  and  hereafter,  denotes  the  stagnation  value. 

By  the  above  definition,  for  any  gas,  V,  =  *j2h,   '  ,  but  for  a  gas 
with  constant  specific  heats,  if  T  is  the  temperature, 


Vt  =  Src-Y?  A(2) 


p 
78 


For  a  perfect  gas,  the  stagnation  speed  of  sound,  a,  ,  is  given  by 

at  «  ^VR*J  =yy  ;yt  A(3) 

where  R  is  the  gas  constant.   [Note  that  for  air,  with  y  =  1.^02  , 
V.  =  109.62  */  T  '  (ft/sec),  where  T,   is  in  °Rankine.] 

If  the  velocity  referred  to  the  limiting  velocity  is  defined  as 

X  =  i  A(V) 

t 

reference  h9   Page  6,  gives  the  following  relationships  (which  are  easily 
derived)  between  the  local  and  isentropic -stagnation  properties  of  a  perfect 
gas: 

f=l-X2  A(5) 

^-=(i_x2)Y-1  A  (6) 

1 

-L=  (l-x2^-1  A(7) 


where  p  denotes  pressure  and  p  denotes  density. 

By  using  these  expressions,  all  compressible  flow  equations  can  be  written 
in  terms  of  pressures,  velocities  and  stagnation  temperature.   (Note  that  in 
the  following  development  the  units  are  not  included.  The  final  expressions 
must  be  expressed  in  consistent  units.) 

The  equation  of  continuity,  giving  the  mass  flow  rate  (m)  per  unit  area 
(A),  may  be  written  as 


£  =  (_£_)  ( JL)  D  V 

Pt 


A   ^'  *▼'  Vt 


=  X  (l-X2)^37  5^/25^7  A(8) 
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RT.  v   pt 


where  we  have  used  the  perfect  gas  equation  of  state,  Eq.  A(k)   and  Eq.  A(7). 


1 


Defining  F(X)  =  X  (l-X2)Y~  A (9) 


and  using  the  relationship  between  the  gas  constant  and  ratio  of  specific 
heats , 


I  -  W  •'<&>•£•     .  Ado) 


which  gives  the  mass  flux  in  terms  of  the  velocity,  the  stagnation  pressure 
and  stagnation  temperature . 

The  equation  of  continuity  may  also  he  written  as 

E  =  JL  (JL).V  =  J2-  (%    (JL)  v 

A   RT  KVtJ      t   RT.  VT;  VV  '      t 

or 

A   ^_x2;    Vlj  Vt  {L1) 

which  gives  the  mass  flux  in  terms  of  the  velocity,  the  static  pressure  and 
stagnation  temperature. 

Equation  A(lO)  and  Equation  A(ll)  can  he  used  depending  on  which 
quantities  are  known  in  any  particular  case.  Equation  A(ll)  can  be  viewed 
as  a  quadratic  for  the  velocity  in  terms  of  the  static  pressure,  stagnation 
temperature  and  the  mass  flux.  Thus  we  write  A(ll)  in  the  form 

X2  +  2BX  -  1  =  0  A(12) 

where  the  dimensionless  quantity,  B  ,  is  given  by 


B  = 


L2L  — £ — 
lv-li  V.  (A/A) 


A(13) 


If  the  mass  flux,  pressure  and  stagnation  temperature  are  known  on  a 
streamline,  then  the  velocity  is  given  by 
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X 


J 


B  +1  -  B 


A(lU) 


(since  the  second  root  is  negative). 

The  measurement  of  mass  flow  rate  by  means  of  an  orifice  plate  or  flow 
nozzle,  depends  upon  measurements  of  stagnation  pressure  and  temperature 
and  static  pressure.  From  Eq.  A(10),  using  Eq.  A(9)  and  Eq.  A(6),  the 
required  form  of  the  continuity  equation  is  then 


hW-M 


vt  w 


.  x 


Y 


A(15) 


Using  the  perfect  gas  equation  of  state  and  the  relationship  of  the  specific 
heat  to  the  gas  constant,  namely  (y-1)  C  =  YR  }  Eq.  A(15)  can  he  written  as 


«-^R 


.   X 


v 


A  (16) 


In  summary,  if  the  local  mass  flux  is  referred  to  the  "stagnation"  value, 
a  "total  flow  function"  is  obtained  which  relates  the  local  mass  flux  to 
the  maximum  value  of  the  mass  flux  attainable  from  the  given  state.  The 
"total"  flow  function,   #.  ,  from  Eq.  A(l6)  is 


m 
$  =  JL  =  -fil.  =  IjL 

*    ptvt    Ptvt    Vp 


A(17) 


Also,  using  Eq.  A(6) 


1 

Y-1 
it  =  X(l  -  X2)    =  F(X) 


A(l8) 
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The  total  flow  function  given  "by  Eq.  A  (17)  can  be  approximated  for  flows 
at  lower  speeds  by  expanding  the  right  hand  side  of  the  equation  in  terms  of 
the  difference  between  the  stagnation  and  static  pressures.  If  A  =  p.  -  p  , 
and  e  =  —  and  is  small,  Eq.  A (17)  may  be  written  as 


-  i-A 


i-i-« 


A(19) 


or 


I  Y  j 


3 


5Y+29 
k 96Y2  l 


2 
e  +  0 


P) 


A(20) 


Using  Eq.  A(20)  in  Eq.  A(l'5)  and  returning  to  the  original  notation, 


m 
A 


'2pt(PfP 
RT^ 


1-^\~h 


96y2; 


'V\     +0()3 


A(21) 


incompressible  r_        ..   ,   '  n 
=      il         LI  -  correction  terms  J 
mass  flux 


where  again,  the  units  of  measurement  must  be  consistent  in  the  use  of  this 
expression. 

These  results  are  summarized  in  Table  Al. 
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TABLE  Al.      COMPRESSIBLE   FLOW  OF  A  PERFECT  GAS   IN  TERMS 
OF  THE    "LIMITING"    (STAGNATION)  VELOCITY 


IEFINITION 


ISENTROPIC  RELATIONS 


(COMPARISON) 


■/ 


2C  T. 
P  t 


(   =  J  2h  '  for  real  gas ) 


X  =  v/v. 


T/Tt  =  1-X2       y 

P/Pt  «  (i-x  )  i 


M  = 


N  = 


V 


v 


/w? 


^Rf^      V  V-l 


FORMS   OF  THE  CONTINUITY  EQUATION ,  KNOWING : 


V,Pt,Tt    [Eq.   A (10)] 


V,p,T.     [Eq.  A(ll)] 


p,Pt,Tt    [Eq.   A(15)] 


w-i-H 


-i 


m 


vt 


x 
i-x2 


-2i.).JL 


v-ii    V. 


($M 


Jfci 

V 


FORMS  OF  THE  "TOTAL"  FLOW  FUNCTION 


1_ 
V-l 


PtVt 


,  ^,F(X,=X(1-Xf  ,   fc) 


SOLUTION  OF  THE  CONTINUITY  EQUATION  FOR  X: 


X  = 


7777. 


B    ,   from     X2   +  2BX  -  1  =  0 
where     B  = 


Y-l)    Vt(m/A) 


APPROXIMATION  OF  THE  MASS   FLUX  FOR  LOW   VELOCITIES: 


m 
A 


'2pt(yp> 

RTl 


v*f«- (3M¥)  ♦*  * 
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APPENDIX  B.   CALIBRATION  OF  THE  FLOW  NOZZLE 

The  flow  measuring  nozzle  in  the  intake  pipe  of  the  compressor  was 
calibrated  in  the  original  installation  shown  in  Figure  5(a).  Because  the 
installation  was  not  a  standard  one,  and  because  changes  in  the  downstream 
geometry  could  be  expected  to  influence  the  flange  tap  pressure  measure- 
ments, a  non-standard  method  of  calibrating  and  using  the  nozzle  was 
adopted. 

The  nozzle  geometry  and  the  installation  of  the  calibration  probe  are 
shown  in  Figure  B  1.   The  nozzle  was  designed  according  to  standards  for 
a  "normduse"  (reference  5).   The  nozzle  throat  to  pipe  area  ratio  is 
small  (.08917),  and  there  is  a  constant  diameter  section  which  is  equal  in 
length  to  the  radius  of  curvature  of  the  entrance.   For  this  geometry 
nearly  axial  streamlines  would  be  expected  in  the  throat,  and  the  static 
pressure  across  the  nozzle  near  the  exit  would  be  expected  to  be  nearly 
uniform.   In  view  of  the  difficulty  involved  in  measuring  the  static 
pressure  reliably  in  such  an  environment,  it  was  felt  that  the  assumption 
that  the  pressure  at  the  throat  taps  was  constant  across  the  flow,  was  to 
be  preferred,  and  this  was  done.   The  probe  shown  in  Figure  B  1,  consisted 
therefore  of  two  flattened  impact  pressure  tips  held  about  one  nozzle 
radius  apart  on  an  aerofoil  strut.   This  design  caused  no  change  in  the 
blockage  while  the  probe  was  transversed.   Two  tips  were  necessary  in 
order  to  cover  a  full  diameter  within  the  constraints  of  the  aerofoil 
movement . 

The  calibration  was  designed  to  obtain  optimum  accuracy  from  the 
measurements.   Furthermore,  the  expression  for  the  flow  rate  was  rederived 


8U 


in  order  to  include  the  effect  of  compressibility  explicitly,  rather  than 
as  an  empirical  correction  as  is  done  in  reference  5. 
B.l.  Analysis 

It  is  desirable  to  measure  the  pressure  drop  through  a  nozzle  as  a 
differential  water  column  measurement  because  the  drop  is  usually  small 
compared  to  the  absolute  level  of  pressure.  An  expression  for  the  flow 
through  the  nozzle  can  be  written  to  include  the  effect  of  compressibility, 
explicitly  in  terms  of  this  pressure  difference.  The  method  makes  use  of 
equations  which  were  derived  in  Appendix  A. 

Observe  first  that  the  upstream  flange  tap  registers  an  average 

stagnation  pressure  for  the  approaching  flow,  p   ,  and  the  throat  taps 

register  static  pressure,  p   .  The  probe  registers  the  local  stagnation 

s 

pressure,  p   ,  and  the  stagnation  temperature  is  assumed  to  be  constant 
at  a  value  measured  on  the  centerline  ahead  of  the  compressor.  Thermal 
expansion  of  the  nozzle  is  not  considered  since  the  entering  air  is  always 
at  close  to  ambient  temperature. 

Using  Eq.  A(19)>  for  any  streamtube,  in  consistent  units, 


*-A^[i -*£*)]•« 


t 

If  we  use  the  upstream  flange  and  throat  taps ,  then  the  coefficient  of 
discharge  is  defined  by  the  expression 


72pf  (pf-p.)    r      ,  /pf-ps\i  ,  ; 

-^-wr   I1  ■  &  (V)]  •  Ao  B(2) 


■t        ^T  N  *t 


which  is  used  to  obtain  the  flow  rate  when  C  is  known.  A_   is  the  area 
of  the  nozzle  throat.   If  R  is  the  distance  of  the  probe  tip  from  the 
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centerline  along  a  radial  traverse  and  R^  is  the  radius  of  the  nozzle 
integrating  Eq.  B(l)  and  equating  to  w  in  Eq.  B(2),  in  any  cousistent 
units, 


R 


C  = 


/ 


0  0 


PPVPS) 

Pf(pf-Ps.> 


0-1 


,p  -p 
3_  rv  *s 


)] 


L1"^ 


3  /Pf"Ps- 


.  R.dR, 


B(3) 


where  we  assume  p_  and  p   to  be  constant  during  a  traverse  of  the 

f       s 

calibration  probe.   In  fact  a  direct  point  by  point  evaluation  of  the 
above  integrand  minimizes  the  effect  of  slow  variations  in  the  flow  during 
the  traverse.   This  was  done  accurately  by  measuring  the  probe  pressure 
with  reference  to  the  flange  pressure  (average  stagnation  pressure). 

The  two  pressure  differences  that  were  measured  in  inches  of  water 
column  are  defined  in  terms  of  the  absolute  pressures,  which  were  measured 
in  inches  of  mercury  column,  by  the  equations 


and 


A  = 


A 


(Pf-P  ) 


(Pf-Ps) 


.07355 


.07355 


■B(k) 


It  is  observed  that  A   is  small  compared  to  p   or  p   and  does  not 

exceed  the  value  of  A 

s 

The  expression  for  the  discharge  coefficient  given  in  Eq.  B(3)  can  be 


expanded  in  terms  of 


and  —  to  give 
Pf       Pf 


C  = 


R 


§-J°I  (Pf,Ap,As)  R.dR 
0  0 


B(5) 
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where  the  function  F  is  given  by 


A_         A    0.07355- A  2n| 

_£L  J.  £_ 


F  =  [l  -  ^  -  ( 


Af     As  *  Pf 


A  A           A 

|1  +  0.03935  •  -£  ( l  -  .03^20  -S.  +  .07355  -£)                     B(6) 

L              P.p  \  .P.p            f  ~~ 

with  A   and  A   in  inches  of  water  and  p.  in  inches  of  mercury 

p       s  f 

column. 

When     C     is  known  as  a  function  of  Reynolds  number  from  calibration 

tests,  the  compressor  mass  flow  rate  is  then  given  by 


w  =  3-3135 


•CV^tT    *     i1-0'021^   •£)  B(7) 


"t0 


where  T,    is  the  stagnation  temperature  measured  in   F  and  where 

t0 
w  is  in  lbs.  per.  sec.  Eq.  B(7)  was  obtained  by  inserting  the  actual 

flow  nozzle  area  into  Eq.  B(2)  and  adjusting  the  units \     p   is  in  inches 

of  mercury  and  A   is  in  inches  of  water  column. 

B.2.   Calibration  Method  and  Data  Reduction 

At  a  fixed  compressor  condition  the  calibration  probe  was  traversed 

between  preset  scale  limits  which  placed  each  probe  tip  at  the  wall.  2k 

points  were  taken,  with  a  spacing  which  varied  from  0.010"  in  the 

boundary  layer  to  .25"  in  the  core.  At  each  point  the  probe  position  was 

read  against  a  vernier  scale  (S  inches)  with  the  centerline  position, 

S.   ,  determined  during  the  probe  installation.  At  each  point,  the 

heights  of  U-tube  water  columns  were  recorded.  The  connections  were  as 


87 


indicated  in  Figure  B  1.   If  H   ,  H   .  .  .  Hn  are  the  readings  of 
the  columns  in  cms.  of  water,  then  the  pressure  differences  are  given  by 


Ap  =  (Hq-H^/2.51* 


A 


A 


(H6-H5)/2.5^ 


(H3-Hu)/2.5^ 


B(8) 


and 


Pf  =  pa  +  0.07355  (H1-H2)/2.51+ 


B(9) 


where  p   is  the  atmospheric  pressure  in  inches  of  mercury,  and  the 
subscripts  1  and  2  refer  to  the  two  probe  tips. 

The  radial  positions  of  the  probe  tips  at  each  data  point  are  given 
by 


\       S0+^S12-S 
R~ 


R. 


and 


R2   S0"*5S12"S 


R 


0 


0 


B(10) 


where  S    is  the  distance  between  the  probe  tips. 

A  computer  program,  "CAL,"  was  written  to  accept  the  primary  data 


(S,  H  ,  Hp...Ho)   for  N  data  points,  with  p   , 


12 


0 


date,  test  number,  compressor  r.p.ra.  and  throttle  setting  entered  for 
the  test.   Equations  B(8),  B(9)j  B(l0) ,  B(6)  and  B(5)  were  programmed  and 
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the  integration  in  Eq.  B(5)  was  carried  out  using  an  overlapping  parabola 
technique  for  unequal  intervals  which  is  described  as  Function  Subprogram 
DAT INT. 

Program  CAL  computes  and  prints  out  a  coefficient  of  discharge  from 
the  data  for  each  probe  tip  separately,  C   and  C   ,  and  also  the 
average  of  C,   and  C   .  From  Eq.  B(7)  the  weight  flow  for  the  test  is 
calculated  using  the  average  coefficient  of  discharge,  and  the  average  of 
all  values  of  A^  and  p   taken  during  the  traverse.   The  Reynolds 
number  based  on  the  orifice  diameter  is  also  calculated  and  printed  out. 
A  listing  of  program  CAL  is  given  in  Appendix  G.2. 
B.3.  Results  of  the  Nozzle  Calibration 

Probe  surveys  of  the  flow  nozzle  were  carried  out  at  compressor 
speeds  of  8,000  r.p.m.  and  10,000  r.p.m.  A  throttle  position  for 
operation  near  to  surge  and  full  open  throttle  were  selected  at  each 
speed,  and  two  intermediate  throttle  settings  were  used  at  the  higher 
speed.  The  values  of  the  flow  coefficient  obtained  from  the  data  using 
program  CAL  are  listed  in  the  following  table;  and  the  average  value  is 
plotted  as  a  function  of  the  Reynolds  number  in  Figure  B  2. 


RPM 

ORIFICE 
SETTING 

Cl 

C2 

CA 

Average 

w 

REYNOLDS  NO. 

8,000 

k2 

.9897 

.9832 

.986U 

1.^937 

350^57 

10,000 

12 

.9901 

.9835 

.9868 

1.78U7 

^20605 

10,000 

15 

.9883 

.9865 

.987^ 

2.159^ 

505530 

10,000 

19 

.9880 

.9873 

.9877 

2.5829 

606825 

8,000 

16 

.9887 

.9883 

.9885 

2.9181 

686883 

10,000 

Qk 

.9900 

.9881 

.9891 

3.5211 

8320J+1 
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Since  the  stagnation  temperature  is  ambient  and  compressibility 
effects  are  very  small,  the  Reynolds  number  is  proportional  to  the  flow 
rate.   Figure  B  2  shows  that  the  discharge  coefficient  can  be  represented 
very  closely  over  the  range  of  the  data  by  the  expression 

C  =  0.98^5  +  0.0013  w  B(ll) 

where  w  is  in  lbs/sec. 

The  flow  rate  is  then  given  by  combining  Eq.  B(ll)  with  Eq.  B(7).   If 
we  write 

C  =  C    +   0.0013  w 

then 

*  "  r — .wo'oi3  .1  B(12) 

where  w'   is  given  by  Eq.  B(7)  when  C  =  C*  =  .98U5   . 

Also  shown  in  Figure  B  1  are  results  from  later  tests  carried  out  to 
check  the  above  calibration.  A  single  pitot-static  probe  was  used  to  trans- 
verse the  exit  of  the  nozzle  with  the  inlet  geometry  as  shown  in  Figure  5(c), 
The  method,  which  is  described  in  reference  6,  again  gave  two  values  of 
discharge  coefficient  which  were  then  averaged.  Figure  B  1  shows  that  the 
later  calibration  agreed  well  with  the  earlier  one,  however  there  were 
greater  differences  between  the  values  of  C  and  C?  and  less  consistency 
between  results  at  two  different  speeds.  The  original  calibration  given 
in  Eq.  B(ll)  was  therefore  not  changed. 
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APPENDIX  C.   THE  MEASUREMENT  OF  TORQUE 

C.l.   The  Initial  Flexure  System 

The  compressor  test  rig,  to  which  the  Hybrid  compressor  was  adapted, 
was  designed  to  have  an  internal,  two  flexure  torque  balance.  This 
method  was  adopted  because  physical  constraints  prevented  either  the 
attachment  of  an  external  arm  on  the  transonic  compressor  (for  which  the 
rig  was  originally  designed)  or  the  use  of  a  commercial  torque  meter.   It 
was  natural  therefore  to  use  the  same  system  for  the  Hybrid  compressor, 
and  merely  to  change  the  flexures  themselves  for  the  range  of  torque  to  be 
measured.   It  is  instructive  to  record  the  difficulty  encountered  with 
this  system. 

A  drawing  of  the  original  flexure  arrangement  is  shown  in  Figure  C  1. 
The  flexures  themselves  could  be  removed.  When  inserted,  they  seated 
against  counterbored  surfaces  in  a  cage  attached  to  the  stator  assembly. 
The  cage  was  free  to  rotate  on  two  ball  bearing  races.   The  inner  end  of 
the  flexure  was  located  at  a  slot  attached  rigidly  to  the  body  of  the 
machine.  Each  flexure  had  four  strain  gauges  in  a  full-bridge  arrangement, 
so  that  in  principle  two  independent  measurements  of  torque  were 
available.   The  two  bridge  outputs  after  signal  conditioning,  appeared  as 
channel  10  and  channel  11  on  the  data  system.   The  channels  were  calibrated 
statically  by  applying  weights  to  a  20"  arm  attached  externally  to  the 
stator  cage.   It  was  found  during  the  early  use  of  the  system  that  the 
calibration  curves  for  "loading"  and  "unloading"  weights  were  not 
separately  repeatable,  but  that  the  sum  of  the  readings  of  channel  10  and 
channel  11  had  considerably  less  scatter  over  the  range  of  torque  to  be 
measured.  Figure  C  2  shows  the  results  of  a  calibration  with  the  straight 
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line  used  initially  in  the  date  reduction  program.   The  line  fits  the 
calibration  data  to  +5$  over  the  range  of  interest.   However  the  departure 
of  the  torque  measurement  of  horsepower  during  a  test  from  the  predictions 
of  two  other  methods  (see  Figure  15  and  Section  3 . U )  was  several  times  this 
amount . 

The  torque  measurement  was  then  examined  in  some  detail.   It  was  found 
from  successive  calibration  measurements  that  the  outputs  of  the  two  channels 
were  sensitive  to  every  detail  of  the  position  and  retention  of  the 
flexures.  For  example,  an  increase  in  the  tension  of  the  retaining  set 
screw  was  found  to  change  significantly  the  output  of  each  channel 
separately.   Figure  C  3  shows  examples  of  calibration  test  results  which 
illustrate  this  sensitivity.  The  curves  for  the  two  channels  are  seen  to 
be  composed  of  straight  line  segments.  However  the  slopes  of  these  lines 
and  their  junctions  depend  on  how  much  of  the  load  is  being  taken  by  the 
flexure  in  question.  Also,  from  the  number  of  separate  slopes  required  to 
fit  the  data,  it  is  probable  that  the  support  or  contact  points  between  the 
flexure  and  mounting  hole  and  slot,  change  with  the  application  of  the  load. 
During  operation,  because  of  thermal  stresses  and  the  distribution  of 
mechanical  stresses  in  the  machine,  significant  changes  in  the  calibration 
of  the  flexures  might  be  expected  to  occur. 

Figure  C  k   illustrates  a  further  difficulty  experienced  with  the 
initial  flexures.  A  recording  of  the  output  of  each  channel  following  the 
loading  of  15  lbs  onto  the  balance  arm  revealed  a  steadily  decreasing  out- 
put.  Channel  11  decreased  by  1%   in  the  four  minutes  following  the  appli- 
cation of  the  load  and  had  not  reached  an  assymptotic  value  after  as  many 
as  50  minutes.  No  reason  for  this  behavior  was  found  although  it  was 
thought  to  be  related  to  the  method  of  mounting  the  flexures. 
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Improvements  to  the  two -flexure  arrangement  that  might  have  overcome 
these  difficulties  were  not  made  because  it  was  chosen  to  design  a  new 
system  to  be  totally  outside  the  machine.   The  desirable  improvements  in 
the  initial  system  would  include  increasing  the  clearance  between  the 
flexures  and  the  holes  in  which  they  were  mounted,  providing  a  cylindrical 
or  spherical  surface  on  the  inner  end  which  contacts  the  side  of  the  slot, 
increasing  the  clearance  in  the  slot,  and  improving  the  method  of 
retaining  the  flexure  against  its  seat. 
C.2.  Design  of  a  Single -Beam  Torque  Balance 

In  order  to  obtain  a  linear  and  repeatable  calibration,  the  torque 
measuring  system  shown  in  Figure  C  5  was  designed  with  the  intent  of 
eliminating  all  redundancy  in  the  transmission  of  loads.  A  single 
continuous  aluminum  beam  is  cantilevered  from  a  frame  which  is  bolted  to 
the  turbine  drive  unit.  The  beam  is  instrumented  with  four  strain  gauges 
in  a  full -bridge  arrangement.  The  free  end  of  the  beam  passes  through  a 
cut-out  in  the  20"  steel  arm  previously  used  to  attach  weights  for  the 
torque  calibration.  The  arm  transmits  a  force  to  the  beam  through  a 
1/V'  diameter  steel  ball  bearing.  A  welded  frame  serves  to  restrain  the 
arm  from  excessive  movement. 

The  beam  was  designed  to  measure  390  in.  lbs  of  torque  but  was 
required  to  withstand  twice  this  amount.  A  maximum  stress  in  the  beam  of 
5,000  psi  and  the  requirement  of  a  deflection  of  less  than  .1"  at  the 
design  point,  resulted  in  a  beam  having  a  section  1.5"  wide  and  0.5"  deep 
when  length  of  10"  was  fixed  by  mounting  requirements.  The  use  of  Alco 
7075  T651  alloy  resulted  in  a  beam  having  a  maximum  endurance  load  of 
150  lbs  from  a  maximum  torque  of  1800  in.  lbs. 

The  calibration  of  the  system  carried  out  by  loading  and  unloading 
weights  on  the  20"  balance  arm  showed  a  hysteresis  of  about  3  in.  lbs. 
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The  flexure  itself  was  then  loaded  directly  with  the  balance  arm  removed. 
From  0  to  50  lbs  load,  the  maximum  departure  from  a  single  straight  line 
was  0.06$>  of  reading.  The  analysis  of  the  following  section  was  carried 
out  to  account  for  bearing  friction,  stator  weight  and  residual  torque 
due  to  out -of -balance  during  calibration  and  operation  of  the  system. 
C.3.  Analysis  of  the  Torque  Measuring  System 


The  above  sketch  shows  the  forces  acting  vertically  on  the  stator  assembly 
during  calibration  or  operation  of  the  torque  system.  The  two  bearings 
are  denoted  by  Bn   and  B   and  the  radii  of  the  bearing  races  by  R^ 
and  R^   .  Distance  from  the  arm  is  given  by  D,   and  D   for  the 
bearings  and  by  D   for  the  C.  G.  of  the  stator  and  arm  assembly.  The 
normal  force  exerted  on  the  stator  by  the  arm  is  N  ,  and  by  the  bearings 
on  the  stator  is  N.   and  Np   .  N   is  positive  downwards  and  N   is 
positive  upwards.  The  weight  of  the  stator  and  arm  is  denoted  by  W 

The  following  sketch  shows  the  vertical  forces  and  torque  acting  on 
the  arm. 
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F  is  the  force  exerted  "by  the  flexure,  W  by  the  applied  weights  and  N 

is  the  net  force  from  the  bearings  and  stator  weight.  R   and  R^  are 

the  displacements  of  F  and  W  from  the  machine  centerline.  The 

residual  torque  due  to  stator  out -of -balance  is  t   and  the  bearing 

e 

torques  due  to  friction  are  given  by  p.-.N-.R^  and  (iJL  ,  and  act  in 
a  direction  opposite  to  the  tendency  to  rotate,  p,  and  ji_  are  bearing 
friction  coefficients. 

In  the  following  analysis  the   -+•   subscript  refers  to  increasing 
W  ("loading")  and  the   —  subscript  refers  to  decreasing  W 
("unloading"). 
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Equilibrium  of  forces  in  the  first  sketch  requires  that 


m2  +  ws  (d2-ds) 

Nl  =      [D-^p 


c(D 


and 


ND,  +  Ws  (Dl-Ds) 
N2  =  [D-^5 


G(2) 


Equilibrium  in  the  second  sketch  requires  that 


N  =  F+  -  W 


C(3) 


where  F  is  the  force  on  the  beam  while  loading  weights ,  and  F  is  the 
force  while  unloading  weights.  Equilibrium  of  the  torque  components  while 
loading  and  unloading  requires  that 


Wfy  =  F+RF  +  OX     +  ^X      -  t 
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Using  Eq.  C(3)  to  express  N  in  Eq.  C(l)  and  Eq.  C(2),  and  using  Eq.  C(l) 
and  Eq.  C(2)  in  Eq.  C(k) ,  there  results 
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We  define  two  "equivalent  frictional  radii" 


C(6) 


"A,    (W   +  ^V    ( W 


*£!  = 


VD2 


so  that  the  relationship  between  the  force  on  the  beam  and  the  applied  weight 
during  calibration  is  given  by 


V,  [l  i  £]  i  WgEs   -A^;[l|21 


Kf 


Ew- 


% 


To  first  order  in  the  small  quantity     —     , 

F 


C(7) 


^.T^lT 


'iVM1*^**; 


+  wsRs  +  \ 


0(8) 


The  following  sketch  shows  exageratedly  the  effect  on  the  calibration  of 
the  presence  of  friction  and  out -of -balance  torque  which  is  predicted  by 
Eq.  C(8). 
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During  the  operation  of  the  compressor,  a  pure  torque  replaces  the  weight 
W  at  radius  R^  .  The  frictional  correction  is  then  changed  because  the 
normal  force  on  the  hearings  is  reduced.   Instead  of  Eq.  C(3),  we  have 


N  =  F, 


C(9) 


and  instead  of  Eq.  C(k)9 


M  =  F+RF  i  ^NR       ±  ,2NR        -  T{ 


C(10) 


With  Eq.  C(l),  Eq.  C(2)  and  Eq.  C(3),  Eq.  C(10)  becomes  using  the 
definitions  in  Eq.  C(6), 


M  =  F+RF  t1  ±  B±  Vb  "T 


c(n) 


100 


However,  during  operation  the  applied  torque  will  vary  at  random  about  a 
mean  value,  and  the  vibration  will  act  to  eliminate  the  effect  of  friction. 
Taking  the  average  of  the  "torque  increasing"  and  "torque  decreasing" 
expressions,  the  torque  during  operation  is  given  by 


M  =  FRF  -  re  C(12) 


To  emphasize  the  implication  of  Eq.  C(8)  and  Eq.  C(ll);  during 
operation,  the  effect  of  friction  disappears  and  only  the  out -of -balance 
torque  is  present.  However,  if  a  calibration  of  the  beam  strain  gauge 
balance  is  carried  out  by  hanging  weights  on  the  arm  then  the  effect  of 
bearing  friction  should  appear,  as  in  the  above  sketch.  The  calibration 
of  the  strain  g|U^ge  balance  as  a  function  of  load  should  therefore  be 
carried  out  by  hanging  weights  directly  onto  the  beam.  This  gives  an 
expression 

F  =  F  (C)  C(13) 

where  C  is  the  number  of  counts  output  for  an  applied  load,  F  .  Then 
the  curves  shown  in  the  above  sketch  should  be  established  by  very  care- 
fully loading  and  unloading  weights  on  the  arm  with  the  system  fully 
assembled.  From  the  intercepts  of  the  loading  and  unloading  curves,  the 

out -of -balance  torque,  t   ,  can  be  obtained  as  shown  in  the  sketch.  The 

e 

torque  is  then  given  by 


M  =  F  (C)R_  -  x  C(lU) 
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C.U.   Calibration  of  the  Single  Beam  System 

The  result  of  a  direct  calibration  of  the  beam  alone  was 


F  =  5§J  0(15) 


where  F  is  in  lbs  and  C  is  the  number  of  counts.  The  maximum  deviation 
from  Eq.  C(l5)  was  2  counts  over  a  range  of  0  to  6,000  counts. 

The  results  of  the  second  stage  of  the  calibration  were  consistent 
with  the  predictions  of  the  above  analysis.  Figure  C6  shows  these  results 
over  successively  smaller  ranges.  Figure  C(6a)shows  the  loading  and 
unloading  curves  to  an  equivalent  full  scale  torque  of  600  in.  lbs.   The 
difference  in  slope  of  the  curves  is  small  but  discernible.   Departure 
from  linearity  is  not  detectible.  Figure  C(6b)shows  the  lower  15$  of  the 
range  and  Figure  C(6c)shows  the  intercept.  From  Figure  C(6c),  the  out-of- 
balance  torque  was  equivalent  to  Ik   counts.  Since  IL,  was  measured  to  be 
II.96  ins.,  the  final  calibration  was 


m  =  (0-110  •  jgl 


or 


M  =  0.09967  (C-lU)  C(l6) 


Finally,  as  a  further  check  of  the  analysis  the  difference  in  the 
slopes  of  the  two  curves  in  Figure  6  was  used  to  obtain  a  value  for  the 
friction  coefficient  of  the  bearings  to  compare  with  published  coefficients. 
The  difference  in  the  slopes  of  the  loading  and  unloading  curves  is 
estimated  to  be  O.M+Mi  counts/lb.  Therefore 
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But,  from  Eq.  C(8), 
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ing  the  definition  of  R^  in  Eq.  c(6)  and  writing 


C(18) 


following  reference  7  >  where  k  is  a  constant  with  the  dimensions  of 
the  radius  of  the  hall  hearings  r  ,  we  obtain  from  Eq.  C(l7)j 


k  = 


0.007^07 


(DrD2) 
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For  the  Hybrid  geometry, 


D   =  7.^8  inches         R_  =  11. 96  inches 
D   =  3.85  inches         R^  =  20.01  inches 


R^  =2.8  inches  r  =  l/8  inch 

R   =1.0  inch  r     =  l/k   inch 

B2  2 


From  Eq.  C(l9)  we  obtain  a  value  of  k  =  O.OO69  inches   .  Reference  7 
gives  a  range  of  k  for  well  oiled  surfaces  of  k  =  0.001  to  0.002 
The  high  experimental  value  is  thought  to  be  partly  due  to  the  difficulty 
in  resolving  the  difference  of  the  slopes  of  the  curves,  and  partly  the 
result  of  additional  real  constraints  on  the  stator  from  instrumentation 
leads . 

Since  instrumentation  leads  were  rearranged  frequently,  and  because  of 
the  importance  of  the  torque  measurement  in  specifying  efficiency,  the 
torque  measurement  was  calibrated  by  hanging  weights  both  before  and  after 
each  test.  The  slope  and  intercept  in  Eq.  C(l6)  were  obtained  by  carefully 
establishing  the  "load  increasing"  and  "load  decreasing"  curves  of  Figure 
C  6  over  the  range  of  operation,  and  calculating  the  mean  line.  An  accuracy 
of  the  order  of  0.1$,  was  thereby  maintained  in  the  specification  of  the 
torque . 
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FIGURE  C.l.   Full  Scale  Drawing  of  the  Original  Flexures  for 
the  Measurement  of  Torque  (see  also  Figure  l) 
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SUM  OF  OUTPUTS  (IN  COUNTS)  FROM  STRAIN 
GAUGES  ON  TWO  FLEXURES  (CH.IO  +  CH.Il) 


FIGURE  C.2.   Calibration  of  the  Initial  Torque  Flexures 
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HYBRID  COMPRESSOR  TORQUE  FLEXURE  CALIBRATION 
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FIGURE  C.3-      Calibration  of  the   Torque   Flexures    Individually 
Showing  the  Effect   of  Set  Screw   Pressure 
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FIGURE  C.4.   Stability  of  the  Original  Calibration  over  Time 
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END  VIEW  (WITH  SECTION  A-A) 


BALANCE 

ARM 

(SYMMETRICAL] 


CALIBRATION 
WEIGHT  PAN 


SAFETY 
FRAME 


SIDE  VIEW  (WITH  SECTION  B-B) 


MACHINE 
FRAME 


ALUMINUM 
FLEXURE 
WITH 
STRAIN 

GAUGES 


X 


J 


i  >. 

i  * 
i' 
i  ■ 
1 1 


:h 


jj 


j 


■> 


FIGURE  C.5.  Single  Beam  Torque  Balance 
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Figure  C.6(a).  Calibration  of  the  Single  Beam  torque  Balance 
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FIGURE  C.6(b).   Calibration  of  the  Single  Beam  Torque  Balance 
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FIGURE  C.6(c).  Calibration  of  the  Single  Beam  Torque  Balance 
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APPENDIX  D.   THE  MEASUREMENT  OF  TEMPERATURES 


The  stagnation  temperature  of  the  flow  entering  the  compressor  was 
measured  using  United  Sensor  Corporation  probes.   The  exit  flow  was 
measured  using  hare  wire  thermocouples  imbedded  in  the  struts  supporting 
the  honeycomb  flow  straightener  (Figure  8).  Check  measurements  in  the 
exit  flow  using  the  same  probes  as  in  the  entrance  invariably  gave 
readings  measurably  higher  than  the  fixed  probes.   Furthermore,  a  lack  of 
steadiness  was  observed  in  the  readings  of  all  probes,  and  two  probes 
mounted  close  together  in  the  flow  approaching  the  compressor  did  not 
give  the  same  reading.  Because  of  these  uncertainties,  all  probes  that 
could  be  removed  without  disconnecting  wires  were  subjected  to  a  static 
comparative  calibration  against  a  secondary  standard,  and  then  measure- 
ments were  made  of  the  temperature  profiles  into  and  out  of  the 
compressor.   This  was  carried  out  after  the  work  reported  in  Section  3. 3. 
and  before  that  reported  in  Section  3.^. 
D.l.  Static  Calibration  of  the  Temperature  Probes 

Five  probes,  identified  by  their  data  system  channel  number,  were 
calibrated  in  three  separate  static  tests .  A  quartz  thermometer  with  an 
accuracy  of  .01  C  was  used  as  the  reference  standard.  The  probes  were 
of  two  types: 

CHANNEL  NO.    PROBE  TYPE 


"\ 


21 

25 
26 

27 
28 


United  Sensor  Corporation  Type  TD 


B 


United  Sensor  Corporation  Type  TC 
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The  three  tests  were  as  follows: 

1)  Room  Temperature  Enclosure.  The  probes  were  taped  together  with 
the  reference  sensor  and  placed  in  a  thick  styrofoam  enclosure 
at  room  temperature  for  several  hours  before  data  was  taken. 

2)  Ice  Bath.   Probes  were  inserted  in  glass  tubes  containing  10-10 
oil  in  an  ice  bath.   The  reference  sensor  was  totally  immersed 
in  the  bath,  a  mercury  thermometer  was  inserted  and  the  bath  was 
stirred.  The  reference  temperature  of  31.83  F  agreed  with  the 
mercury  thermometer. 

3)  Oil  Bath.  The  probes  in  glass  tubes  and  reference  sensor  were 
immersed  in  an  oil  bath  heated  from  below  and  stirred.   Insula- 
tion was  added  until  the  bath  stablized  with  a  very  low  heating 
rate. 

The  results  of  these  tests  are  shown  in  Figure  Dl  in  which  the  deviation 
of  the  channel  output  from  the  reference  probe  output  is  shown  for  the 
three  tests.  In  this  plot,  1,000  counts  =  1  millivolt,  and  the  reference 
temperature  was  converted  into  counts  using  a  computer  subprogram  TEMPIC 
which  fits  the  ISA  standard  tables  for  iron-cons tantan  thermocouples  to 
the  accuracy  of  the  tables.  A  listing  of  TEMPIC  appears  in  Appendix  Gl. 
The  conclusion  drawn  from  Figure  Dl  was  that  the  reference  junction 
box  was  in  error  by  approximately  1.5  F  (or  k6   counts),  and  that  four  of 
the  probes  were  in  good  agreement  over  the  range  of  temperatures  from 
32  F  to  150  F.   Probe  28  appeared  to  be  intermittent.   The  application  of 
a  correction  of  +  k6   counts  to  account  for  reference  junction  error  leaves 
the  following  limits  of  accuracy: 
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PROBE 

ICE  PT      ROOM  TEMP       ~  150°F 

21 

25 

5  Counts 

+  5  Counts 

26 

+  15  Counts 

27 

28   - 

-  To  within  the  above  limits  when  not  intermittent. 

The  higher  deviation  in  the  oil  bath  results  was  thought  to  be  due  to 
temperature  gradients  in  the  bath,  since  the  probes  were  immersed  further 
than  was  the  reference  sensor. 
D.2.  Temperature  Measurements  in  the  Entering  Flow 

The  unsteadiness  in  the  readings  of  all  thermocouples  in  the 
compressor  intake  and  exit  flows  was  traced  to  large  fluctuations 
occurring  in  the  atmospheric  stagnation  temperature.  The  coastal 
location  of  the  laboratory  is  responsible,  and  hence  the  amplitude  and 
nuisance  of  the  fluctuations  are  a  function  of  the  date  and  time  of  day. 
An  example  is  shown  in  Figure  D2  which  shows  the  output  of  one  thermo- 
couple placed  in  the  free  atmosphere  close  to  the  compressor  intake. 
Readings  were  taken  by  continuously  rescanning  only  channel  21  on  the 
data  system,  giving  one  point  approximately  every  2  l/2  seconds.  The 
temperature  of  the  air  entering  the  throttle  is  seen  to  fluctuate  +  2  F 
about  an  average  value.  Because  of  the  thermal  inertia  of  the  metal 
plates  in  the  settling  chamber  and  of  the  intake  pipe  itself,  the 
fluctuations  are  damped  but  not  eliminated  before  the  air  enters  the 
compressor.  The  lower  frequency  fluctuations  remain,  so  that  variations 
of  greater  than  1  F  can  occur  in  the  inlet  temperature  channels .  For 
this  reason,  the  temperature  at  the  inlet  and  exit  of  the  compressor  were 


115 


recorded  on  adjacent  channels  of  the  data  system  and  a  direct  measurement 
of  the  difference  between  outputs  of  the  inlet  and  exit  probes  was  added. 

Temperature  measurements  of  the  flow  into  the  compressor  were  made 
in  the  19  l/2"  I.D.  pipe  and  in  the  7"  approach  pipe  at  10,000  r.p.m.  and 
with  full  open  throttle.  For  these  surveys  (using  channel  numbers  to 
identify  probes),  21  was  the  average  output  of  2  probes,  held  fixed  in 
the  larger  pipe,  27  was  a  single  probe  held  fixed  in  the  larger  pipe,  and 
26  was  the  survey  probe  which  was  read  with  respect  to  27  as  reference. 
The  results  of  three  radial  surveys  plotted  against  the  normalized  radial 
displacement  are  shown  in  Figure  D3.   In  viewing  this  figure,  note  that 
the  surveys  were  made  along  different  radial  lines,  as  shown  in  the  sketch, 
but  have  been  plotted  together  for  convenience.  The  bars  indicate  the 
spread  in  six  readings  taken  approximately  30  seconds  apart,  about  an 
average  value  represented  by  the  symbol.  The  figure  shows  that  the  air  is 
heated  by  as  much  as  2  F  in  the  settling  chamber  and  induction  pipe.  The 
wall  at  the  survey  station  is  seen  to  be  hotter  than  the  flow,  since  the 
compressor  exhausts  into  the  test  cell. 

Several  conclusions  were  drawn  from  these  data.  First,  the  fixed 
probes  (21)  measured  the  average  temperature  of  the  inlet  flow,  indicated 
by  the  probe  survey,  to  within  l/2  F.   Second,  the  recovery  factor  of  the 
temperature  probe  was  not  measurably  different  at  15*+  ft/sec  than  at 
23  ft/sec.  Finally,  the  specification  of  the  inlet  temperature  from  the 
output  of  a  number  of  probes  in  the  inlet  flow  cannot  be  done  to  better 
than  +  l/2  F  without  resort  to  a  heat  exchanger  or  thermal  capacitor. 
D.3.  Temperature  Measurements  in  the  Exiting  Flow 

Preliminary  measurements  of  temperature  were  made  in  the  flow  leaving 
the  honeycomb  flow  straightener  using  a  hand  held  probe  (26)  on  the 
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center line  at  various  circumferential  positions.  At  10,000  r.p.m.  and 
with  a  fully  open  throttle,  the  mean  of  the  readings  of  four  probes  in 
the  support  struts  of  the  honeycomb  (22)  was  almost  3  F  lower  than  the 
mean  of  the  readings  of  probe  26.   It  was  also  found  that  a  fixed  probe 
(28)  on  the  centerline  at  a  particular  circumferential  location  at  the 
exit  (station  5  in  Figure  l)  registered  the  average  of  the  distribution 
measured  by  26  very  closely.   These  measurements  were  repeated  for  a 
throttle  setting  close  to  design  at  12,000  r.p.m.  and  the  results  are 
shown  in  Figure  D^. 

In  viewing  Figure  DU,  it  should  be  noted  that  the  distribution  refers 
only  to  probe  26  which  was  moved.  The  readings  of  the  four  fixed  probes 
(22)  and  the  fixed  single  probe  28  are  shown  adjacent  to  the  corresponding 
reading  of  probe  26  for  convenience.  The  conclusion  was  that  the  average 
of  the  distribution  was  about  1.5  F  higher  than  the  average  of  the 
readings  of  the  fixed  probes.  That  probe  28  gave  a  reading  very  close  to 
the  average  of  the  distribution  was  also  confirmed  for  this  throttle 
position. 

Since  the  absence  of  scatter  in  successive  readings  on  the  fixed 
probes  (22  and  28)  compared  to  that  on  probe  26  suggested  a  peaked  profile, 
the  radial  distribution  was  measured  at  one  station.  The  result  in 
Figure  D5(b)  shows  that  the  profile  was  flat  and  was  not  the  cause  of  the 
scatter.  Figure  D5(a)  shows  a  partial  distribution  of  temperature 
increase  measured  radially  at  the  exit  of  the  rotor  with  the  face  of  the 
probe  (26)  set  at  the  design  flow  angle  of  65  to  the  machine  axis.  The 
profile  appears  to  be  more  peaked,  however  the  data  is  preliminary  in 
nature.  A  special  probe  is  required  to  resolve  the  distribution  of 
temperature  in  channels  which  are  0.6"  deep. 
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The  error  in  the  reading  of  the  fixed  probes  (22)  was  only  partly 
explained  by  considering  the  effect  of  compressibility  on  the  temperature 
recovery  factor  of  an  unshielded  thermocouple.   From  the  definition  of 
the  adiabatic  recovery  factor,  t  =  (T  -T)/(T  -T)   ,  where  T    is  the 
sensor  temperature  with  zero  heat  loss,  and  T  and  T   are  the  static 
and  stagnation  temperatures  respectively,  the  difference  between  the 
sensor  and  true  stagnation  temperature  is  given  by 

*r  "  VTaa  =  (1't)  (Vt> 

and  using  Eq.  A(l),  this  gives 


%  -   (1-)  i-  D(D 

p 


Taking  t  =  .85   ,  the  adiabatic  recovery  factor  for  a  flat  plate,  to  be 
as  small  a  value  of  t  as  could  be  expected,  b-     was  calculated  as  a 
function  of  the  weight  flow  rate  using  Eq.  D(l).  The  result  is  compared 
in  Figure  D6  with  the  test  data  for  the  deviation  of  the  fixed  probes  (22) 
readings  from  the  readings  of  probe  28.   It  was  concluded  that  the  error 
in  the  fixed  probes  was  larger  than  could  be  explained  by  compressibility 
alone,  and  conduction  loss  while  difficult  to  estimate  would  probably 
account  for  it. 

Following  the  above  tests,  the  four  unshielded  probes  were  replaced 
by  eight  shielded  probes  as  shown  in  Figure  19.   Tests  were  conducted  in 
which  the  probe  outputs  were  recorded  separately  as  a  function  of  time 
for  a  fixed  operating  condition.   The  time-averaged  outputs  were  then 
examined.   It  was  found  that  seven  probes,  which  were  nominally  identical, 
gave  a  spread  of  -  3J?fo   in  the  indicated  temperature  rise  at  all  throttle 
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settings.   (One  probe  was  found  to  be  broken).   Furthermore,  fluctuations 
in  the  temperature  rise  (as  in  Figure  D  5)>  originating  from  the  atmospheric 
temperature  variation  (Figure  D  2)  were  not  eliminated  by  recording  the 
rise  differentially.  Consequently,  as  can  be  seen  in  Figure  17,  the  degree 
of  agreement  between  the  different  measures  of  horsepower  (and  therefore 
efficiency)  could  be  adjusted  by  including  in  an  averaging  process  selected 
numbers  of  probes.  However,  the  fewer  the  number  of  probes  included,  the 
larger  was  the  uncertainty  resulting  from  time -dependent  fluctuations. 

It  was  concluded  that  the  disagreement  between  the  separate  measures 
of  power  could  easily  be  explained  in  view  of  the  difficulty  of  establishing 
an  accurate  value  for  the  temperature  rise  by  probe  measurements.   The 
torque  measurement  of  power  was  therefore  accepted  and  no  further  effort 
was  made  to  improve  the  measurement  of  temperature. 
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COUNTS  DEVIATION  FROM  THE  REFERENCE  TEMPERATURE 

FIGURE   D.l.      Results   of  Static  Calibrations    of  Five   Thermocouple   Probes 
sing  a  Platinum  Resistance   Thermometer  as   Reference 
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UNITED  SENSOR  CORP.   PROBE  TYPE  TD   (IRON-CONSTANTAN  THERMOCOUPLE) 
EXPOSED  TO  THE  ATMOSPHERE  OUTSIDE  THE  TURBOPROPULSION  LABORATORY 
AT  THE  S4VAL  POSTGRADUATE  SCHOOL. 
DATE        :  21  AUGUST  1972 
TIME        t  11*1*5  HRS. 

WEATHER:   OVERCAST  WITH  LIGHT  BREEZE 

FIGURE  D.2.      Fluctuation  of  the  Atmospheric   Temperature   in  the 

Neighborhood  of  the  Compressor  Test  Rig 
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FIGURE  D.3.      Temperature   Distributions   Measured  Ahead   of  the  Compressor 
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FIGURE  D.k.      Distribution  of  Temperature  on  the 
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FIGURE  D.5.  Radial  Distributions  of  Temperature  Downstream  of 

the  Rotor  and  at  the  Exit  of  the  Flow  Straightener 
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APPENDIX  E.  ANALYSIS  OF  THE  LABYRINTH  LEAKAGE 


The  above  schematic  shows  the  leakage  paths  through  labyrinth  seals. 
In  the  following  approximate  analysis,  the  notation  shown  above  is  used 
with  the  direction  of  the  arrows  showing  the  positive  direction  for  the 
mass  flow,  m  .  All  temperatures  are  stagnation  temperatures,  T   ,  and 
the  following  table  lists  the  use  of  subscripts: 
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MASS  FLOW  RATE      STAGNATION  TEMPERATURE  DESCRIPTION 

m   T   Flow  into  the  inlet  pipe 

mn  TQ  Flow  through  the  nozzle 

riL   T   Flow  through  the  rotor 

in       T   Flow  through  the  stator 

m   T   Leakage  through  outer 

labyrinth 

nv   T   Leakage  returned  ahead 

of  the  flow  nozzle 


m   T0  Leakage  returned  into 

c  d 

inlet  of  the  rotor 


ra   T   Leakage  into  the  stator 

through  rear  labyrinth 


It  is  assumed  that  the  gas  leaking  into  the  rotor  housing  from  the  rear 
labyrinth  enters  at  the  same  temperature  as  the  gas  entering  the  rotor. 
It  is  further  assumed  that  it  is  heated  during  passage  through  the  housing 
to  mix  at  the  exit  of  the  rotor  with  gas  at  the  same  stagnation  tempera- 
ture, T   .   It  is  assumed  that  the  leakage  through  the  inner  and  outer 
labyrinths  is  adiabatic  at  the  rotor  exit  temperature.  We  then  have  the 
following  relations  from  mass  conservation: 
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m.    : 


a         c 


%  =  &«  +  *b  E(l) 

n^  =  mQ  +  mc  E(2) 

E(3) 

n^  =  n^  +  m^  E(U) 

Conservation  of  energy  requires  that 

&0T0  =  *bT2  +  AcoT.  E(5) 

*1T1  =  *0T0  +  AcT2  E(6) 

fi^T     =  &1T1  +  ifa T     +   (Rate   of  Work  Input)  ~  E(7) 

P 

and  we  note  that  from  measurements ,  the  quantities  T   ,  T   ,  T   . 

CO         0         c. 

and  m..  are  known,  m   ,  m.  and  m   can  be  estimated  from  measure - 
0  a  '  £  c 

ments  of  pressures  and  the  labyrinth  geometry.  The  measurement  of  torque 
and  speed  of  the  compressor  is  a  direct  measurement  of  the  "Rate  of  Work 
Input."  From  Eq.  E(7), 

Rate  of  Work  Input  =  C  [nuT  -  (m+m )T  ]  , 


=  Cp  (m1+mi)[T2-T1]  using  Eq.  E(U)  , 

-C      (mn+m+mjrTo   -       °.°+ .c  2  1 
p       0     c     fL   2  m  +m         J 

using  Eq.    E(2)  and  Eq.    E(6)   , 
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(mn+m  +m  ) 

p   (m  +m  )  0J  L  2  0J 

0  c 


Cp"0  ^W 


1  +  t-4t- 

m  +m  . 

0  c 


E(8) 


Therefore,  since  the  left  hand  side  is  porportional  to  the  measured  horse- 
power based  on  the  torque  and  the  speed,  ^r-\      >  and-  C  m.  [T  -T  ]   is 
proportional  to  the  horsepower  based  on  compressor  flow  measurements  if 
the  leakage  is  ignored,  HP.,-   ,  it  follows  that 

r  \    1 

^Cl  =  HPC2  L1  +  m^Tj  E(9) 


The  efficiency  which  is  of  interest  is  one  that  evaluates  the  compressor 
from  the  inlet  face  of  the  rotor  to  the  exit  of  the  stator.   In  principle, 
this  eliminates  the  effects  of  leakages  through  the  inner  and  outer 
labyrinths  which  are  only  there  to  enable  the  torque  to  be  measured.   In 
actuality  the  presence  of  the  inner  labyrinth  will  cause  some  disturbance 
to  the  entering  flow.  The  leakage  through  the  rear  labyrinths  must  be 
included. 

Thus  we  are  concerned  with  the  same  rate  of  work  computed  above.  We 
consider  the  mass  flow  that  goes  through  the  rotor  and  is  raised  in 
enthalpy  from  temperature  T   to  temperature  T  =  Ti   ,  and  the  mass 
flow  that  leaks  in  through  the  rear  labyrinth  and  is  raised  in  enthalpy 
from  temperature  T,   to  T.   .  Hence, 

the  actual  rate  of  work  =  C   (m+m  )(T2~T  )  ,  as  given  by  Eq.  E(7)> 
and 
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the  ideal  rate  of  work  =  C  (m_+m  )  (T,   -T  ) 

is 

so  that  the  efficiency  of  interest  is  given  by  T]  ,  where 


Tl+.  "Tl 

Tj  =  -§V  E(l0) 

2  1 


m  T  +m  T„ 
0  0  c  2 

J.  I        ""       .    ~~i 

4 .      m  +m 
is 0  c 

2  "   V*c 


using  Eq.  E(6)  and  Eq.  E(2), 


IS    \      C  /   IS 


T  -T  /   m  V  T  -T 
2  0  '    0  N   2  0 


m 


where  7|,   is  the  efficiency  based  on  the  measurements  of  impact  pressure 
and  stagnation  temperature  at  station  0  and  station  k.     Finally, 

c  /    '1  , 


*  -  %  L1  -  ^  (— ;]  e(ii) 


so  that  the  correction  to  the  measured  efficiency  to  account  for  the 
labyrinth  leakage  is  at  least  as  low  as  second  order. 

We  would  like  to  find  the  effect  of  the  labyrinth  leakage  on  the 
efficiency  based  on  the  torque  measurement.  From  Eq.  E(8), 
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M  to  =  C  m_    (T0-Tj    (1  +    i    f>     ) 


0     c 


E(12) 


and  substituting  into  Eq.   E(l0), 


C  m_    /  &„ 

a  «  JL£  (i  +  ^JU 

1         M  to   \         m~+m 
0     c 


T4.  "To 

L_  IS 


where  Eq.  E(2)  and  Eq.  E(6)  have  been  used.  From  the  definition  of  TL 
the  efficiency  based  on  torque  measurements,  and  using  Eq.  E(l2),  we  have 


m. 


mn+m  y'3   Vm.+m 
0  c<     x  0  c 


m. 


1  =  i1  +  Azr\  ' 


E(13) 


or  rewriting 


l-Jj 


1  +  -i—  A    .  ,£  A 

An+rfL    V  m£   Tl3 


0     c 


E(lU) 


This  equation  is  interesting  in  that  it  says  that  the  correction  to  the 
efficiency  based  on  the  torque  and  speed  to  give  the  true  efficiency  of 
the  compressor,  depends  on  the  relative  rate  of  leakage  through  the  rear 
and  through  the  inner  labyrinths . 

There  is  a  relationship  between  the  efficiencies   TL   and  TL  gi 
by  equating  Eq.  E(ll)  and  Eq.  E(13): 


;iven 


V-K'-A),-!  ( 


m 


1  +  -r 


0  c 


m 


m^+m 
0  c 


Therefore, 
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*» 


m 


m_+m 
0  c 


-1 


m.+m 
0  c 
i1 


m      m 


\  +   m0+mc   m0 


E(15) 


or, 


\=(m4r) 


0  c 


m.   v      m      m 

i  +  *-L-\  -     °-  +  ,£ 


E(l6) 


We  now  rewrite  the  above  expressions  for  efficiency  using  the  following 
definitions : 


m 
c 

e   =  «r- 
c   m„ 


h  =  VAo 


E(17) 


j 


In  terms  of  the  quantities  defined  in  Eq.  E(17)  which  are  very  much 
smaller  than  unity,  the  Eqs .  E(ll),  E(lU),  E(15),  and  E(l6)  can  he 
expanded  to  give,  correct  to  second  order  in   e  , 


Tl 


=  \  ~    ec  i1  -   \) 


E(l8) 


=  ( 


i  -v  +  s  -  'i'cA  -  ^ 


-  e 


E(19) 


\  -   I1  +  ec  "  h  +   h)\   ■ 


E(20) 


\   "  l1  -  ec  +  h  ~   2  Vc 


+  £c2)^  + 


E(21) 


If  the  labyrinth  leak  rate  is  calculated  from  measurements  of  pressure 
and  the  geometry,  corrections  can  be  applied  to  TL   and  TL  using 
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Eq.  E(l8)  and  Eq.  E(19),  to  obtain  redundant  measurements  of  the  actual 
efficiency  of  interest. 

There  follows  an  approximate  analysis  of  the  flow  through  the 
labyrinths : 
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yy 
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Assume  that  the  labyrinth  behaves  as  a  series  of  square-edged 
orifice  plates  in  series.   Furthermore  it  is  assumed  that  the  stagnation 
pressure  downstream  of  one  tooth  of  the  labyrinth  is  equal  to  the  static 
pressure  at  the  gap.   The  latter  is  equivalent  to  assuming  a  sudden 
expansion  of  the  flow  to  infinite  area.  The  flow  is  considered  to  be 
adiabatic  and  with  the  low  pressure  drops  found  in  the  hybrid  labyrinths, 
the  flow  can  be  assumed  to  be  incompressible. 

We  can  use  the  expression  for  the  mass  flow  through  an  orifice 
derived  in  Appendix  B,  Eq.  B(2)  with  the  compressibility  correction 
omitted.   Thus  for  the  nth  tooth 


m  =  0.5385  CA. 


K   (Pn  -  Pn+l) 
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where  the  pressures  are  in  inches  of  mercury  absolute,  the  temperature  T 
is  in  Rankine  and  the  area,  A  =  nDt,  is  in  square  inches.   Taking  the 
coefficient  of  discharge,  C   ,  to  have  a  value  C  =  0.95  , 


or,  on  squaring  and  rearranging, 


A  =  1.607  -2-r  •  Vp.  /i>  -  P^  E(22) 


n  l^n   *n+l 


/   m   \     T 
Vl  =  pn  "  Vl^OTDt )     •   5"  E(23) 


The  method  of  obtaining  m  using  Eq.  E(23)  is  an  iterative  one. 
First,  assuming  a  value  for  m   ,  Eq.  E(23)  is  applied  across  each  tooth 
of  the  labyrinth  in  succession  until  a  value  of  p     is  computed.   This 
is  compared  with  the  measured  pressure  and  depending  on  the  difference,  a 
correction  is  applied  to  the  assumed  value  of  m   .   This  is  repeated  until 
the  error  in  the  final  pressure  is  less  than  the  inaccuracy  in  the  pressure 
measurement . 

A  card  program  was  written  for  the  Monroe  Calculator  1655  (designated 
Program  207)  to  perform  the  above  iteration  for  an  arbitrary  labyrinth 
(convergence  was  easily  obtained  in  less  than  two  minutes  for  each  point). 
The  leakage  rates  were  calculated  for  the  fully  open  and  near -surge  throttle 
conditions  for  preliminary  data  at  12,000  r.p.m.  and  with  more  recent  data 
at  15,000  r.p.m.   The  results  were  used  in  Eq.  E(l8)  and  Eq.  E(19)  to 
determine  the  corrections  to  the  measured  efficiencies.   These  data  are 
given  in  the  following  table,  from  which  it  may  be  concluded  that  in  the 
range  of  interest  for  the  efficiency  (away  from  wide  open  throttle)  the 
effect  of  labyrinth  leakage  is  negligible . 
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CALCULATED  EFFECT  OF  LABYRINTH  LEAKAGE 


DATA: 


N 
D(ins.) 
t(ins.) 


INNER  LABYRINTH 


9 
.01 


REAR  LABYRINTHS 
(OUTER)    (INNER) 


.006 


5 

1.5 
.006 


TEST  RESULTS 

SPEED  (RPM) 
RUN  NO. 
POINT  NO. 


12000 

23 
1  13 


33 


15000 
37 


k6 


INNER  LABYRINTH: 

P1(ins.  Hg) 

Pn  +  -^ins.  Hg) 

T(°F) 


m  (lbs/sec) 
c v   ' 


29.69 

25.  Ok 

2U.8U 

22.99 

27.52 

2^.02 

22.52 

21.29 

120.1+6 

123.00 

1^8.62 

151.65 

.01123 

.00710 

.01030 

.00852 

REAR  LABYRINTHS: 

P1(ins.   Hg) 

Pn  +  1(ins.   Hg) 

T(°F) 


m  (lbs/sec) 


\ 

S 
\  -  \ 

\  -  \ 


31.12 

29.85 

29.85 

27.95 

120. k6 

71.05 

.00725 

. 0090U 

.00300 

. 00302 

.0019H 

.00385 

.V73 

.732 

.i+52 

.7^3 

.00158 

-.00081 

.00387 

- . 00016 

.00229 

-.00065 

30.01 

30.01 

28.32 

27.22 

67.18 

70.  ok 

00860 

.01093 

00237 

. 002  in 

00198 

.00309 

.7^8 

.760 

.727 

.730 

00060 

-.00058 

OOO93 

-.00016 

00033 

- . 000U2 
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APPENDIX  F.  ANALYSIS  OF  COMPRESSOR  PERFORMANCE  TO 
ACCOUNT  FOR  NON -UNIFORM  EXIT  FLOW 


In  engine  component  testing,  a  variety  of  methods  have  been  devised 
to  account  for  non-uniformity  in  flow  profiles.  In  many  cases,  the 
different  methods  yield  very  similar  values  for  efficiencies.   This  was 
not  true  in  the  case  of  the  Hybrid  compressor  because  the  flow  velocities 
at  the  exit  measuring  plane  were  large,  and  because  the  outlet  geometry 
did  not  allow  a  very  complete  experimental  definition  of  the  profile.   The 
following  paragraphs  describe  methods  which  were  applied  at  various  stages 
of  the  test  program.   The  discussion  is  limited  to  accounting  for  non- 
uniformities  in  stagnation  pressure.   The  outlet  stagnation  temperature  is 
assumed  to  be  constant  throughout  the  flow- -a  reasonable  assumption  if  the 
actual  horsepower  is  derived  from  the  measurements  of  torque  and  speed 
rather  than  temperature  rise  and  flow  rate.  Also,  since  the  ratio  of  the 
depth  of  the  annular  passage  to  the  mean  diameter  is  small  (0.0555)  the 
radial  component  of  velocity  is  taken  to  be  zero.  Consistent  with  this 
assumption,  the  variation  in  static  pressure  radially  from  hub  to  tip  is 
estimated  to  be  less  than  0.07$>  for  exit  flow  angles  of  21  to  the  axis. 
The  static  pressure  is  therefore  assumed  to  be  constant  at  the  value 
measured  on  the  outer  wall. 
F.l.  Efficiency  Based  on  Weighted  Averages  of  the  Impact  Pressure 

The  simplest  approach,  though  not  an  analytically  correct  one,  is  to 
define  a  suitable  average  of  the  measured  impact  pressures  and  use  the 
expressions  for  efficiency  derived  for  uniform  flow. 

For  n  equally  spaced  impact  probe  measurements 


Va 
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defines  the  "area  averaged  impact  pressure,"  and 

(O  -^ F(2) 

l 

defines  the  "mass  averaged  impact  pressure."  In  terms  of  the  notation 
introduced  in  Appendix  A,  using  Eq.  A(l7)  and  the  equation  of  state  for  a 
perfect  gas,  Eq.  F(2)  becomes 


iKx* 


k  "4 


(%)  =^—   —  F<3) 


1 


\\ 


where  the  flow  function,  $    ,  at  each  probe  is  given  by 


(?)  A1 "  G 


§t.   Vp.  /  7L*   Vpt  > 


f(U) 


The  average  exit  impact  pressure  given  by  Eq.  F(l)  or  by  Eq.  F(3)  is 
used  in  Eq.  (10) ,  Eq.  (13) ,  and  Eq.  (15)  to  obtain  7]   ,  7)   ,  and  7U 
respectively. 
F.2.  Efficiency  Based  on  an  Average  Exit  Velocity 

With  the  assumptions  stated  above  the  flow  at  station  k   is  as  shown 
in  Figure  Fl.  Writing  Eq.  A(ll)  in  differential  form  for  a  streamline, 
using  the  notation  of  Appendix  A, 
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iim  .  J>      (_£*_) 


i_  f^u\  ii. . 


V  \ 


dA 


XU  I     2v     \         PU 


(w)  •  vT  cos  *\  ■  Mk  F(5) 


1  -  \2  '% 


The  integral  of  Eq.  F(5)  over  the  annulus  defines  a  relationship 
between  an  average  velocity  at  station  k,   which  is  consistent  with 
measured  static  pressure,  and  the  product  of  an  area  "blockage  factor" 
and  the  cosine  of  an  average  flow  angle. 
Thus, 


*/g  =  J  d(w/g)  -  (-^—)   (-^)  .  ^  •  %   cos  ak 


F(6) 


'k 


The  representation  separately  in  terms  of  blockage  factor  (which 
allows  for  separation)  and  flow  angle  (which  might  be  expected  to  be 
nearly  constant)  is  physically  realistic.   However,  without  a  knowledge  of 
the  profile,  the  product  itself  can  not  be  determined,  and  so  it  is 
sufficient  to  define  the  "blockage  factor,"  C.   ,  as 


C^  =  C'u  cos  o^  F(7) 


Then,  Eq.  F(6)  gives 


Xk2  +   ZB^  -1  =  0  F(8) 
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where 


BU  = 


fe*l) 


V.   (w/g) 


F(9) 


The  average  exit  stagnation  pressure  p   is  defined  in  terms  of  the 

% 
average  exit  velocity  X.  by  Eq.  A(6),  so  that 


|)-.(x-^)-{#. 


2BA 


-  1 


F(10) 


Using  Eq.  F(10)  in  Eq.  (15))  an  efficiency,  TV,  is  defined  thus 


*  9    Tt 

TV  =  0.89108 = 2. 

l(r  MN 


Y-l 
V    2B,X 


-  1 


!*■% 


F(ll) 


which  can  be  evaluated  in  terms  of  the  available  measurements,  using  Eq. 
F(8)  and  Eq.  F(9)>  if  a  value  is  assumed  for  the  blockage  factor,  C, 
F.3.  Efficiency  Based  on  an  Assumed  Form  for  the  Distribution  of  Impact 

Pressure 
F.3.1.   Ideal  Horsepower  as  an  Integral 

The  actual  horsepower  absorbed  by  the  compressor  is  given  accurately 
by  Eq.  (U)  if  the  tangential  velocity  is  everywhere  zero  ahead  of  the 
rotor  and  also  at  the  exit  of  the  flow  straightener  in  the  stator  section. 
However,  the  ideal  horsepower  is  given  by  Eq.  (9)  only  if  all  particles  of 
the  flow  arrive  at  the  same  final  state  (in  Figure  10).   If  they  do  not' 
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because   of  gradients   in  stagnation  pressure,  Eq.    (9)    is  replaced  by  an 
integral  over  the  mass  flow,  thus : 


y-1 

IHP 


1.4137  J        C  T       [( Jt)  '       -  l]   •    dw  F(12) 


With  the  assumptions  made  above  that  result  in  the  flow  description  shown 
in  Figure  Fl,  Eq.  F(5)  can  be  used  to  express  the  mass  flux  over  an  area 
element  at  station  k.     Eq.  F(5)  is  first  rewritten  using  Eq.  A(6)  and 
Eq.  A(l8)  as 

P     Pt 
«*/«>  =  (7^1)  y^  •  (  pj  ■  \)    '  cos  «k    •  dA^       F(l3) 

which  in  Eq.  F(l2)  gives 


pt,,     "  P.\f 


0  t4(i+)   4    4      0 


Also,  Eq.  F(13)  can  be  integrated  to  give  the  total  flow  rate,  which 
is  known  from  measurements :  thus 

\   (h)        *        4 

Equation  F(l4)  can  always  be  written  in  terms  of  an  average  flow 
angle , 


y-1 

P*  V4 


IHPC  =  1A1376  CpT      (^_)     "L.cob^.J-    (^.t)[(-|A)       -lK    F(l6) 

0  t^  (1+)  4         4  U 


lUO 


But,  if  it  is  then  assumed  that  the  variation  in  the  exit  flow  angle  is 
small,  cos  ov   is  given  also,  to  a  reasonable  approximation,  by  rewriting 
Eq.  F(15)  as 

*/«  =  (7^1)  v^  "  cos  %  *  I  ( 5T  §tkK  F(17) 

Substituting  for     cos  cvj,      in  Eq.   F(l6), 

,;,  §  ■  \)  k^  -  * 

IHP  =  1.1+1137  wC  T   -^ F(l8) 

C  "^0  ^t 

In  order  to  evaluate  the  integrals  in  Eq.  F(l8),  in  which  only- 
variations  in  impact  pressure  appear,  complete  surveys  must  be  carried  out 
at  station  k}   or  an  analytical  form  for  the  distribution  must  be  assumed 
in  which  free  parameters  can  be  evaluated  from  more  limited  measurements. 
F.3.2.  Distribution  Assumed  in  the  Peripheral  Direction 

The  following  method  was  a  first  attempt  to  evaluate  the  integrals  in 
Eq.  F(l8)  from  very  limited  data.  Four  impact  pressure  probes  at  station 
k   had  indicated  a  highly  distorted  profile.  Since  all  probes  were  on  the 
centerline  of  the  annulus,  two  located  peripherally  a  distance  b/U  from 
a  pressure  surface  trailing  edge  (type  A)  and  two  at  b/k     from  a  suction 
surface  trailing  edge  (type  B)  at  the  exits  of  different  blade  passages, 
the  distribution  shown  in  Figure  F2  was  the  most  general  that  could  be 
determined.  Here,  b  is  the  blade  spacing  measured  peripherally;   s  is 
the  normalized  distance  measured  peripherally  between  blades  and  h  is  the 
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normalized  distance  from  the  inner  surface  of  the  annulus.   It  is  assumed 
that  along  the  centerline  of  the  annulus ,  at  radius  R   ,  the  impact 
pressure  is  represented  by 


m- =  s(s  -  1)  (k  s  +  k  )  F(19) 


where  k,   and  k?  are  determined  by  one  probe  of  type  A  and  one  of 
type  B.   The  radial  distribution  of  impact  pressure  is  replaced  by  a 
uniform  condition  over  the  fraction  C'   of  the  annulus  height.   C   is 
then  a  blockage  factor  to  be  determined  by  continuity. 

From  one  pair  of  impact  probe  readings,  indicated  by  subscripts  A 
and  B   ,  Eq.  F(l9)  gives 


pU 


and 


.  ^  f  A  ..  i~l 


k2     3  Lvk 


T 


F(21) 


The  element  of  area  in  Eq.   F(l7)   and  Eq.   F(l8)    can  be  written  here  as 


C5 
dA.    =  2ttR      •  -^  •    ds  F(22) 

k  m        n 


where  n(=38)   is  the  number  of  blades.   The  integrands  are  then  functions 
of  s  only . 


1U2 


Writing 


I,,  =  — -  §, 


and 


V  Y 


Ju  -  ^  •  [(pr)    ■  XJ 


F(23) 


Eq.  F(l8)  becomes 


IHP  =  1.U137  wC  T. 

c  p  t 
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ds 


II 


F(2U) 


where  I  and  II  each  denotes  a  pair  of  impact  probes.  Eq.  F(l7)  gives 
for  the  product  of  the  area  blockage  factor  C'   and  the  average  flow  angle, 


*/g       .     V4 


C._  =  C'    cos 
2         P 


-   -  (%  -  1N\ A  

%  -  \  2Y  ;    P.  a,      ,.  1 


Vif 


F(25) 


[I  <V"  +  J  (VTTds 

0  I  0 


While  the  integrands  in  Eq.  F(2U)  and  Eq.  F(25)  are  analytic 
functions  of  s  ,  the  integrals  were  not  obtained  in  closed  form.   In 
applying  the  method,  for  each  pair  of  probes  I  and  II  ,  values  of  k 
and  k2  were  obtained  from  Eq.  F(20)  and  Eq.  F(2l).  From  Eq.  F(l9)  in 
Eq.  F(23),  (using  Eq.  F(U)),  the  integrands  I.   and  J.   were  calculated  at 
intervals  of  s  of  0.05  for  each  of  I  and  II  .  The  four  integrals 
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were  then  evaluated  using  subprogram  DATINT  which  is  described  in 
Appendix  B  of  reference  k   and  listed  in  Appendix  H  of  this  report. 
Finally,  the  total -to-total  efficiency  was  evaluated  from 


IHP   (Eq.  P(2U)) 
\   =        HP^  (Eq.  (k))  F<26) 


F.3«3.  Distribution  Assumed  in  the  Radial  Direction 

The  degree  of  uncertainty  in  the  efficiency  of  the  compressor  was 
clearly  a  function  of  the  distribution  of  stagnation  pressure  at  the  exit 
of  the  stator.  Twelve  holes  were  drilled  radially  to  permit  surveys  to 
be  made  at  6  different  stations  with  respect  to  the  blade  surfaces.  The 
first  surveys  showed  that  the  distribution  of  impact  pressure  assumed  in 
the  peripheral  direction  in  Section  F.3.2.  reasonably  represented  the 
actual  distribution  in  the  radial  direction.  By  inserting  12  Kiel  probes 
in  the  radial  survey  holes ,  which  were  distributed  almost  randomly  among 
the  38-"blade  passages,  the  equivalent  distribution  of  measurements  over  a 
single  blade  passage  (representing  an  average  condition)  shown  in  Figure 
F3  was  obtained.  The  impact  pressure  was  then  represented  by  the  cubic 
equation  in  the  radial  direction, 


Pt,  "  PU 

=  h  (h  -  1)  (k2h  +  k2)  F(27) 


PU 


and  was  taken  to  be  locally  constant  in  the  peripheral  direction  over  a 
fraction,  Cj-  of  the  quadrant  represented  by  a  pair  of  probes.  The 
blockage  factor  C^  was  again  determined  by  requiring  that  the  assumed 
distribution  satisfy  the  continuity  equation. 
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With  the  notation  shown  in  Figure  F3j  it  is  first  observed  that  an 
integral  of  any  function  of  the  impact  pressure  over  the  exit  can  be 
written  as 


J   (  -  )dAk   =  Ah   •  J1  [J1  (  -  )R  .  dh]ds  F(28) 

(k)  0   0 


where,   if     tu     is  the  blade  height, 


-        R  ^ 

R  =  r =  1  "  3T  (1  "  2h)  F(29) 

m  m 


If  the  function  in  the  integrand  goes  to  zero  when  p    is  zero,  then 


Eq.  F(28)  can  be  written  as 


J   (  -  )dAu  =  Cfa   £  £  [f(  -  )Rdh"j  F(30) 


W 


where  i  denotes  a  pair  of  probes  at  s  =  s . 


Using  the  definitions  given  in  Eq.  F(23)>  we  further  define 
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F(3D 


Jl+R  "  JhR 


and 


E 
m 

0 

E 
P 

■  J1  Jte  * 

0 

F(32) 


lk$ 


With  this  notation,  Eq.  F(l8)  becomes 
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F(33) 
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Using  Eq.    (4),  we  define  the   efficiency  based  on  the  assumed  impact 
pressure  distribution  as 


IHPC   (Eq.    F(33)) 
^6  =     HPcl   (Eq.    (I*)) 


F(3M 


or 


where 


6 

Is 


E 


m 


i=l 


F(35) 


H  = 


1.^137       r 
1.5865   '     p 


T     w 
t0 


■M   (n/ICT) 


F(36) 


The  blockage  factor  is  obtained  from  Eq.    F(l7),   giving 


Vs  V. 


k  6 


C^   =  c;   cos  a,    =  ft   "  1S)  : — -  •   -r 

6         6  4       V     2Y     /     p^  £ 


E 
L      m 

i=l 


F(37) 


1U6 


The  total -to -total  pressure  ratio  for  the  compressor  can  then  be  defined, 
using  Eq.  (9)>  as  that  pressure  ratio  that  would  be  obtained  for  the  same 
ideal  horsepower  if  the  flow  were  uniform.  Equating  Eq.  (9)  with  Eq.  F(33), 


F.(38) 


In  evaluating  TV   ,  Cs     and  -= —  from  the  probe  measurements,  first 
the  constants  in  Eq.  F(29)  were  calculated  for  each  pair  of  probes   (i) 
from  Eq.  F(20)  and  Eq.  F(2l).   I.    and  J.    were  calculated  from  Eq.  F(3l), 
using  Eq.  F(29)  and  F(23)  at  each  of  20  equal  increments  of  h  in  the 
distribution  given  by  Eq.  F(27).  E   and  E   were  then  obtained  for  each 
value  of  i  from  Eq.  F(32)  using  the  integration  subroutine  DATINT.   TV  , 
C,     and  p  /P   were  then  evaluated  using  Eq.  F(35),  Eq.  F(37)  and  Eq.  F(38) 
respectively. 


Xkl 


///// 


///// 


h 


777777 


s> 


9 


^VA 


STATION       (U) 


////// 


VELOCITY  IN  A 


RADIAL  PLANE 


VELOCITY  IN  A 


CIRC  INFERENTIAL 


PLANE 


FIGURE  F.l.   Velocity  Profiles  Following  the  Given  Assumptions 


At  S=5, 


BLADE  TRAILING 
EDGES  AT  STATION  k 
PERIPHERAL  SEP.-  b 


^f-^-  =   s(s-i)(k,s  +  lc,) 


FIGURE  F.2.   Peripheral  Distribution  of  Impact  Pressure 
Assumed  in  the  Analysis  of  Section  F  3-2 
lU8 


RADIAL  IMPACT  PRESSURE 


DISTRIBUTION  AT  S 


R* 

i 

t 

-■k- 

1 

■-♦-- 

r 

U- 

i  1 

1      L 

— 

r 

«3 

S4 

s 

BLADE  TRAILING  EDGES 
AT  STATION  k 


1.0 


.75  - 


.25  - 


=  h(h-i)(kfh+kt) 


POSITION  OF  12  KIEL  PROBES 
AT  STATION  h  SHOWN  WITH 
RESPECT  TO  A  SINGLE  STATOR 
BLADE  PASSAGE .(PROBES  WERE 
LOCATED  INDIVIDUALLY  IN 
DIFFERENT  PASSAGES  TAKEN 
AT  RANDOM  AROUND  THE  EXIT 
ANNULUS) 


ASSUMED  PERIPHERAL  DISTRIBUTION 


AT  ANT  VALUE  OF  h 


■i-C 


S  =  S;- 


ft.- ft 


FIGURE  F.3.     Radial  Distribution  of  Impact  Pressure 
Assumed  in  the  Analysis   of  Section  F  3.3 
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APPENDIX  G.  ANALYSIS  OF  MEASUREMENTS  OF  THE  FLOW  FROM  THE  ROTOR 

The  following  analysis  applies  to  probe  surveys  carried  out  at  station 
2,  the  exit  of  the  rotor  of  the  Hybrid  compressor.  However,  the  subscript 
2  has  been  omitted  everywhere . 

The  following  sketches  introduce  the  notation: 


PERIPHERAL 
DIRECTION  (u) 


'  AXIAL 

I  DIRECTION  (a) 


The  following  expressions  were  programmed: 


U  = 


TT 
3&0 


.  N  .  R 


V  =  V  .  cos  a 
a 


W  =  V  sin  a  -  U 
u 


W  =  /v  2  +  W  2 
-v/  a     u 


P  =  cos 


R  = 

R. 
1 

+  h' 

R 

R.    +  h' 

1 

R 
m 

R 
m 

h  = 

h' 

H 

G(D 


dA  =  2i#dR  =  2T^dh' 


Slip  factor  (p)  is  defined  as 


W 


U 


p,  *  ~  +   1 
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The  probe  measurements  can  be  reduced  to  profiles  of  mass  flux  and 

integrated  to  establish  a  blockage  factor: 

w  =  J  p  V  cos  a  dA  G(2) 

A 

Using  Eq.  A(2),  Eq.  A(l7)}  and  the  assumption  of  a  perfect  gas, 


«-(&*:j 


p    t 


cos  a   . 


dA 


*-'  "t  A"  l. 
Using  the  above  notation  and  introducing  conversion  factors 


w 


_  /gy\  32.172 
'  \y-lj    2.036    * 


pA 


V 


t 


'<&  0  • 


cos  a  . 


_R_ 

R 

m 


dh 


G(3) 
2 


where  the  static  pressure,  p,  is  in  inches  Hg.  absolute,  A  is  in  (inches) 
and  V,  is  in  feet  per  second. 

Defining  the  blockage  factor  by  referring  the  mass  flow  rate  to  the 
centerline  properties  (denoted  by  subscript  m), 


w  =  0  V  .  cos  a     .  k  A 
rm  m       mm 


G(h) 


or,  similarly  to  above, 


-(&>-•*■■ 


cos  a    [  —  $, 

p  t 


.  k 


m 


G(5) 


m 


Defining, 


Ik  = 


cos 


m 


G(6) 


and  using  Eq.  G(5)  in  Eq.  G(3)>  the  blockage  factor  is  given  by 


k  =  -~-       f1  .  I,  .  dh 
m   (^m  J 


G(7) 


0 


Equation  G(6)  was  programmed  with  the  relationships  in  Eq.  G(l)  as 
Card  Program  210  for  the  Monroe  Model  1655  calculator.   The  program  instruc- 
tions and  program  listing  follow: 
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MONROE  CARD  PROGRAM  210,  -  OPERATING  INSTRUCTIONS 


This  program  calculates  the  velocity  triangle  for  the  exit 
flow  (station  (2)  )  of  the  Hybrid  compressor  from  cobra  probe  data. 
A  mass  flux  integrand  is  also  calculated. 


Inputs  t  T  - 
N 

h«  - 

*2  " 

P2  " 

Pt2  * 


(°F) 

(rpm) 

(ins)  -  probe  position  above  the  hub. 

(°) 

(any  units)  -  wall  static  pressure (absolute) 
(same  as  p2)  -  impact  pressure (absolute) 


Outputs:  See  list  below. 

TO  LOAD  :  PRESS  (RESET ),(T0(0)),( LOAD) 
OPERATE  :  PRESS  (T0(0) ), (RESUME) . 


LOAD  k   CARDS 


-RELEASE  (LOAD). 


L- PRESS  (RESUME) - 

OR 

—PRESS  (T0(0))— ' 
FOR  NEW  Tt2 


ENTER  T 

■   N 

i 

"   h 


t2 


i 


't2 


RECORD 


THEN 


(0) 


(RESUME) 
(RESUME) 
(RESUME) 
(RESUME) 
(RESUME) 
(RESUME) 
h 


U. 


(1) 
(2) 


't2 


a  2 


W 


u2 


t(*) 

t(6) 

t(7) 

(6)  +  f(U)  +  T-~^ 

(5)  t  f(7)  -/sin^'Y2ndV^0V  -p* 

V  cos  y  ^FUNcy^  )        v 


(3)  (Ik) 


W. 
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APPENDIX  H.     COMPUTER  PROGRAM  FOR  THE  HYBRID  COMPRESSOR 

H.l     Notation 

There  follows  a  list  of  the  most  important  symbols  and  their  meaning 
when  used  in  the  on-line  reduction  program  given  in  Section  H.2,  in  the 
batch-process  reduction  program  given  in  H.3>  and  in  the  text  of  this 
report. 


FORTRAN  SYMBOL 

TEXT 

VARIABLE 

ARRAY 

SYMBOL 

Input  Data: 

RUN 

SRUN 

- 

PT 

SPT 

- 

DAY 

- 

- 

MONTH 

- 

- 

YEAR 

- 

- 

CHN(I,J) 

- 

- 

CH1(I) 

CN1 

- 

CH2(I) 

CN2 

- 

CH(J) 

CN 

- 

NDIFF 

- 

- 

GAMMl 

- 

Yo'i 

GAMM2 

- 

Y()'2 

SL0PE1 

- 

- 

CEPT1 

- 

- 

SL0PE2 

- 

- 

CEPT2 

_ 

_ 

UNITS 


DESCRIPTION 


(°) 
(°) 


Test  run  number 
Data  point  number 

^Date 

No.  of  counts  on  channel  I  scani- 
valve  port  number  J 

Data  value  on  port  I  of  channel  1 

Data  value  on  port  I  of  channel  2 

Data  value  on  channel  J 

=  1  for  Diffuser  A 
=  2  for  Diffuser  B 

First  diffuser  row  stagger  angle 

Second  diffuser  row  stagger  angle 


Slope  and  intercept  of 
*  torque  calibration 


Slope  and  intercept  of  flow 
'nozzle  AP  transducer 
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FORTRAN  SYMBOL 

TEXT 

VARIABLE          ARRAY 

SYMBOL 

Physical  Constants: 

GAMMA 

V 

A 

Ao 

\               '       " 

\ 

C 

c 

UNITS 


DESCRIPTION 


Ratio  of  specific  heats 

(ins)    Flow  area  at  station  (0) 

2 
(ins)    Flow  area  at  station  (k) 

Value  of  the  flow  nozzle  coeffi- 
cient at  zero  flow  rate 


Pressures : 
PA 
PI 

POC 


SPA 


POT 
PTUCl 
PTUC2 

PTUC3 

PTl4CU 

PTU(I) 

BW 


Pr 


'fl 


PIC 

SP1C 

Pl 

P2C 

SP2C 

P2 

Pi+C 

SP^C 

?k 

PT1C 

SPT1C 

Ptl 

PT2C 

SPT2C 

Pt2 

ptUc 

SPT*4C 

(ptU} 

m 


(ptA>I 

(ptB}I 

(ptA}II 

(ptB}II 

*W2 


ins .  Hg   Atmospheric 


Upstream  of  the  turbine 
flow  orifice 

Upstream  flange  of  the  compressor 
flow  nozzle 

Wall  at  station  (l) 

Wall  at  station  (2) 

Wall  at  station  (h) 

Stagnation  at  station  (l) 

Centerline  stagnation  from  cobra 
probe  at  station  (2) 

Mass  average  of  12  Kiel  probes  at 
station  (k) 

Stagnation  ahead  of  the  turbine 


Fixed  impact  probes  at 
*  station  (h) 


Ith  Kiel  probe  measurement 
ins.  HpO   Drop  across  turbine  flow  orifice 
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FORTRAN  SYMBOL 
VARIABLE    ARRAY 


SYMBOL 


UNITS 


DESCRIPTION 


Pressures  (con't) 
HWC 


As 


Drop  across  compressor  flow  nozzle 


Pressure  Ratios 


PRTST 
PRTTC 

PRTTPR 
PRT1C 

ptUr(i) 

PTT^R(l) 
PRA1C 

prUtc 


SPRTTS 
SPRCTT 

SPRCP 


Total-static  across  turbine 

Total-total  across  compressor 
(mass -averaged  p  ,  ) 

Total-total  across  compressor 
from  radial  profile  method 
(Eq.  F(k2)) 

Static-total  across  compressor 
from  station  (l)  to  station  (k) 

Ith  Kiel  probe  to  station 
(h)   static 

Ith  Kiel  probe  to  station 
(l)  stagnation 

Atmospheric  to  station  (l)  static 

Station  (k)   static  to  station  (l) 
stagnation 


Temper atur 

ss : 

Tl 

- 

5t 

°I 

TIN 

STINT 

T 
tTl 

TOUT 

STOUTT 

T 
XtT2 

TINC 

STINC 

Tto  ■ 

TOUTC 

S TOUTC 

TtU 

DELT1 

SDELTC 

- 

DELT2 

SDELTT 

_ 

Flow  approaching  turbine  flow 
orifice 

Flow  into  the  turbine 

Flow  out  of  the  turbine 

Flow  into  the  compressor 

Flow  at  station  (h) 

Rise  through  the  compressor 

Drop  through  the  turbine 
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FORTRAN  SYMBOL      TEXT     UNITS 
VARIABLE    ARRAY    SYMBOL 


DESCRIPTION 


Flow  Rate: 

WT1 

- 

wm' 

T 

WT 

SWTT 

T 

WC 

SWTC 

w 

lbs/sec   1st  approx.  for  the  turbine  flow 
Turbine  flow  rate 
Compressor  flow  rate 


Forces,  Speed  and  Power: 
FORCE 

TORQ  STORQ       M 

RPM  SRPM       N 

HPT  SRPT2      HP, 


HPC 


HPC2 


T 


SHPTQ     HPC1 


SHPC2     HP, 


lbsf 
in-lbs 
rpm 
h.p. 


C2 


Axial  force 

Compressor  stator  moment 

Rotational  speed 

Horsepower  to  the  turbine  from 
flow  measurements 

Horsepower  of  the  compressor  from 
torque  and  speed 

Horsepower  of  the  compressor  from 
flow  measurements 


Referred  Quantities: 
REDTQ     SREFTQ    M^ 


REDRPM 
REDHPT 
REDHP1 

REDHP2 

REDWC 


SREFN 


SREFHP 


SREFW 


N. 


REF 


HP. 


REF 


w_ 


REF 


in-lbs_   Referred  torque 

r.p.m.    Referred  speed 

h.p.     Referred  turbine  horsepower 

h.p.    Referred  compressor  horsepower 
based  on  torque  and  speed 

h.p.    Referred  compressor  horsepower 
based  on  flow  measurements 

lbs/sec   Referred  compressor  flow  rate 


Efficiencies : 


ETATO 
ETAO 


SETATS 
SETAO 


\ 


Total-static  turbine  efficiency 

Overall  compressor  efficiency 
defined  in  Eq.  (20) 
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FORTRAN  SYMBOL      TEXT     UNITS 
VARIABLE    ARRAY    SYMBOL 


DESCRIPTION 


Efficiencies  (con't): 
ETA1      SETA1      T^ 


ETA2 


SETA2 


\ 


ETA3 

SETA3 

h 

ETAU 

SETA3 

\ 

ETA5 

SETA  5 

\ 

ETA6 

SETA6 

\ 

ETAD 

SETAD 

\ 

Compressor  efficiency  defined  in 
Eq.  (10) 

Compressor  efficiency  defined  in 
Eq.  (13) 

Compressor  efficiency  defined  in 
Eq.  (15) 

Compressor  efficiency  defined  in 
Eq.  (16) 

Compressor  efficiency  defined  in 
Eq.  (17) 

Compressor  efficiency  defined  in 
Eq.  (18) 

Efficiency  of  the  stator  defined 
in  Eq.  (2l) 


Other  Variables 


RD 

SRFD 

h 

CCOSTH 

SC6 

c6 

CB 

SC5 

S 

Ck 

scU 

ck 

VZERO 

_ 

V, 

B 
VZEROU 


to 


B 


'tk 


B, 


Static  pressure  recovery  factor 
defined  in  Eq.  (23) 

Blockage  factor  from  Eq.  F(37) 

Blockage  factor  from  Eq.  F(25) 

Blockage  factor  from  Eq.  F(7) 

ft/sec   Limiting  velocity  at  entrance 
temperature 

Defined  in  Eq.  A(l3)  for  flow  at 
station  (0) 

ft/sec   Limiting  velocity  at  exit 
temperature 

Defined  in  Eq.  A(l3)  for  flow  at 
station  (k) 
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H.2  Online  Reduction  of  Performance  Test  Data  -  PROGRAM  TCRIG1 

Data  from  the  B  &  F  data  scanning  system  (Figure  9)  is  output  on  paper 
tape.   The  format  received  by  the  computer  via  the  teletype  input  is  a 
series  of  data  statements  that  are  equivalent  to  data  cards.   The  following 
shows  an  example  of  a  data  scan  output  on  the  teletypewriter: 


RUN  PT. 

NO.  NO.  DATE 

2}  16  2p087g  1ST  LINE  FORMAT 

ooi  oi  +02905  o        Remainder"  " 
ooi  02  +OU5U6  o 

/up  to  2U  port\ 
Inos.  ON  CH.l  J 

001  2k   +03000  0 

002  01  +02985  0 

/up  TO  2U  port  \ 
INOS.  ON  CH.2    J 

002  2k  +03U27  0 

010  00  +038U5  0 

Oil  00  *03Ul6  0 
■ 

/UP  TO  90\ 
yCHANNELSJ 

090  00  001201  0 

I        I 
CHANNEL     DATA 
NO.      VALUE 

SCANI-     (NOT 
VALVE       USED) 
PORT  NO. 


The  tape-computer  interface  is  described  in  references  8  and  9, 
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H.3  Batch  Process  Reduction  of  Performance  Test  Data  -  PROGRAM  TCRIG4 

Original  data  that  is  input  on-line  in  TCRIG1  via  paper  tape  is  output 

on  nine  cards  for  each  point  scanned.  Cards  from  up  to  50  points  can  be 

processed  using  TCRIGU  to  obtain  original  data  and  reduced  data  in  tabulated 

form  (and  also  computer  plots  from  library  subroutines.) 
A  listing  of  the  program  follows: 
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H.U  Computer  Program  for  the  Calculation  of  Flow  Nozzle  Coefficient  from 
Survey  Data  -  PROGRAM  CAL 
This  program  accepts  raw  data  from  a  two-probe  survey  of  the  Hybrid 
compressor  flow  nozzle  and  outputs  the  discharge  coefficient,  the  weight 
flow  and  the  Reynolds  number.   The  manometer  hook-up  is  shown  in  Figure 
Bl.  The  method  is  described  in  Appendix  B. 

Application  to  other  nozzles  would  require  the  following  changes: 
(i)  The  first  data  card  (statement  2) 
(ii)  The  nozzle  radius  (statements  15  and  2U) 
(iii)  Manometer  arrangement  (statements  7  to  10) 

There  follows  a  listing  of  the  notation,  followed  by  a  listing  of  the 
program: 


FORTRAN 
SYMBOL 


TEXT 
SYMBOL 


UNITS 


DESCRIPTION 


NDAY 

NMTH 

NYR 

NSCAN 

NRPM 

ORIFICE 

PATMOS 

TATMOS 

N 
SZERO 

S12 


N 


■A 


"to 


r .p.m. 
ins.  X  10 
ins .  Hg 


12 


inches 


inches 


DAY 

MONTH 

YEAR 

IDENTIFICATION  #  OF  SURVEY 

COMPRESSOR  SPEED 

Overlap  of  largest  holes  in  orifice  plate 

Atmospheric  pressure 

Flow  stagnation  temperature  (taken 
to  be  atmospheric) 

Number  of  survey  points 

Scale  reading  for  probes  at  equal 
radii 

Distance  between  probe  centers 
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FORTRAN 

TEXT 

SYMBOL 

SYMBOL 

S(I) 

S 

H1(I) 

Hl 

H2(I) 

H2 

H3(I) 

H3 

*(I) 

Kk 

H5(I) 

H5 

H6(I) 

H6 

H7(I) 

«7 

H8(I) 

H8 

PT(I) 

pf 

D1(I) 

*1 

D2(I) 

fiP2 

F1(I) 

•  F 

F2(I) 

F 

B1(I) 

Rl 

R2(I) 

R2 

FRl(l) 

F.R. 

FR2(l) 

F.R. 

UNITS  DESCRIPTION 


inches      Scale  reading  for  Ith  point 

cms.  HpO     Height  of  1st  manometer  column  for 
Ith  point 

Height  of  2nd  manometer  column  for 
Ith  point 

Height  of  3rd  manometer  column  for 
Ith  point 

Height  of  Uth  manometer  column  for 
Ith  point 

Height  of  5th  manometer  column  for 
Ith  point 

Height  of  6th  manometer  column  for 
Ith  point 

Height  of  7th  manometer  column  for 
Ith  point 

Height  of  8th  manometer  column  for 
Ith  point 

ins.  Hg     Pressure  at  upstream  flange  tap 

ins.  HpO     Difference  between  upstream  flange 
and  probe  1  pressure 

ins.  HpO     Difference  between  upstream  flange 
and  probe  2  pressure 

Function  defined  in  Eq.  B(6)  for 
probe  1 

Function  defined  in  Eq.  B(6)  for 
probe  2 

Normalised  radial  displacement  of 
probe  1 

Normalised  radial  displacement  of 
probe  2 

Integrand  in  Eq.  B(5)  for  probe  1 

Integrand  in  Eq.  B(5)  for  probe  2 
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FORTRAN 
SYMBOL 


TEXT 
SYMBOL 


UNITS 


DESCRIPTION 


AINT 

BINT 

CI 

C2 

C12 

WTC 
R 


C. 


'Average 

w       lbs/sec 


fFRdR  between  1  =  1  and  I  =  N  for 
probe  1 

jFRdR  between  1  =  1  and  I  =  N  for 
probe  2 

Orifice  flow  coefficient  from  probe  1 


Orifice  flow  coefficient  from  probe  2 

(C1  +  C2)/2 

Compressor  weight  flow 
Reynolds  Number 
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APPENDIX  I.      TABLES    OF  PERFORMANCE  AND  ROTOR  EXIT  FLOW 
DATA  AT   1^,000  RPM 
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TABLE  I-I.   TABLE  OF  DIFFUSER  PERFORMANCE  RESULTS 


The  table  is  in  six  sections: 


Section 

Run  No. 

Diffuser 

(i) 

h9 

A 

(ii) 

50 

B  1 

(iii) 

51 

-J 

(iv) 

52 

B 

(v) 

53 

A 

(vi) 

5U 

A 

same  configuration 


Diffuser  "blade  angles  are  given  in  the  table. 
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TABLE  I-II.   TABLE  OF  OVERALL  PERFORMANCE  RESULTS 


The  table  is  in  six  sections 


Section 

Run  No. 

Diffuser 

(i) 

^9 

A 

(ii) 

50 

B   1 

(iii) 

51 

*   i 

(iv) 

52 

B 

(v) 

53 

A 

(vi) 

5^ 

A 

same  configuration 


Diffuser  blade  angles  are  given  in  the  table. 
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TABLE  I-III.   TABLE  OF  MEASURED  DATA 

The  table  is  in  six  sections  as  follows: 
Section        Run  No.         Diffuser 


(i) 

h9 

A 

(ii) 

50 

B  ] 

(iii) 

51 

-  J 

(iv) 

52 

B 

(v) 

53 

A 

(vi) 

5h 

A 

same  configuration 


Diffuser  blade  angles  are  given  in  the  table. 


200 


I 

CJ— . 
UJ2- 
UJQ. 

o-o: 

OO  — 
I 

.2  CO 


Si  1 

az 

2.*-i 

w 

1 

O 

OoO 

J\ 

u_cr> 

l  — 

UJ 

or:  — 

on 

I—  • 
UJQ 
LU 


CJ   O   O   tJ   O   O   O   O   O   CJ   O   CJ   O   CJ   CJ   O 

a>ooo^c^o>Of>oocT>oooo>o 

c^r-cctvO(\jfrisrsto^"ino>4-oor*f-H 
j •••• 

K     f*     N     rt     in     K     w     i4     ^j     08     ifl     Of*     O     ^     W     SI 


vO 

—i 

CO 

m 

«4> 

1-1 

o 

>t 

r- 

o 

■4- 

CO 

CO 

UN. 

r- 

CM 

•o 

ro 

ro 

o 

-r 

o 

in 

•4- 

rH 

Is- 

«4- 

ro 

r-t 

cj> 

Is- 

>o 

co 

in 

•M 

Cvl 

00      CO 


on     oo 

ui     -4* 
Is-     r- 


^cOeOcocOcMcMCMCMCMrococMCMcor^ 


Is-       OvfrCM.-l^i-lO^^'CMCOCMOi-'CMOO 

h-     >rcocjrOv0^r-ooNp^oorfiooo-<vO 

3    •••••••••••••••    t 

ouors-is-NO>d'^"m^oo*otf>fs-mm>4-\C 

♦—      'Nj(Nj(M<M(\JCM(M(\|(\jrsjfvJ<\JfMC\j(\jCM 

c\ic\jo^vO>i-ooo>>tr*-^-icoinroomfM 
ZLnmooroa*'-tcof\ih-ois-r«.oro-4D     4* 

■nmiftmtnmmmmm^mtnmmm 


Z> 

a: 

a: 

o 
oo 

00 
UJ 
Cf 

a. 

s: 
a 
o 


at: 

CO 

> 
x 
i 
i 
i 

< 
i- 
< 
o 

o 

LU 

or 

00 

< 

UJ 

u. 

o 

UJ 

-I 
co 
< 


H 
H 
H 
I 
H 

3 


0* 

(M 

o 

r- 

in 

o 

00 

CO 

in 

—1 

^ 

m 

CM 

CM 

CO 

m 

•4- 

r\j 

o 

CNJ 

00 

>4- 

ON 

in 

CO 

00 

CO 

*o 

vO 

Is- 

Is- 

o 

o 

1— 

Q- 

>0 

LO 

ST 

CO 

co 

CM 

CM 

CM 

— i 

CM 

fO 

CO 

CM 

CM 

m 

LO 

ro 

CO 

CO 

CO 

CO 

rn 

CO 

CO 

CO 

CO 

CO 

co 

CO 

rO 

(O 

rO 

00 

O 

r-4 

o 

00 

•4- 

O 

co 

I— 1 

CM 

o 

i—i 

r-i 

CM 

CM 

CO 

fO 

UJ 

cv 

«* 

i-4 

cr> 

ro 

00 

CO 

o> 

r- 

«fr 

p- 

CM 

CM 

T 

r-\ 

1-1 

fH 

Qi  — 

1— 

Ijoo 

o_ 

o 

00 

«o 

»0 

in 

in 

•4- 

>4" 

«* 

>t 

•O 

nO 

in 

in 

00 

00 

OOCQ 

«* 

ro 

CO 

CO 

CO 

CO 

CO 

CO 

fO 

ro 

ro 

ro 

CO 

ro 

CO 

ro 

o0< 

LU    • 

QfO 

O.X 

LO 

o 

>o 

o 

in 

O 

r- 

t-i 

>* 

in 

CO 

-H 

CO 

O 

t-4 

o 

• 

I—I 

i-t 

r- 

ON 

s0 

CM 

ON 

m 

LO 

CO 

>r 

•4" 

»r 

P- 

Is- 

o 

o 

—J  00 

H 

<2 

a. 

<3 

«* 

ro 

CO 

co 

CM 

CM 

CM 

CM 

CM 

CO 

CO 

CM 

CM 

<f 

<t 

K»-h 

CNJ 

CNJ 

CM 

f\i 

rg 

CM 

CM 

CM 

CM 

CM 

CM 

r-4 

CM 

CM 

CM 

(NJ 

O  — 

h- 

Is- 

O 

Is- 

00 

00 

<> 

o 

*o 

h- 

h- 

r- 

O 

ON 

00 

Is- 

r- 

00 

00 

00 

00 

00 

00 

CO 

CO 

00 

00 

00 

00 

oo 

00 

00 

00 

<t 

O- 

o* 

CJN 

0 

ON 

ON 

0> 

ON 

ON 

ON 

o> 

ON 

0> 

ON 

ON 

ON 

ON 

<N 

CNI 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

1 

o 

00 

LO 

<f 

i-H 

r- 

in 

in 

O 

>t 

CM 

CM 

ON 

00 

O 

CJN 

o 

ON 

ON 

ON 

ON 

00 

00 

00 

00 

00 

ON 

ON 

00 

00 

O 

CJN 

■4- 

00 

Q- 

o 

0> 

o 

O 

ON 

0s 

o 

CJN 

ON 

0> 

ON 

CJN 

0> 

o> 

O 

ON 

LU 

CO 

CNJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

CM 

ac— 

300 

OOoQ 

00< 

•4- 

Is- 

O 

CO 

CM 

oo 

in 

m 

.-1 

co 

CO 

*-t 

CM 

»H 

CM 

CM 

LU    t 

•4- 

i-l 

O 

>* 

CM 

o 

cr 

ON 

CJN 

ON 

CO 

co 

O 

o 

f-l 

—I 

CH3 

rg 

Q-X 

0l 

0* 

00 

r- 

r- 

r- 

r- 

>o 

O 

>o 

nO 

Is- 

p- 

Is- 

p- 

00 

00 

• 

CNJ 

CNJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

OoO 

h*Z 

!-•-« 

<«- 

CM 

O 

in 

CO 

CM 

in 

00 

in 

>4" 

o 

Is- 

r- 

i—i 

00 

O 

ON 

r- 

ro 

o 

CO 

o 

Is- 

4- 

t-i 

i—i 

O 

1—1 

00 

00 

CO 

CM 

o> 

CO 

oo 

rH 

a. 

m 

«4- 

en 

c0 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ro 

CO 

CM 

CNJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

z.  • 

hQ 

a. 


CO 


in    ^o 


00         CJN 


o-*cMco>frmvors-co 

CMCMCMfMCMCMCMCMCM 


201 


I 

C_>  — . 

■ua. 
ujcx 
act 
co«— 

I 

t— ~. 

UJ-J 
S    I 

DZ 


OOCJOtJOOO'JCJOCJOO 
^■cOrO^cOrO^rOro^rO^frrOrci 

roocr>ror-4-mror\irMroco3j.-H 


I 

O 
3:lu 
OW 

LL.CO 

—I 

\  — 


CO 

UJ 

cr-» 

h-  • 
<o 

QLdJ 

ujq 
Q_~ 

'<- 

LU 


<N 

0> 

O 

r- 

in 

r- 

*0 

<f 

>o 

>r 

(N 

nO 

<»- 

in 

rM 

ro 

00 

o> 

ro 

>* 

CT> 

<-* 

CO 

o 

ro 

r~ 

-0 

•—I 

ro 

•H 

0> 

r- 

m 

ro 

O 

Ch 

r- 

r- 

nO 

m 

o 

r-H 

o 
in 


rOrorororo?sir\)r\jr\jc>gc\iro 


moooo^o^cM^tovOOvor— 
k  tncMcr»>rcr<rocor»'--*mcor-r\jm 
Zl       ••••••,.       t 

O<-*r0OOC7>OCTc0OC0C0J3OCT> 
h-rorOrorocMcorMCMCMCMCMCMcMCM 

■4"inro>0'-<mmoNOOroro  .— i  co 
zrvji^^oom^oorocoh-ocNjom 

»_      r- *oCT,oooooocooocOh-ooco<»oo 
-Ovommmmmmmmmmmm 


2: 
z> 
a: 

en 
o 

CO 

CO 

LU 

a: 

o_ 

s: 
o 


a: 
co 

> 
x 

I 

l 

I 
< 

< 

a 
a 

LU 

r> 

CO 

< 

LU 


O 

LU 

-I 
CO 

< 


m 

in 

co 

O 

r-» 

r- 

O 

in 

O 

P«- 

CO 

«t 

CO 

CO 

** 

CM 

>o 

CNJ 

CO 

CM 

r- 

f-t 

00 

O 

** 

ro 

rM 

ro 

i— ( 

H- 

O- 

vO 

in 

in 

>* 

•4" 

ro 

rO 

CM 

CM 

CM 

CM 

CM 

CM 

rO 

ro 

ro 

ro 

rO 

ro 

ro 

ro 

CO 

ro 

CO 

ro 

ro 

CO 

ro 

CO 

CJ 

CO 

i— i 

—* 

o 

CM 

CO 

t-l 

o 

1— 1 

i-h 

CO 

>t 

<o 

LU 

CM 

o 

0> 

CM 

O) 

CO 

O 

CO 

u\ 

p> 

CM 

i— i 

o 

r-t 

CO 

QC  — 

H- 

3co 

0. 

r- 1 

o 

0> 

CO 

r«- 

s0 

in 

in 

in 

in 

m 

in 

in 

in 

COCO 

-4" 

ro 

CO 

ro 

ro 

CO 

ro 

CO 

CO 

ro 

ro 

ro 

CO 

CO 

CO< 

UJ    • 

oca 

O.X 

<\J 

CO 

o 

r- 

CO 

r» 

>* 

o 

vO 

r- 

vO 

rO 

r- 

o 

• 

F-« 

0> 

r\J 

r* 

o 

•4- 

o 

in 

co 

r-4 

o 

o 

O 

o 

*o 

-J  CO 

K 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

t 

<Z 

o_ 

>o 

vO 

in 

in 

•4" 

>r 

rO 

CO 

CO 

CO 

ro 

rO 

ro 

co 

h-—> 

(NJ 

cm 

CM 

CM 

CM 

rM 

(N 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o~ 

h- 

h- 

>o 

o 

in 

>4" 

P* 

>o 

O 

r*~ 

r*- 

r- 

in 

in 

vO 

O 

O 

a 

o 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

< 

a 

o 

O 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

ro 

ro 

CO 

CO 

CO 

CO 

r0 

CO 

CO 

CO 

ro 

rO 

CO 

CO 

1 

CO 

CO 

CO 

r- 

-o 

o 

CO 

CO 

>o 

in 

i— i 

O* 

f-l 

00 

o> 

a> 

CT> 

o> 

o 

o 

CF 

0^ 

CT> 

o> 

0> 

00 

cr> 

0> 

«fr 

CO 

a. 

0> 

o> 

O 

0> 

o 

o 

0> 

o 

0> 

o> 

0> 

O 

0> 

0> 

LU 

CM 

CM 

CM 

CM 

CM 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ac— 

3co 

COCO 

CO< 

>t 

r^ 

00 

>* 

CM 

>fr 

ro 

O 

c- 

r- 

r- 

—* 

vO 

r- 

LU    • 

sO 

o 

in 

—> 

r- 

m 

ro 

CM 

CM 

CM 

CM 

ro 

CM 

ro 

C*0 

CM 

0-X 

a. 

0> 

0^ 

co 

00 

r- 

r- 

r- 

r»- 

r- 

r- 

h- 

r- 

r^ 

l»- 

• 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rM 

CM 

CM 

oco 

HZ 

»-•-• 

<>■» 

in 

>o 

CM 

o 

o 

O 

m 

O 

nO 

CJ^ 

.—1 

O 

O 

r-l 

1- 

o 

•4" 

O 

CO 

ao 

>T 

o 

CO 

P~ 

<? 

r* 

vO 

O 

-4 

co 

i-4 

Q. 

>0 

in 

>* 

<!• 

co 

co 

ro 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

CM 

cm 

CM 

CM 

CM 

CM 

r\l 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

4-   • 

oz 
o. 


-h     cm     m 


in 


oo     cr>     o     — i     cm     co     in 


202 


Q  — 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

l_> 

O 

JJZ 

*£ 

r-t 

vO 

o 

r- 

cr- 

u\ 

C\J 

cr« 

m 

—1 

r- 

CO 

LUO. 

o 

0> 

o 

o> 

on 

o 

O 

o 

i— • 

O 

cr> 

0s 

Q.CX 

-4- 

co 

vfr 

co 

rO 

4- 

sf 

CO 

st 

•4" 

CO 

CO 

CO'—' 

1 

1—  -» 

CO 

>*• 

CO 

CO 

CNJ 

CNJ 

o 

CO 

o 

ON 

r— 1 

CO 

zcu 

C* 

LU-J 

a 

in 

O 

CT- 

ir\ 

sO 

r-l 

CO 

in 

vO 

r- 

in 

o 

51  1 

u 

Is- 

.0 

ro 

i— t 

CO 

Is- 

m 

r-t 

cr- 

co 

r-t 

>£> 

QZ 

I— 

4" 

4- 

4- 

-4- 

cO 

ro 

co 

ro 

r\j 

CNJ 

CO 

4- 

£»-h 

"■* 

1 

CO 

cr> 

cr 

ro 

r-l 

-*• 

•4- 

o 

O 

fO 

in 

CNJ 

o 

CM 

ao 

co 

Is- 

i— 1 

r» 

00 

Is- 

O 

CO 

Is- 

r- 

JiLU 

m 

r-l 

cr 

Is- 

UN. 

ro 

o 

00 

r- 

o 

co 

r-l 

Oi/) 

3 

• 

• 

• 

• 

• 

• 

t 

• 

• 

t 

• 

• 

_JS 

>* 

«* 

rO 

co 

ff) 

co 

CO 

<\j 

CNJ 

<\J 

<M 

vf 

LLCD 

_J 

1  — 

co 

uj 

0d~* 
3U- 

t-     • 

aCOJ 

UJO 


cooco^-^-mOvO^co^oo 
H-Is-0s'CNjrHCNjc0*OCNJr-rs-inr-* 

3       •'..*■     t       t       •••••••• 

oooxOi-ii-ic^oo,'Oxoor^^ 

I—      cOC\|rOrOroc\|rOC\lc\|CNjc\lcO 

o"^"<-*ro>4"in^coNCcooo 

t—      oor-O'O'^r-O^cococOr^O 
mmtn«ovOinmmm,    mm>o 


in 


■z. 

CO 

sO 

CO 

CO 

•—I 

r-t 

0s 

in 

SO 

o 

ro 

CNJ 

z> 

<r 

cnj 

p- 

CNJ 

r- 

CNJ 

ao 

r-t 

Is- 

>t 

CO 

Is- 

Is- 

CC 

h- 

a. 

sO 

m 

in 

*• 

>t 

ro 

ro 

CNJ 

CNJ 

CNl 

CNJ 

m 

CC 

CO 

CO 

CO 

CO 

CO 

ro 

CO 

C<1 

CO 

CO 

CO 

ro 

u 

CO 

CO 

UJ 

CO 

•4- 

r-t 

h- 

u> 

o 

<* 

<7> 

<f 

Is- 

m 

00 

Is- 

a: 

UJ 

ini 

r-t 

CNJ 

CNJ 

CNJ 

CNJ 

Is- 

ao 

>r 

CNJ 

r— • 

ro 

Is- 

a. 

a:—. 

i— 

S 

X>co 

a. 

r-i 

O 

ON 

CO 

r- 

O 

in 

in 

in 

in 

in 

0N 

O 

coco 

•* 

>fr 

ro 

ro 

co 

ro 

ro 

co 

m 

ro 

ro 

ro 

o 

co<r 
UJ    • 

o 

al(j) 

t— < 

o-x 

CO 

r- 

cr> 

r- 

CNJ 

Is- 

O 

o 

co 

ro 

r-t 

CNJ 

CL 

• 

.— t 

CO 

CO 

-O 

o 

*t 

o 

in 

CNJ 

o 

o 

ro 

CO 

ao 

—ICO 

H 

• 

# 

t 

• 

• 

• 

• 

• 

• 

• 

• 

• 

> 

<;z 

Q. 

•O 

o 

s\ 

in 

•4" 

4- 

ro 

.-o 

CvJ 

ro 

CO 

o 

X 

1 

1 
■ 

1— 

r\j 

CNJ 

CNJ 

CNJ 

(NJ 

CNJ 

CNJ 

C\! 

CNJ 

CNJ 

CNJ 

CNl 

< 

CO 

s0 

sO 

r~ 

CO 

ON 

ao 

ON 

o 

o 

O 

o 

1- 

o 

o 

o 

a 

o 

o 

O 

o 

r-t 

r-t 

i-i 

r-l 

< 

< 

a 

o_ 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

a 

ro 

CO 

CO 

ro 

ro 

CO 

ro 

ro 

CO 

CO 

CO 

ro 

Q 

UJ 

CC 

1 

ID 

CO 

r- 

CO 

0> 

o 

a- 

r-t 

ON 

sO 

o 

o 

CNJ 

CO 

o 

0> 

0> 

0> 

o 

cr 

O 

0* 

<> 

ON 

o 

o 

<t 

4< 

• 

• 

• 

• 

• 

• 

• 

■ 

• 

f 

• 

• 

UJ 

CO 

o_ 

0s 

& 

0> 

0^ 

o 

0> 

O 

0> 

cr> 

ON 

o 

o 

SL 

UJ 

or— 

(\J 

CNJ 

CNJ 

CNJ 

ro 

CNJ 

CO 

CNJ 

CNJ 

CNJ 

ro 

ro 

LL 

=3  CO 

o 

CO  CO 

co< 

CO 

CNJ 

GO 

CNJ 

CNJ 

ro 

in 

ON 

in 

1— 1 

O 

00 

UJ 

UJ    • 

in 

r-t 

in 

r-t 

r» 

in 

co 

CNJ 

CNJ 

CO 

CO 

O 

_J 

QCO 

N 

t 

• 

• 

• 

• 

• 

• 

• 

t 

• 

• 

• 

CO 

O-X 

Q. 

o> 

0> 

CO 

00 

r>- 

Is- 

Is- 

Is- 

Is- 

Is- 

r~ 

ON 

< 

t 

cm 

CNJ 

CNJ 

CNJ 

cni 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

I- 

OcO 

1 

<-*z 

*-*- 

^-^ 

<*- 

co 

in 

ro 

^O 

0- 

O 

Is- 

CO 

o 

Is- 

ON 

O 

•H 

1- 

o 

in 

ON 

co 

r- 

in 

o 

00 

vO 

sO 

CO 

m 

•H 

CO 

•-4 

•H 

o. 

vO 

in 

-t 

>4- 

ro 

ro 

CO 

CNJ 

CNJ 

r\j 

CNJ 

m 

H 

CM 

CNJ 

CNJ 

CNJ 

CM 

CNJ 

CNJ 

CNJ 

CN) 

rvj 

CNJ 

CNJ 

H 

H 

I 

H 

1- 

H 

z  • 

OZ 
a. 

r-t 

CNJ 

ro 

<■ 

in 

>0 

Is- 

CO 

ON 

o 

r-l 

T-t 

CNJ 

r-l 

203 


I  ••• 

l_J-»  OOOCJCJOOOOOO<-> 

uji  Z     h     ui     ^     rti     Is     'O     h     rg     r»     h     ^     m 

U-IO-  OOOsO^^OOoO^OOO^ 

Q.a:  xt^rOrocOcO-^-J-rO^cOCO 

I 

Z.CO  0«         ••/•••••.••'•"• 

uj_j  a.r-<-4omcMrnf\isOUNr-,^0N<J 

5"  I  O00if>(r)t-400v0rfi^(?00(MO 


CM 

in 


a: 

a 

oo 

UJ 

a 
a 

s: 
o 


a: 

CD 
> 
X 

I 

I 

I 

< 
I— 
< 


a: 

oo 
< 

UJ 

2: 

LL 

o 

UJ 


0 

:<uj 

Ooo 

~s. 

-J-v 

U-GO 

_J 

1  ~ 

\- 

00 

~D 

UJ 

O 

a:-* 

1- 

23  U- 

H-    • 

<o 

CtlUJ 

UJO 

z 

a«- 

►— < 

s. 

1- 

UJ 

CO 

t— 1 

00 

0* 

sO 

CM 

CO 

0 

On 

o> 

co 

r~ 

0 

co 

CO 

in 

cO 

0 

0 

r- 

m 

en 

0 

c* 

rr\ 


CO 


^^m^m^^NN 


o> 


CM 


lf\ 


Is- 
on 


o 


in 


•0 
o 

t 


o* 

(O 


in 

• 
c\j 


O 
0> 


—I  CM 

0*  CO 

•  • 

co  cm 


no 


o 
in 


on 

on 

CO 

CO 

Is- 
co 

-0 

m 

en 

m 

<0 

co 

Is- 
co 

in 
co 

co 

0 
in 

CM 

0 

<0 

r- 

1— t 

in 
en 

00 

in 
o* 

0^ 

CM 

in 

0 
in 

0 

1— 1 

5. 

H- 

r- 

0 

co 

o> 

^0 

O 

•J- 

-4" 

•4- 

vO 

>r 

in 

UJ 

sO 

r- 

O 

vO 

O 

vO 

0 

O 

O 

nO 

^0 

0 

1— 

1 

0 

■O 

r-i 

GO 

in 

CM 

in 

O 

rO 

<r 

co 

s«- 

<T 

1— 1 

co 

o> 

*fr 

CO 

in 

0 

CO 

in 

CM 

Is- 

O 

(- 

Q. 

<o 

in 

>fr 

»t 

co 

CO 

CO 

CM 

CM 

CM 

%t 

CM 

m 

co 

CO 

rO 

CO 

CO 

co 

CO 

<o 

co 

CO 

CO 

00 

vO 

i-i 

0 

vO 

vO 

*£> 

CM 

u> 

in 

ON 

in 

0 

UJ 

CM 

in 

0 

CO 

vO 

>0 

1— 1 

in 

CM 

0 

CO 

0 

1— l 

ac- 

1— 

Doo 

a 

0 

0> 

CO 

Is- 

•0 

vO 

m 

in 

in 

>4- 

CO 

m 

00  CO 

>t 

CO 

CO 

rO 

CO 

CO 

co 

CO 

rO 

CO 

co 

CO 

00< 

UJ    • 

oiO 

O.X 

sf 

in 

vO 

Is- 

in 

>fr 

0 

CO 

r- 

CO 

vO 

0^ 

• 

»H 

O 

r- 

CO 

CO 

CM 

CJ^ 

in 

CO 

r-l 

0 

CM 

CM 

_ioo 

t— 

• 

# 

• 

• 

• 

• 

• 

• 

• 

• 

t 

• 

<Z 

Q. 

■O 

in 

in 

•4- 

>*- 

CO 

CO 

CO 

r*\ 

co 

in 

rO 

1—  HH 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

O-" 

H- 

O 

in 

r- » 

O 

O 

CO 

•4- 

CO 

•4" 

V* 

CM 

in 

0> 

o> 

0> 

0> 

O 

co 

O 

0> 

On 

ON 

0> 

cr> 

< 

• 

• 

• 

• 

t 

• 

• 

• 

• 

• 

• 

• 

a. 

O 

o> 

cr 

0* 

0> 

0> 

Cj> 

0* 

o> 

ON 

O 

ON 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

1 

in 

CO 

.—1 

vO 

«4- 

CM 

in 

^f 

f-l 

in 

— 1 

CM 

0 

0 

0 

CJ> 

CJ> 

O 

0> 

0> 

0^ 

CO 

O 

On 

>t 

• 

• 

• 

• 

• 

t 

• 

• 

• 

• 

• 

• 

00 

Q. 

0 

0 

0 

0> 

CT» 

c^ 

O 

ON 

ON 

0^ 

O 

ON 

UJ 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CNJ 

CO 

CM 

oi-» 

ID  00 

OO  CO 

O0< 

in 

CO 

0 

O 

CO 

in 

in 

co 

CM 

00 

in 

in 

UJ    • 

o> 

CM 

0 

r- 

CO 

CM 

r-» 

rH 

1— 1 

r-t 

ON 

t-» 

aCO 

CM 

O-X 

a 

00 

CO 

00 

r- 

P» 

r- 

r- 

P* 

Is- 

Is- 

Is- 

Is- 

• 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

(M 

CM 

CM 

CM 

CM 

Ooo 

-<z 

!-"• 

<w 

f-< 

CO 

CO 

CO 

CO 

>o 

0 

r* 

-O 

CM 

.—1 

in 

1— 

r- 

U> 

m 

—* 

■O 

CO 

0 

00 

r"- 

Is- 

in 

CO 

00 

1— 1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o_ 

in 

«* 

<f 

<f 

co 

co 

(O 

CM 

CM 

CM 

Nt 

CM 

cm 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

mO 
OZ 


•-<cMco>4-mNors-cocjN 


20U 


cm     rq 


I 

Q— 

ijjQ. 
O.0C 
co— > 


I 


2CU 
OJ-J 

s:  l 

QZ 


^cooutcocoir\>^r«-r"-rocoo»Ni«- <      in      o      >$•      u»      >fr 
COOO>co^r>OCOCOOOOOaOOCoCOO 

o        •        •        ■ •        ••••• 

acro>4-<j>^coro^cococoinco>r'-*c_jcor"->*«'co 
l—     ^romrnrorOrncOroc\jrvj<\jc\jcM(\JC\i(\j<Nco 


I 

LLCO 

—I 

I  — 


CO 

LU 

t—    • 

<o 

QC  LU 
LUQ 

z: 

LU 


CO 

—t 

CO 

o» 

nO 

m 

cm 

CO- 

ro 

rc\ 

in 

in 

CNJ 

—j 

>T 

ro 

— 1 

CO 

>o 

in 

o 

•4- 

00 

in 

CNJ 

co 

;\i 

>t 

i-i 

o-» 

c\j 

CO 

-o 

r- 

r- 

<f 

ON 

co 

o* 

in 

«t 

CNJ 

t~t 

o 

CO 

co 

00 

'** 

m 

m 

^r 

ro 

i\i 

rH 

-T 

«*J 

CO 

rorrifOrofflrnrgfMMMfMNMfxIMlMMMm 


soo^NOO^m^cNjr^^-h-f^-inomocoocNj 
»—  ror-^m(^i-too(\jo>t^'4'>0'-<mmcoi-i 

=>••••••••••• •• 

ocM>rtMmcMNfr'nfr>.4-invfr>frm^inmcom>t 

»—      rorororOrOfOfOrorororOfOrororororOrOfO 

ir\(D5\r\lX(OlOM(\lH^C0^QHC01DN^ 

2?:>tOs0a3'-t^o>mcoco>oo>oco'-'o>o^m 

nO^OnOvOnO-0,OvO>0»0,OvOvO>0>0>OsO>00 


m 


ex 

a: 
o 

CO 
GO 

I  LI 

a. 
s: 

o 
o 

o 

I— I 

Ot 
CO 
> 
I 

I 

I 

I 

< 
r- 
< 
Q 

O 

LU 

cc 

CO 
< 


LL 

O 

LU 

—I 
CO 

< 


>t 

CO 

o- 

O 

r- 

r-t 

r- 

<i 

o 

in 

ON 

m 

r*- 

ro 

f~t 

m 

in 

CO 

rH 

* 

O 

<f 

i—i 

CO 

ro 

o 

r^ 

h- 

in 

CNJ 

CO 

CO 

h- 

in 

rO 

o 

vO 

r-i 

CO 

Cl 

o 

in 

in 

<f 

>*• 

<J- 

co 

cO 

CO 

ro 

rNI 

CNJ 

CNJ 

CNJ 

rvj 

CM 

CNJ 

CO 

4- 

m 

m 

ro 

m 

rO 

CO 

r<l 

CO 

ro 

ro 

CO 

CO 

ro 

CO 

rO 

CO 

ro 

ro 

CO 

CO 

ao 

in 

o 

>fr 

00 

CO 

<o 

ro 

rH 

■4- 

ro 

o 

r- 

r-4 

CO 

rH 

o 

r-4 

r— 

UJ 

CNJ 

CNJ 

i— « 

in 

ou 

CNJ 

CO 

u\ 

in 

ro 

O 

CO 

r- 

o 

in 

>T 

<r 

vO 

o 

CO 

cx-» 

K 

ID  CO 

Q. 

CO 

CO 

r- 

n£> 

O 

in 

in 

in 

in 

in 

<f 

>4- 

vt 

>fr 

-t 

<!■ 

>4- 

in 

o 

1/>C0 

m 

en 

ro 

CO 

fO 

ro 

ro 

ro 

m 

ro 

ro 

ro 

CO 

ro 

rO 

ro 

rO 

ro 

ro 

co< 

LU    • 

QCO 

U.X 

co 

o 

rH 

^ 

O 

o 

CO 

0> 

rNi 

r- 

CNJ 

m 

r- 

r- 

(NJ 

o 

(NJ 

i— t 

o 

• 

rH 

o 

CO 

CO 

>3- 

■— t 

r- 

in 

in 

>r 

rH 

O 

CO 

CO 

r- 

r- 

r- 

CO 

f-^ 

ir* 

-J  CO 

H 

<Z 

a. 

vO 

in 

vf 

*r 

«*• 

ro 

ro 

ro 

ro 

ro 

CO 

(NJ 

rNj 

CNJ 

CNJ 

CNJ 

rNI 

rO 

>r 

j— »— • 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

rNj 

CNJ 

Cn! 

rvj 

rvj 

CNJ 

CNJ 

CNJ 

CNJ 

(NJ 

CNJ 

Ow 

r- 

r- 

CO 

CO 

00 

o- 

co 

CO 

r- 

r«- 

CO 

00 

r«- 

co 

ON 

CO 

00 

CO 

CO 

CO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

O 

O 

o 

< 

o. 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

CO 

co 

CO 

co 

fO 

CO 

ro 

ro 

ro 

ro 

ro 

ro 

rO 

CO 

ro 

CO 

ro 

ro 

ro 

1 

v0 

CNJ 

>r 

>0 

-t 

o 

n0 

m 

•O 

in 

>t 

(O 

rO 

CNJ 

CM 

co 

ro 

in 

vO 

CO 

co 

o> 

o 

r> 

0> 

O 

c* 

O 

0> 

ON 

ON 

CO 

ON 

CO 

CO 

CO 

o 

CO 

>*- 

CO 

CV 

CO 

go 

co 

co 

co 

co 

O 

o 

& 

o 

CO 

ON 

O* 

0> 

CO 

CO 

CO 

ON 

CO 

LU 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

(NJ 

CNJ 

CNJ 

(NJ 

CNJ 

CNJ 

CNJ 

rNI 

(NJ 

CM 

CNJ 

CNJ 

a— 

3cO 

co  en 

CO< 

rH 

h- 

r- 

*o 

si" 

rH 

O 

rH 

ro 

CNJ 

in 

CNJ 

o 

r- 

«*■ 

CNJ 

CNJ 

00 

o 

LU    • 

00 

— i 

CO 

m 

rO 

rH 

O 

o 

o> 

CO 

r- 

h- 

r- 

o 

vO 

vO 

0- 

o- 

s0 

a:o 

CNJ 

Q-X 

a. 

0> 

CO 

CO 

CO 

CO 

00 

CO 

CO 

r- 

h- 

r- 

r- 

r» 

r- 

0- 

r- 

0- 

r- 

CO 

• 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNI 

CNJ 

CNJ 

(NJ 

CNJ 

OOO 

~z 

r-»-« 

<— 

rH 

CO 

rH 

CNJ 

in 

00 

h- 

vj" 

O 

Co 

r- 

r^ 

in 

CO 

in 

vO 

rH 

ro 

>!■ 

1— 

>3- 

O 

ro 

<T> 

»o 

CNJ 

o 

r-^ 

O 

0- 

<j 

vO 

m 

•J- 

Nf 

>r 

vO 

r* 

a» 

CO 

rH 

CL 

in 

>fr 

<4" 

co 

ro 

CO 

CO 

fO 

CO 

CNJ 

CNJ 

CNJ 

(NJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

ro 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

rNj 

CNJ 

rNj 

N3 

OZ 

a 


CNJ      rc> 


in     *o 


CO      o> 


(NJ       ro 


in    no 


CO        Q\ 


205 


I 

o— 

O.OC 

WW 
I 

(— — 

ZCO 
LU-J 
21   I 

QZ 


^cu^j-r—vo^vrvor-     — '(Mr—     is-vOcmis-<m«-'o<-> 

oincorOvOincM»-«cM»0'"",coois-rococ\jrMcj 
C?        ...... •         • 

£x«->r,*-rnms0>TCvcr'tjoi^"O<-<>i'Lr)>Tc\j<«r>u 


I 

3EUJ 

ILCO 
—I 


CO 

LU 
OS  — 

t-  • 
<o 

uCUJ 
lilQ 
□_  — 
X" 
UJ 


in 

CO 

cm 

O 

CM 

co 

in 

co 

0 

cr- 

r- 

m 

cm 

O 

cr. 

>o 

in 

CM 

o> 

CSI 

co 

CM 

1— » 

o> 

r- 

>0 

-4- 

^0 

~* 

0 

cr- 


>* 

ti 

CM 

co 

O 

>o 

ro 

0 

nO 

Is- 

rO 

00 

>* 

0 

ro 

vO 

Is- 

m 

«$• 

vO 

O 

t-l 

in 

r\l 

in 


mrorOrnrr>mrOCsj(\j(\jc\jiNjrrirorri 


H-m<7>co{\jor*-.-<rs-oorr>inmr*-!-ir-'-tr»-mo 

O        .........        < .        . 

0^-ini^-»oaoc^r*-p,-NOis-Os0^mvor~oooo 

2<NjrotMvto^asmomO'-«rMrs-oo(Nmr-4>rfo 
r-m>j-rs-stis-sOrs-OvOis->ominsf     o     %o     r-     00     a* 


x> 

QC 

u 

CO 

CO 

LU 

a 

Ql 

2: 
o 
o 


a. 

CD 

> 

X 

I 
I 
I 
<r 

< 
O 

Q 
iu 
a: 

CO 

< 
LU 


LL 
O 

UJ 

-I 
CO 
< 


<r  00 
<t     <3r     a 


r\j  in 


CO 

NO 


in 


in 

vO 


CM 


(X>    in 


Is- 


sO 

cr 


r- 
m 


cm 


CO 


in 


0  00 


Ol 

h- 

Is- 

vO 

sO 

in 

in 

«t 

vt 

ro 

co 

ro 

CM 

CM 

CM 

CM 

ro 

CO 

>*• 

Is- 

ro 

CO 

CO 

co 

co 

ro 

m 

m 

ro 

co 

ro 

cr, 

CO 

CO 

co 

rO 

CO 

ro 

ro 

CO 

co 

m 

CT> 

sO 

cr. 

CM 

NT 

CM 

^f 

>4- 

>t 

Is- 

CO 

1— t 

m 

CT- 

00 

ro 

LU 

IM 

-J- 

m 

co 

cr- 

CM 

r- 

V 

in 

i-H 

r*- 

>r 

r-t 

CO 

Is- 

0 

r*- 

Is- 

O 

in 

a:— 

1— 

3  CO 

Cl 

■-i 

0 

cr. 

CO 

CO 

r- 

-0 

-O 

O 

in 

in 

in 

>1- 

>t 

in 

m 

m 

r*- 

0 

l/>CO 

-4- 

>t 

co 

co 

co 

n"> 

ro 

m 

CO 

CO 

ro 

co 

co 

CO 

CO 

ro 

CO 

CO 

vl- 

cO< 

LU    * 

a:  CD 

O-X 

-r 

r- 

CM 

ro 

>*• 

CM 

0 

-O 

CO 

CO 

CO 

<r 

0 

in 

>*• 

r- 1 

r- 

— 1 

Is- 

• 

r-< 

—1 

vO 

CM 

vO 

CM 

CO 

«4- 

cr< 

r- 

in 

CO 

1— • 

0 

CO 

0 

in 

vO 

m 

Is- 

_JCO 

H 

<z 

Q. 

Is- 

vO 

O 

in 

in 

>J" 

>4- 

m 

CO 

CO 

"0 

CO 

co 

CM 

ro 

ro 

co 

•sT 

0 

H-»-t 

«NI 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rM 

CM 

CM 

CM 

O  — 

(— 

CO 

Is- 

Is- 

vO 

in 

>t 

in 

in 

CO 

ro 

<f 

CO 

<*■ 

CM 

co 

ro 

CO 

vO 

CO 

O 

0 

0 

0 

0 

O 

O 

0 

0 

O 

O 

O 

0 

O 

0 

O 

O 

O 

O 

< 

Q. 

O 

0 

O 

O 

0 

O 

O 

O 

Q 

O 

O 

O 

0 

O 

O 

O 

O 

O 

O 

ro 

m 

co 

co 

m 

rn 

m 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

rO 

CO 

CO 

CO 

1 

•— 1 

CM 

i— 1 

CO 

r» 

«t 

0 

vO 

CO 

— 1 

1-4 

o> 

0 

CO 

Is- 

.—1 

CM 

O 

Is- 

C\J 

CM 

CM 

«-4 

.-1 

1— 1 

^H 

0 

O 

O 

O 

cr- 

cr- 

cr- 

0 

O 

O 

•—1 

i-l 

>d- 

CO 

Q. 

O 

O 

O 

O 

O 

0 

O 

0 

O 

O 

O 

0> 

o> 

cr. 

0 

O 

O 

O 

O 

LU 

ro 

co 

co 

co 

ro 

ro 

CO 

m 

CO 

co 

rO 

CM 

CM 

CM 

CM 

ro 

CO 

CO 

CO 

CL~ . 

DCO 

CO  CO 

CO< 

CNJ 

co 

in 

CO 

1— 1 

m 

in 

0 

1—t 

r-l 

-t 

vO 

CM 

Is- 

CM 

O 

r- 

CM 

CM 

LU    • 

O 

in 

1— 1 

sO 

>*• 

1-* 

CO 

0 

m 

«* 

rO 

CM 

CM 

i-4 

CM 

>fr 

»t 

O 

vO 

aCO 

cm 

o-x 

O- 

O 

o> 

0 

CO 

CO 

00 

r- 

r- 

r- 

Is- 

r- 

Is- 

Is- 

Is- 

Is- 

Is- 

Is- 

Is- 

O 

• 

m 

CnI 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

i_)cO 

»-*z 

h-»-H 

<w 

in 

.-« 

r- 1 

Is- 

CM 

-r 

cr- 

CO 

in 

cr. 

co 

nJ" 

m 

'T 

-O 

CM 

Is- 

CO 

O 

H 

CM 

CO 

>t 

CO 

in 

<— • 

h- 

CO 

CM 

0 

00 

Is- 

0 

in 

O 

O 

^H 

00 

cr 

co 

1— 1 

vO 
CM 


in  m 

CM   CM 


CM 


>J->rcOfOCOrOCMCMCMCMCMrO 
CMCMCMCMCMCMCMCMCMCMCMCM 


CO      CO 
CM      CM 


in 

CM 


•-a 

oz 


CM      CO 


in    vO 


00     cr-     o 


cm     co     >r     in     O 


00     0> 


206 


TABLE  I-IV:  RESULTS  OF  COBRA  PROBE  SURVEYS  AT  THE  EXIT 
OF  THE  HYBRID  COMPRESSOR  ROTOR 


The  table  is  in  12  sections: 
Section      Pt.  No. 


(i) 

23 

(ii) 

25 

(iii) 

27 

(iv) 

11 

(v) 

12 

(vi) 

12 

(vii) 

13 

(viii) 

18 

(ix) 

19 

(x) 

16 

(xi) 

18 

(xii) 

19 

-\ 


Run  No 


h9 


51 


52 


53 


5h 


Diffuser 


B 


B 


Diffuser  blade  angles  and  other  test  parameters  are  given  in  the 
previous  tables . 
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